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ESTIMATING THE EFFECT OF QUANTIZATION BY LEVEL 
ON THE PROCESSES IN AUTOMATIC DIGITAL SYSTEMS 


Ya. Z. Tsypkin 


(Moscow ) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 281-285, March, 1960 


Original article submitted July 28, 1959 


An estimate is given of the effect of the level-quantization interval in automatic digital systems whose con- 
tinuous portions contain both constant and variable parameters. 


INTRODUCTION 


In automatic digital systems (ADS) at least one quan- 
tity is subject to quantization by level and by time. The 
quantization process consists in fixing previously given 
discrete levels at discrete moments of time. Level-quan- 
tization interval og is determined by the accuracy, and 
time-quantization interval T by the speed of action, of 
the digital devices which belong to the automatic system. 

The effect of the time-quantization interval was 
studied earlier[1,2]. An estimate of the effect of the 
level- quantization interval was made in [3] on the basis 
of the canonical forms of equations. Due to this, the 
results of the estimation process were obtained in an un- 
necessarily complicated form. Moreover, it is not pos- 
sible to concur with all the conclusions of [3]. 

Below, we present a simple way of estimating the 
effect of the level-quantization interval which is not- 
related to the canonical forms of equations, and which 
is valid for ADS with constant and with variable param- 
eters. 


1. Equivalent Scheme of an ADS 





We consider the ADS scheme shown on Fig. 1. The 
functional schematic of the device (AD) for analog-to- 
digital conversion (Fig. 2, a) may be presented in the 
form of a series connection of a pulse element and a 
relay element with a multistep characteristic (Fig. 2,b). 

In its turn, the relay element can be presented in the 
form of a parallel connection of an amplifying element 
with gain k = 1 and a nonlinear element (Fig. 2, c) whose 
characteristic ¥ equals the difference between a linear 
and the relay characteristic. 





f Continuous zd 


AD portion 
[ap 




















Fig. 1. 


























Fig. 2. 


By using this representation of the AD converter , we 
can replace the ADS scheme by the equivalent scheme 
of a limiting sampled-data system in which the effect of 
the level quantization is taken into account by the 
stimuli applied to it via the nonlinear element with 
characteristic ¥ . In Fig. 3, a, this system is represented 


























as a sampled-data system with a linear feedback path; 
on Fig. 3, b, the feedback path fs included in the limiting 
system. 

We denote the impulse characteristic of the closed 
linear sampled-data system, which can also contain variable 
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parameters, by k,{n, m,€]}. Methods of determining 


parameters, then 


k,[n, m,e] =k,[n—m, ¢]. (1.1) 


By using the impulse characteristic, we can write 
the equation of the ADS of Fig. 3, b in the form 


natn; a]00 Dy ken, m, stm) — 


¥ (z4[m, 0))] + f{m)). 


2. Deviations Engendered by Level Quanti- 
zation 

A system which differs from an ADS by its absence of 
level quantization we call a limiting automatic sampled- 
data system (ASDS). The equation of the limiting ASDS 
can be obtained from the equation for the ADS (1.2) by 
setting in it 





¥ (/{m}) = ¥ (calm, 0}) =0. (2.1) 
Thus, for the limiting ASDS we will have 
zs[n,e®] = >) ke{n, m, e] /[m] (2.2) 


m=0 


or 
z,[n,e]= >) k(n, n—-m,e]/[n—m]. (2,3) 
m=o 


The difference between zgn, €] and z,.[n, €], 
defining the effect of level quantization, may be presented 
on the basis of (1.2) and (2.2), in the form 


z[n, e] = z4g[n, e] — 2, [n, e] = s k. [n, m, ] 
m=0 


(¥ ({{m])—¥ (z,[m, 0))}. — 


The deviation characterizes the effect of level 
quantization. The next problem consists in estimating 
this deviation. 


3. Estimate of the Maximum Deviation 





We now find an upper bound for |62{n, €]| By 
replacing, in the left member of (2.4), the absolute mag- 
nitude of the sum by the sum of the absolute magnitudes, 
we get 


[az in, l)< SF [kgln, m,e]| {| ¥ (/ fm) | + 


m=0 


+ | '¥ (glm, 0}) |}. 


(31) 
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But the difference between the linear and the relay 
k{n, m, €] were described in[4], If there are no variable characteristics or, what amounts to the same thing, be- 


tween the nonquantized and the quantized quantities , 
may not, in absolute magnitude, exceed half the quan- 
tization intervalo/2. Therefore, 


[¥ (flm)) |<, | ealm, |< 5. 2) 


By using these inequalities, and replacing the left 
member of (3.1) by the maximum value, we obtain 


max |éz[n, €]|< >) lAcIn, m, €]|¢ 


m=0 


(3.3) 


or, after replacing the variable of summation m by mm, 


max | 6z[n, €]|< >>> lk.[n, n—m,e]|o. (3-4) 


m=0 


On the basis of (1.1), for an ADS with constant param- 
eters we obtain from (3.3) and (3.4), respectively, 


n 


max | 6z[n, e]| < >) |ke[nm — m, e]|c (3.5) 
m=0 
or 
max | 6z[n, ¢] |< ir |Ae[m, e]|c. (3.6) 


m=0 


Thus, the process’ maximum deviation engendered 
by level quantization does not exceed the sum of the ab- 
solute values of the impulse characteristic of the limiting 
ASDS multiplied by the quantization interval. 

An estimate of the steady-state deviation is found 
from (3.4), (3.6) if the limit of summation n tends to @. 

For an ADS with constant parameters, the estimate 
of the steady-state deviation is the greatest possible one. 

This estimate corresponds to the maximum possible 
accumulation in the limiting ASDS for stimuli which are 
modulus-limited, and is the discrete analog of the ac- 
cumulation problem of B. V. Bulgakov. 


4. Comparison of Some Properties of the 
Limiting ASDS and the ADS 








If a step function stimulus 
1 for n>O, 
fInj=j oo 
0 for n < 0, 


is applied to the limiting ASDS, its reaction ts determined 
by its time or, more accurately, by its transient charac- 
teristic. By using (2.3), for example, we can present the 
transient characteristic h[{n,¢€] in the form 


h[n, e] = zs[n, e] = >} ko[n,n—m,e]. (41) 
m=0 








2) 


.3) 


4) 


ram- 


3.5) 





The steady-state characteristic will be 


co 
hstein, #] = ¥) k,[n, n—™m, &]. 
m=0 


(4-2) 


It is a function ofn andé . Foran ADS with cons- 
tant parameters, 


he [n, €] = h(e) = p> km, ] 
m=0 


is a function of € only. If the ASDS is designed to re- 
produce the stimuli applied to it then, in general , 


hye (n, ©] =h(e)<1. 


(4.3) 


(404) 


We assume that the impulse characteristic of the 
limiting ASDS, k,[n, n-m, €], for givenn andé€ and 
arbitrary m , is nonnegative, i.e., 


ko [n, n — m, e] >0, (4.5) 
so that we shall have, from (3.4) and (4.1), 
n 
max |6z[n, e] |< be k,[n, 
m=0 (4.6) 


n—m,e]¢o=h[n, e]c, 

Thus, in this case the maximum deviation at the 
given moment of time which is engendered by the level 
quantization does not exceed the quantity equal to the 
product of the quantization interval by the value of the 
characteristic at that same moment of time. 

The estimate of the steady-state deviation will have 
the form 


co 
max |5z,,.[”, €] |< > k,[n, n— m, 4.1) 
m=0 


&]o = hye [n, €] , 


Let us consider in somewhat more detail an ADS with 
constant parameters. For it,the [( 4.5)] that k,f{m, €] = 0 is 
equivalent to the condition that h[n,€] be monotonic 
for € = const. By taking into account that h[n, € ]= hgte- 
fn,e} ©h (e€), we find that 


max |8z [n, e]| Ch[n, e]o < hye [n, ] ¢ = A(s) 2 (4-8) 
or, by taking (4.4) into account, 


max | 6z[n, e]| <ce. (409) 


Thus, for a monotonic characteristic, the deviation due 
to level quantization does not exceed the quantization 
interval. In this case, and only in this case, we can 
deduce the stability of the ADS from the stability of the 
limiting ASDS. If the characteristic is not monotonic, 


i.e., the impulse characteristic changes sign, the sum of 
the absolute values of the impulse characteristic may 
significantly exceed unity, and instead of (4.6) we shall 
therefore have 


max |z(n, ¢]| Cac. (4-10) 


where a> 1. Under these conditions, it is possible 

to have, in an ADS with constant parameters, an auto- 
oscillatory mode with an “amplitude” not exceeding ao. 
The greater the oscillation of the characteristic h[n,€ ] 
for € = const, i.e., the greater the oscillation of the 
limiting ASDS, the greater will be the “amplitude” of 
the autooscillatory modes possible in the ADS. 


Thus, in spite of the assertion made in [3], it is, 
strictly speaking, impossible in this case to conclude that 
the ADS is stable from the fact that the limiting ASDS 
is stable. 


The result just derived does not depend on the size of 
the quantization interval. However, as the quantization 
interval o decreases, i.e., as the number of places in the 
digital devices increases, max |Sz[n, €]| will increase 
and, despite the presence of the possible autooscillations, 
their amplitude will be so small that the processes in the 
ADS and the ASDS will differ less and less from each 
other. 


The estimates given above allow one to determine 
the maximum quantization interval for which the dif- 
ference of the processes in the ADS and the ASDS will not 
exceed a given magnitude. 
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DETERMINATION OF THE OPTIMAL IMPULSIVE RESPONSE FUNCTION 





IN THE PRESENCE OF INTERNAL NOISE 


P, S. Matveev 


(Moscow ) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 286-292, March, 1960 


Original article submitted May 8 , 1959 


The problem of determining the optimal transfer function [1, 2] is generalized to the cases when the input 
stimuli are applied ton different elements of aservosystem and to systems with variable parameters. 


In [1, 2] a solution was given to the problem of 
determining the optimal impulsive response function for 
cases when input stimuli are applied to two different 
elements of a servesystem. The aim of the present work 
is to generalize the results obtained in [1, 2] to the case 
when the stimuli are applied ton elements of a servo- 
system, and to systems with variable parameters. 


1. Determination of the Impulsive Response 
Function in the Class of Systems with 
Constant Parameters 











We assume that the servosystem has n elements to 
which stimuli are applied, whereby a controlling signal 


y (t) = g(t) +mi(t), (1) 


is applied to one basic element, where g (t) is a given 
function of time, with 


























In addition, noise n(t), in the form of a stationary 
random function, is superimposed on the controlling signal. 
To the remaining n— 1 elements of the servosystem (Fig, 
1), stimuli are applied in the form of stationary random 
functions uy (t), Ug(t), «+ + » Up-y(t) We assume for 
simplicity that m(t), m(t), ug(t), Ug(t), « + +, Un=g(t) have 
zero mean values and are mutually uncorrelated. We also 
assume that the transfer functions W,(p), W2(p),-+-+, 
Wn-1 (p) [or the corresponding impulsive response func- 
tions b,(t), b(t), « « +, bn= (t)] are known. What is not 
known is the transfer function W{p) (impulsive response 
function) of the correcting device. We initially determine 
the optimal impulsive response function k(t) of the closed 
system from which, by using well-known methods [2, 3], 
we can easily determine the transfer function Wp) of 
the correcting device . The problem may be formulated 
in the following way; from given correlation functions 
Rin (T)s Rn (T)> Rut)» RugT «+++ Rup, (7) [spectral 
densities Sm (w), S,(), S,,(¥). Sy(H)» coe, Sun-1 (w)], 
given duration T of the transient response, given re- 
production operator H(p), error coefficients C; and trans~ 
fer functions W,(p), Wop), «+ + » Wy-3 (Pp), to find the 
impulsive response function k(t) such that the mean- 
square error will be a minimum. 


pst) 
ho 


u,{t) 








z(t) 





W,/p) 


























r 
g(t)=>ig,t', (2) 
i=o 
and m (t) is a stationary random function. 
glt}+mitt)+ntt, we 
+m(t)+n(t) 
Yip) Wp) 
Fig. 1. 


It is easily shown that the scheme shown on Fig, 2, for n = 3, is equivalent to the scheme of Fig. 1. This 
scheme will be used for the subsequent argument. Then [4], the system's output quantity can be written in 


the following way : 


T 
z(t) = (Ig (¢—2) + m(t—1) + n(t—x)) k(x) de — 


0 
T 


—(k(ar | b, (c) do i 
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(3) 


Uy (t —t —o —p) bg (u) dp — 























4 © co co 


—lk(e)de \ u(t —t — cc) b, (0) do + \ b, (x) dx ( u, (t —t — 0) b, (0) de 


0 —oo -<o —co 


oe ( Ug (t — t) bg (t) dt, 
where 


k(t) = N(t) + s E; 8 (t) + 5 sit—T) O<t<7). (4) 


jo jo 


Let the ideal system be required to reproduce at its output the signal 












































h(t) = H,(p) g(t) + H(p)m(t), (°) 
where 
nal. "A, ro; 
H(p)= Di 7 H,(p)=H(r)—- SP ie”. (6) 
_ j=o 
al W(B) Mo (t)+ Mofo) (t) 
(t)+m{t) +n(t) mt zit) 
. 25-1 e(p) |->| (0) [>] Wolo) | 7-88 

ine Wp) Wal) u(t) + Molo) up) 
2s Fig. 2 
‘ Here, for convenience, we have introduced the operator Hg(p) which simultaneously takes into account the 
| operation of reproducing g(t) and the errors with which it is reproduced. 
“4 Then the reproduction error can be written in the form 


s(t) =h(t)— x(t). 


If it now be required that the mean error of reproduction equal zero, then [2, 4, 5] we are led to limitations 
of k(t) of the form 


T 
thy = (— 1)"| 2" (2) de (v=0, 41, 2,...,7) (7) 
0 
or 
T 7 
H, =(— 1)'\<'k (x) de, —C;=(— 1) |<! k(t) de (j sei = 0, 1,2,..., 7), (8) 
0 0 


If (8) holds, the reproduction error can be written in the following way: 


7 


e(t) = \ m (t — 2) x(x) de —| k(t) dt[m(t —t) +n(t—*)— 





mt b, ()dt ( ni(t—+—0— wh wda— | us (t —t— 0) by (2) do | — (9) 


oa b, (t)dt \ ui(t—+—o)by(o)do— | Ug (¢t — t) by (t) de, 


—0oo —oo 











x(t) = = \ H (jw) e do, 


—oo 


Further, by using the same arguments as in[2, 4, 5, 6], we obtain the necessary and sufficient condition for the 
minimum average value of the square of (9), when (7) holds, in the form of an integral equation 


(Rn (c — 6) + RA, (t — 8) + ( by (s)de ( bg (4) du \ by (A) dh 


x | Ra (et+o+¢—9—2— x) ba(n) dn + 


rs ( b (0) de | Ru, (¢ + 9 —0 — 2) by (4) dd] k (0) a = 
—0o —0oo (10) 
- ( Rim (t — 8) x (8) d0 + ( by(2)de | bs (u) du ( 6, (0) 40 | Ry, (t+0+ 


) co r 
+p—0—2)d, a+ | dy (o)de | Ru(e+e— 9) ba (0) d+ Sx, 
—oo —oo i=0 
where the y; are Lagrange multipliers. 


We now introduce the notation 
co 


Ry, (t — 9) = \ by(a)de | ba w)dy ( by (A) dh x 


—oo 


(11) 


—oo 


x \ Ry, («+0 +p —0 —2— 4) by (x) dn, 


Ri, (« — 9) = \ b,(2) do \ Ry, («-+ ¢ —8—2)b, (A) dd. (12) 


With the use of the notation of (11) and (12), (10) can be presented in the following form: 


T 
| [Rm (« — 8) + Rn(t — 9) + Ru, (« — 8) + Ri, (« — 9)] k (8) dd = 


4 r 
= | Rm(s—0)x(0)d0 + Ri, (2) + Ral) + Hut O<t<7). (13) 
—oo i=0 
Equation(18) does not differ in form from the integral equations obtained in [2, 4, 5, 6] and, therefore, the 
same methods of solution may be used for the determination of k(t) Obviously , integral equation (13) is easily 
generalized to the case of n inputs. For this, integral equation (13) can be written in the form 
T 
| [Bim (t — 8) + Rn (t— 9) + Re, (t— 9) + Ra, (2 — 9) + 


422+ Rayy yy (#— 9] & (0) a8 = \ Rin (* — 9) x (0) d+ oi 
+ Ri, (2) + Ru,(*) +... + Rag» (t—-9+ Hyt O<t<7). 
i=0 


If the correlation function R,,(T) + Ry (T)+R* (7) tees AR 1) (7) corresponds to a rational-fraction 
spectral density . 
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ao + ayo? +...+ 4,0” M (wo) M* 
Whe) os Kk? _ _-M (w) M*() 
(o) bo + byw? +... + b,07 L(w) L*(w) ’ (18) 
then, by using the same arguments as in[5], we can shown that the impulsive response function is defined by the 
expression * 
e r 2k q . q . 
k(t) = Sy At’ + SY) Bie! + S £89 (y+ HY DEM —T) + 
i=0 i=1 j=0 j=0 
co 
+ L(p) L*(p)M-*(p) M™*(p)| \ Rim (t — 2) x(2) de + 
—oo 
+ Ri, (t) + Rig (t) +... + Raga (| (0<t<7), (16) 
where q= I--k-1 and L{p) and M(p) equal L (w) and M(w) for p = jw, which are easily determined from the spectral 
density of (15). 
The unknowns Aj, Bi, Ej. Dj in (16) are determined as in [5]. Further, by the usual methods [3], the 
) transfer function of the closed system was determined and later, also, the transfer function of the correct- 
ing device. 
2 Determination of the Impulsive Response Function in the Class of Systems with 
Variable Parameters 
We assume that the system (Fig. 3) is acted upon by stationary random noise n(t) and signal y(t) = g(t) +mt) 
consisting of the given function of time 
= 1» gt 
g(t) = >} 8:8 ’ (17) 
i=0 
where the g; are unknown coefficients and the t! are known functions of time, plus the stationary random function 
l) m(t)} To the system's second element there is applied a disturbing stimulus in the form of the stationary random 
function u(t) 
(t)+n(t) ies 
t)+m(t) +n rt 
g(t) va wip) bet! 
2) 
Fig. 3. 
The problem solved below may be formulated as follows; from given correlation functions R,,(7), Rn (7), and 
Ry(T ) (spectral densities) and given object transfer function W,(p), to find the servosystem's impulsive response 
function such that the system will provide an exact reproduction of g(t) at each moment of time t > 0 and a 
minimum dispersion of the random process at the system's output at each moment of time. 
. smttntt)-Wyput) youl) 
+n{t)-Wfpjult r(t)=2,(t)+W,(p) ult, 
att mp) ofa [a ={t)+Mipult 
Fig. 4. 
It may be shown that the scheme of Fig, 4 is equivalent to that shown on Fig, 3. Henceforth, the scheme of 
Fig. 4 will be used. 
The reproduction error of the system (Fig. 3) can be given in the form 
e(t) = g(t) +m (t)— x(t) = [1+ Wily) Walp, OF le) + mH) — OMI— ae 
4) —[i + W, (p) W, (P, ty" WwW, (p) W, (Pp, t) n(t), 
where Q(t) = W,(p) u(t). 


The reproduction error of the system given by Fig. 4 may be written in the form 


_— 


* An analogous generalization can also be made for the case when g(t) is a harmonic function [4] . 








&, (t) = g(t) + m(t)—Q(t)—2,(t) = 
=[1+ W,(p) Walp, Ole) +md—Q@()] — 
—{1+ W,(p) W2(p, £)1-* Wy (p) W2(p, t)n(t). (19) 
By comparing (18) and (19), we find that € (t) = € x(t) and, consequently, just as in[4, 7], the reproduction 
error can be written in the form 


e(t) =a (t) + m()— Q(t) —2(0) = 8) + mM — OW — 
— | te) + mH)—Q@) + nk, Hae «>>. 


ki (20) 
The requirement that there be no bias, i.e., zero dynamic error, leads to r + 1 conditions of the form 
tx \ vk (t, t) dr. (21) 
0 
If (21) holds, the expression for the error has the form 
t 
e (t) = m(t)—Q(t)— | [m(z)— Q(z) + ne) k(t, tar. (22) 


—oo 


It may be shown [7] that the necessary and sufficient condition for a minimum mean value of the square of (22) 
amounts to this-that the impulsive response function k(t, T) satisfy the integral equation 


| (Rm (¢—9) + Rn (s— 8) + Ro(t— 8) k(t, 8) a 


ne ‘ (23) 
= Sut + Rn(l) + Rol), 
where +: 
Ro(« —8) = | b(0)do | Ru(e+o—8—2)bA)a, (24) 


and bfr ) is the impulsive response function corresponding to W,(p). 
For the class of stationary processes for which spectral density corresponding to Rm(r + R,, (T) + Ra(T) can be 
given in the form of (15), the solution of (23) can be written [7] in the following form; 


k(t, =F Be + 3 Aloe + BE 8) (t — x) 
f=} 
+ L(p) L’(p) M*(p) M** (p) [Rm(t) + Ro(t)] (>>) (25) 


The rest of the argument is the same as in Section1 of this paper. 
In conclusion, we consider two simple ew 


Example 1, Let 
2aM* k 
Su, (w) = N?, Su, (w) ~o+at ’ W, (p) = sts =’ We wrx 
H (p) = H,(p) =1, g(t) =80 + fit. 
Then [3], 
k2 2 NB 





Su (0) = 1 Wa (io) 1 Wa Gio) Su, (0) = ape * 


2ak2 M? 


re =|W.(i 25 = 
"MAGA ged tein (w? +- 22) (w? + a?) 








i. k2 -(2a.M? +- k® N®) w? + 2ac? M2 + 22k? N?2 
S — " s* = 2 1 1 1 a 
5) = Su, (o) + Su, (a) | (w* a?) (w* + 22) | 


w? + 2 














19) 


20) 


21) 


22) 


2) 


23) 


24) 











and 


On the basis of (16), the impulsive response function will have the form 





k (t) = Ag+ Ayt + Bye™! + Bge™ + £58 (t) + £3’ (t) + Dot (t—T) + 
+ Dy’ (t—T) + L(p) L’ (p) M~ (p) M* (p) [Ri (t) + RY, 1, 


where 


Example 2. Let 


Then, 


M (p) M*(p) = k2 [— (2aM? + k2 N*) p? + 2ax2M +. 2a%k2N2}], 
L (p) L* (p) = (x3 — p*) (a* — p*) (t? — p*) 


Aq and Az, are the roots of the equation M(\) M*(A)=0. 
The expressions for R* u,(t) and R® y(t) are determined from (11) and (12). 


g(t) = go + git, S, (@) = N*, Wi(p) = 


S., (w) = | Wi (jo) |? S,, (w) = 


On the basis of (25), we have 

k (t, t) = Ao (t) + Ai (t) t+ Eo (t)® (¢ —t) + Do (t) 8 (¢—*) 
+ L(p) L*(p) M-* (p>M*~ (p) Ri, (0), 

where M (p) M*(p) = ki N*, L(p) L* (p) = ** — p*. 
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ON THE AUTOOSCILLATIONS OF A DYNAMIC TABLE PLATFORM 





Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 293-300, March 1960 


Original article submitted December 1, 1959 


Determined in this paper are several properties of the motion of the platform of a single-axis dynamic table 
with an asynchronous motor with a tubular rotor, as well as the conditions under which autooscillations will arise when 


there is dry friction in the bearings. 


In the design and testing of autopilots, one uses tables, 
that is, dynamic models comprised of a device with a 
platform on which the autopilot’s sensitive element is 
placed. The platform turns about its center of mass in 
accordance with the aircraft's law of motion. 


Tables may be divided into two classes. 

a) Tables whose drives work like a servosystem at 
whose input a signal from an electronic computer is 
applied. The errors of such tables, given a small coef- 
ficient of friction, are basically determined by the inertia 
of ie drive and the platform . 

b) Tables which operate on the principle of dynamic 
similarity. The output element of such a table is the 
platform on which torques act which are proportional to 
the aerodynamic torques acting on the aircraft. The 
coefficient of similarity equals the ratio of the inertial 
torques relative to the corresponding axes of the platform 
(with the autopilot's sensitive element mounted on it) and 
the aircraft . In this case, the platform's moment of 
inertia is, in some scale, the aircraft's moment of inertia. 
The table's viscous friction can be taken into account by 
the proper decrease of the aerodynamic damping. In 
this case, the table's errors are principally determined 
by the presence of dry friction, which may introduce 
significant distortions in the simulation of certain modes, 
for example, motion at high altitudes where the damping 
is small. 


Therefore, the study of the effect of friction and of 
other factors causing distortions in dynamic similarity, 
including the autooscillations to be studied here, is of 
practical interest, the more so since tuning of the coef- 
ficient of static stability a, is carried out by measuring 
the frequency of the platform's autooscillations. 


A single-axis mechanical dynamic table was used 
by Oppelt [1]. The theory of such tables was developed 
by him. 


In the present paper, the nonlinear problem is solved 
by a method similar to that of A. A. Andronov. 

The single-axis dynamic table is represented by a 
mathematical two-dimensional model leading, in one 
case, to Milne's equation and, in the other to Liénard's 
equation. 
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Fig. 1. 


The table (Fig. 1) consists of a platform turned by 
an asynchronous motor witha tubular rotor whose transfer 
function (if the motor is considered as a linear link) fs 
kmo/P(1 + pPTmo), where kro and Tmo are the motor's 
gain and time constant, equal to the aircraft's time 
constant, and p is the Laplace operator. 


On the motor's shaft are fitted the slider of potentio- 
meter # and a tachometer generator with gain k,,- The 
signals are applied to the motor’s control winding via 
dividers a, and ay and a summing device. However, in 
actuality the motor has a nonlinear characteristic. In 
the absence oi any signal to the control winding, with a 
sufficiently weak impulse about its axis, the rotor slowly 
turns and stops (due to the dry friction), and for a suf- 
ficiently strong impulse, the rotor begins to rotate more 
rapidly. For certain values of ay and a, after a shock 
(an impulse about the rotor's axis) there is set up an 
autooscillation whose frequency can be varied by varying 
ae 

For investigating the autooscillation, one must find 
the nonlinear motor equation. 


The adjustment of the table (Fig. 1) under the system 
to be studied consists of a sequential regulation of the 
table parameters to quantities defined by the equations 
of motion of the aircraft relative to the axis under 
consideration. 


Let us consider, in the Liénard plane {¥,w} , where 
w =®, the phase trajectory and its mathematical ex- 
pressions, which correspond to the nontrivial steady-state 
platform modes. It is known [2] that, for a two- 
dimensional system, such modes can only be periodic. 











lo- 
he 





where 
x — 
state 


We denote by I the total moment of inertia of the 
platform and load with respect to the motor axis, by 
Myo the rotational torque applied to the motor's rotor 
and by Mg, (w) the dry frictional torque in the bearings 
and platform reducer. By using Newton's second law, we 
can write the platform's equation of motion in the form 


Iw = M,, +My. (1) 


If we express the dry frictional torque in terms of the 
platform's weight G and the constant frictional coefficient 
with respect to the motor's axis k >0, we get 


—kG fo w>O0, 
M;, () = Miro for eo= 0, (2) 
kG fr wo< 0. 


By applying asignal depending on w to the motor's 
control winding, we can obtain the artificial damping 
torque ayw, where a, is the damping coefficient. A torque 
depending on w also arises as the result of the voltage 
induced in the control winding when the rotor moves, 
from viscous friction in the bearing , from air friction of 
the rotating parts, etc. 

By denoting this torque by f(w)w, and by ignoring 
the lag in the motor's control winding, we write the 
expression for the rotational torque ; 


1.9 = [41 + 9 (@)] © + asd, (3) 


where a,;>0 and a, >0, and the function ¢ (w) is even, 
bounded and continuous, together with its first two 
derivatives, in each finite segment {wy, we} .« 

It has been experimentally established that there 
exists an angular platform velocity wy such that 


9 (w) < 0 for |wl <p, p(w) = 


= 0 for |@ |= @p, o(#) >0 for |@| > @p. (4) 


Moreover, ¢ (w) does not decrease for w > 0 and 
does not increase for w < 0. 


ee fox. 
mAh fad 
Th 


sry a 
(0) 
Fig. 2. 
By adding and subtracting ¢ (0) in (3), and by 


denoting [(Fig. 2)] ®(w) = y (w)~ ¢ (0), Ay = ay + 
+ 9(0), we obtain 


M_, =(Ai+ D(@)} @ + agh. (5) 


It is essential to note that the function © (w) is even 
and positive-definite, i.e., & (Ww) > 0 for w > 0 and 
®(0)=0. At the same time, Ay> 0 for | a,| >| 9 (0) 
Ay< 0 for | as! < | ¢(0)|and Ay> ¢ (0) since a, > 0. 


Ply, pie) 











> 





By substituting (5) in (1), we get 
Iw + [A, + D(o)jo + a, = Mg (). (6) 
The investigation “4 (6) breaks down into two cases; 
1) Ay2 0, ie., a, 2 -— (0), (a,> 0); 
2) A; < 0, ie., ay <-— (0). 
Case l, (A,20). If, with this, Mp (w) = 0, then 
(6) reduces to a form considered by Milne [3] and Sansone 


[4]. 
#0) 


Fig. 3. 








We now consider the phase trajectory for Me, (w)=0. 
If D = Ay —4Ia, = 0, then the solution » = (t), 
[ ¥ (0) #0] cannot have more than one Zero in the time 
interval t[{0, oo}, which corresponds to a stable node 


(Fig, 3). 
w 
0) . ad ? 
Fig. 4 
In the opposite case for D < 0, all solutions oscillate, 
and the phase trajectory corresponds to a stable focus 
(Fig. 4) With this, if we denote neighboring zeros of 


¥ (t) by t, amdts4, , ts< t,,,, then the Milne-Sansone 
inequality holds: 








oni eee ; 


tots y 


where v? = -D /41*> 0 and Q» is an arc between 0 and 
m /2 such that tan a =vI/A;- 

Consequently, if Ty = 2 1 /y is the period of the 
platform's free oscillations for @(w)= 0, and T =t,,, 
“ts. is the period of free oscillations for 6(w) # 0 
then, since ty, ~ts-, = Atsy, ts), we obtain the 
ine quality 

r<7,(i—%). (7) 

For example, for ay = - g{0) (i-e., for Ay = 0), 

a, = 1/2, and the period T satisfies the obvious re- 
lationship T < T», since the case T = Ty, occurs only 
for @(w)=0. 


We now consider the phase trajectory corresponding 
to (6) for Mg, # 0. Motion is possible if the angle of 
the platform's initial position »(0)[ » = 0] Mes outside 
the angle of stagnant dry friction — the segment 
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{— ver ver }, where p ¢, = KG/ag , that is, if 
|¥(0)|> % ~~ In this case, for constructing the trajec- 
tory we can use (6), in which the point {— w¢,, 0} is 
chosen for the origin of coordinates for the upper half- 
plane (since for w > 0, = ¥¢, = const) and, for the lower 
half-plane, the point {yp , 0} is chosen (since for 

w < 0, wf = const). 





$, | GJ * 


Fig. 5. 


p40) 





Then, on the basis of [3] and [4], in the case D = 0 
the trajectory (Fig. 5) will consist of not more than two 


pieces, the second of which abuts the segment (— vf,, W¢,). 


The initial conditions will be the point {y(0)+ v-, 0} 
for the upper plece and the point {y (t,)}- vr, 0 } for 
the lower piece, where t,; is the moment of time cor- 
responding to the intersection of the » axis by the 
trajectory. For D < 0, the trajectory (Fig. 6) will 
consist of a finite number of pieces, the last one of 
which abuts the segment {- ¥¢,. ¥,}. 





dy 


Fig. 6. 

For each piece of the trajectory in the case under 
consideration (Mg, (w) # 0), (7) is valid. Consequently, 
the dry friction does not change the period T of free 
oscillations, but increases the decrement of damping in 
such fashion that the radius vector p (p? = y? + uw). 
during one period T is additionally decreased by the 
amount 2%, Thus, if the radius vector is decreased by 
SP = 0 for Mfr ( #)= 0, then for Mg, (W) # 0, it is 
sa by Spme. 9 0 =? tor Pred “Ape, = 0 
“206 
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From the last relationship (for arbitrary D) and from 
the Archimedean axiom follows the finitude of the 
number of previously mentioned arcs, the last of which 
abuts the segment { -¥¢,, %¢,} 

Case 2. (A;< 0). If we introduce into (6) the 
function 
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F (w) ==1A, + D(o)]o, 


where p is the angular frequency of the undamped free 
oscillations, u* = a,/I, we get 


wo + pF (w) + p%) = ; Mg (). (8) 
For Ms, = 0 we obtain, by differentiating (8), 
o + pf(o)o + po = 0, (9) 


where OF 


1 ’ 
{(o) =F Jo =n [a, + —() + 9’ (@) @). 
For sufficiently small w, w y’ (w) *Oand 
1 
{() = {la + 9 ()L, ” 


but with our condition that — W9< W < We, ¥ (w)< 0. 








\ Fw) p(w) Fier) spur) 
en 8 a Re 
cA 
#10) 
Fig. 7. 


Since A; = a,+¢9(0)< 0D ,and ¢¥(w) fs contin- 
uous,then, for neighboring values of w, a; +¢(w) < 0 
and, consequently, f(w)< 0 for some w (Fig. 7). 

Since, by our assumptions, g(w) and g‘(w) are 
continuous and bounded, f (w) is continuous and satisfies 
the Lipschitz condition (guaranteeing existence and uni- 
queness) and, moreover since 9(w) satisfies (4), then f(wJ, 
for increasing w, intersects the axis of abscissas at some 
point 6 and, for decreasingw, at the point -5, Le., 


{(@) <0 for —8<o<4, 
f(o)=0 for  w=', o=—d a) 
f(w)>0 for —S >>. 


Consequently, the function satisfies the Sansone- 
Massera condition [4, 6], asserting that to (9) there 
corresponds a unique nontrivial mode, stable according 
to the Lyapunov criterion — a periodic solution — a 
cycle whose region of attraction [7] coincides with the 
entire { w, &} plane with the exception of the origin 
of coordinates. 

From whence it follows [5] that the homogeneous 
equation cofresponding to (8): 


@ + pF (w) + uw = 0 (12) 


also has a unique nontrivial mode — a stable cycle. The 
phase trajectory of (10) contains a cycle which intersects 
an arbitrary half-line —the ray OA — in the unique point 
B (Fig. 8). 
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Fig. 8. 


Let Mj be an interior point of the segment OB of ray 
OA. Then the trajectory passing through M, at time 
= &= 0 intersects OB in the point Ny at time t'> 0, where 
OR, 20M}. The last inequality follows from Bendixon's 
theorem [8]. Analogously , if Me is an exterior point 
of ray OA, i.e., OM, > OB, then the trajectory defines 
another point Ne, where OM, > ON, . 

We consider the phase trajectory of (8) for a suf- 
ficiently small coefficient k of dry friction. We assume 
that the stagnant angle ~segment {—ye,,y5,} — lies iriside 
the region bounded by the arcs of trajectory MjNj and by 
the segment MiN; + As in the first case, each trajectory 
consists of a segment , or segments, lying in the upper 
(w > 0) and lower (w < 0) half-planes. The equations 
of these segments are led to the form of (10) by changing 
the variable ¥ — translating the origin of coordinates to 
the point { ¥,,0} for the transition from the upper to 
the lower half-plane, and to the point {0,-¥} for 
the converse transition. We assume that the friction angle 
2, is so small that the interior MjN, and exterior M, 
trajectories are so little deformed, that the inequalities» 
OM, < ON; and OM, > ON, retain their validity. 

Since the continuous dependency of the trajectory on 
the initial conditions is retained, each point of the closed 
segment MjMg can be put into correspondence with a point 
of segment NjNe , lying strictly interior to MjM,- On the 
basis of the Bolya -Brauer theorem [9], with such a corr- 
espondence there is at least one point of the segment 
MM, which corresponds to itself. It defines the periodic 
mode described by (9) for sufficiently small k. However, 
it is impossible to assert in this case that such a modew 
be unique. This is borne out by the following example: 
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Let there exist an interval of variation of the angular 
velocity {- %, %} such that ®(w)=0 for|w|<| |, 
We then consider the trajectory defined by the equation 


Iw + Aw + agp = Ms, (w) (13) 


and passing at time t = &% = 0 through the point 

{{U,+ ¥¢,), 0} (Fig 9), where Uy>0 and vf = kG/ ag. 
At time ty = Ty/2> 0, the trajectory intersects the 

axis of abscissas in the point {U, + ¥¢, 0} . Here, 

Ty = 27 /v ds the period of the platform's free oscillations. 
Solving (13), we get Up, + 2, = Uy exp (~Ay7/ 2Iv), 

where A,=a,+9(0)<0, , _ V—D and D = Aj — 


nA 
- 4la, <0. whence ,U, = 0,04 (— 75) - — 2H... 





Motion is possible for U, >0, L.e., for 
/ 


U, exp i ra) > 2b¢ 





or ats 
U,> 2, exp (=) . 

Since the laws of motion in the upper (w > 0) and lower 
( w <0) half-planes are identical with respect to their own 
centers of rotation, My (for w>0) and M, (for w< 0), it is 
then necessary and sufficient that, for the representative 
point of the trajectory moving in the lower half-plane 
and passing through the point { U, + ¥,,, 0} at time 
t; = Ty/2 to fall on segment NyMj, the following two 
inequalities hold: 0 < U, < U;. 

From the inequality 


a U; =U, exp (— 7) — 2H, = 


s [Us exp (— Fay — 244, |oxp(— z=) 
— 2. <U, 


follows 
zy fr 


exp (-75) —1 


From the inequality U, > 0 follows 





U,< 


hr or 
2p [: +ex0( 75) | 
nA, . 
#p (— =p) 


{ a. 
+ exp (F55 de ye. 





U,> 




















exp (— =a) exp i 75) —i1 (14) 











mah 
nA, 
exp * 75) —1 


point falls on the segment { - fr ¥frpt+ and motion 
ceases. 





For U/, > 


the representative 





r+ (45) 
(75) 


point falls on the ray defined by the inequality 
v < - (Uy + ve. ); the system will oscillate. 

From the single-valuedness and continuous dependency 
on the initial conditions of the solution of (13) there flows 
the existence of a nontrivial steady-state mode-cycle with 
— frequency f = v/2m and with oscillation spread of 

= I¥min | +1¥ max!» where 


For U,< 





the representative 

















2 
I tnin! = U1 + oe 3 7 + Ver» 
=a 
2 
| nay |= Ua + ty =——*# | — — & - 
‘+x (75) 


From whence, 
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This mode is unstable in the strong sense, according 
to the Lyapunov criterion , i.e., both on the interior and 
the exterior sides of the region, the orbits of this mode 
are limited. The topological picture of the interior 
portion of the domain is a set of trajectories consisting of 
a finite number of rays, the last one of which abuts the 
resting segment MyM The exterior side consists of 
trajectories swept out with the cycle's orbits. 

In the example considered, it is natural to assume that, 
for sufficiently largevalues of the distance R of the repres- 
entative point from the origin of coordinates (R*= y*+ w”), 
there will appear additional dissipative torques applied 
to the table's platform. The mathematical form of these 
torques can be found if we start from the assumption that 
coefficient A; of viscous friction in (13) varies as a func- 


tion of R according to the law 

? a,+9(0)<0 fr R<LV+e, 
, = 
E>0 for R>>¥+e, 


where € and E are certain fixed positive numbers. 








We consider two trajectories (Fig. 10), the first of 
which passes through point Sy,S,N, =€, at time t= % =0 
and the second through point S, for decreasing time 
parameter t, l.e., fort< 0. These trajectories define two 
points S; and Sj, where point S$ Mes to the left of point 
S;, and point 5, lies in the interval S,N; . 





Fig. 10. 


Based on what has been said, we conclude that there 
exists a periodic mode which lies in the hashed portion 
of the figure. Since all the R< ¥/2+ € define points 
lying on trajectories swept out with unstable orbits, and 
since all R>W¥ + € define points corresponding to damped 
motion (A; > 0), then R = ¥ + € defines a point lying on 
the trajectory corresponding to a stable cycle with spread 
¥,; = ¥+ 2€ and frequency f. 

The author wishes to thank B. V. Shirokorad for his 
aid in the preparation of this paper. 
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PHASE-LOCK AUTOMATIC FREQUENCY CONTROL OPERATION 


IN THE PRESENCE OF NOISE 
V. I. Tikhonov 
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Original article submitted May 20, 1959 


It is established, by means of the two-dimensional Fokker-Planck equation, that due to external and internal 
fluctuations the mean frequencies of the fundamental and the synchronized generators in a standard phase-lock 
automatic frequency control scheme do not equal each other. An approximate method is suggested for computing 
the mean difference of the generators’ frequencies and its dispersion. 


INTRODUCTION 

In [1] there was considered the effect of noise on 
the operation of a simplified circuit of a phase-lock 
automatic frequency control scheme (a phase-locked 
oscillator, henceforth denoted by PLO) in which the low- 
frequency filter was considered to be ideal. In this work, 
the circuit was described by a nonlinear first-order 
differential equation. 

In actual PLO circuits, one may use, as a low- 
frequency filter, either a low-frequency amplifier or a 
series-connected RC network. The voltage from this 
latter filter is applied to a reactive tube. 

If it is possible to approximate the low-frequency 
amplifier's frequency characteristic by the frequency 
characteristic of the RC network, then the operation of the 
PLO scheme will be described by the two equations [1] 


! = A,—su +4, (1) 
RCu +u= 5 Ag t(Ad + E,) sing — £,cos 9]. (2) 
We now adopt the following nomenclature; 
= @,—,, b = 6, — 6,, A, = WJ, — Wo; (3) 


u is the transformation coefficient of the phase detector, 
3 is the reactive tube's transconductance, w 7, 1s the mean 
frequency of the locked (synchronized) generator for 
u(t)= 0. With this, the expressions for the oscillations 
of the fundamental and the locked generators take the 
form 

u(t) = A, cos DM, = A; Cos (wot + 6,), 

(4) 
uz (t) = A,sin DM, = A, sin (wt + 0,), 


where Ay and A, are constant amplitudes , wy and w are 
the mean frequencies, 6 (t) and @,(t) are the random 
phases of the generators’ oscillations, resuiting from the 
noise of the generators’ tubes and the impedance losses 
of their elements. The normal fluctuation noise n(t), 
whose spectral density is symmetric about frequency Wp, 
is described in the form of a quasi-harmonic oscillation 
n(t) = E(t) cos (wot + 9), Ey (t) = E(t) cos (6 — 9), 
E, (t) = E (0) sin (0 — 9,). (5) 


By differentiating (1) with respect to time, and then 
substituting the expressions for i from (2) and for u from 
(1) , we obtain 


gt+ap+adsing=ab,+&(t, ), (6) 
where 


i 
a= po: A= 7 ps,A,, E(t, 9) = 


=o+ap+ a (Ey cos p — E ,sing). 
Nonlinear second-order differential equation (6), with 
a specific random function in the right member, describes 
the operation of a PLO scheme when one has taken into 
account the natural instability of the generators’ 
frequencies and the presence of fluctuation noise in the 


fundamental generator's channel. 
It should be mentioned that the investigation of (6) 


is of interest apart from its connection with the concrete 
problem being considered here. Such equations describe 
many physical systems. Thus, for example, an analogous 
equation holds for a pendulum acted upon by random 
disturbances and, under certain conditions, for a gyroscope 
as well [2]. 

The basic purpose of the present paper is to compute 
the mean difference 9 of the fundamental and locked 
generators’ frequencies , which is of the greatest practical 
interest. Here and in the sequel, a superscript bar denotes 
the operation of averaging over the set of different 
realizations. 

The Two-Dimensional Fokker-Planck Equa- 





tion 

~~ We use the Fokker-Planck equation to solve our 
problem. It is well known [3] that for the investigation of 
the behavior of a system with time constant tT, on which 
there acts a stationary random disturbance ¢ (t) with 
correlation time T, * T, , one can use the apparatus of 
Markov processes and, in particular, the Fokker-Planck 


equation. By the correlation time we understand the 
quantity T, = \ | Rz()|dt, where Ry (rT) is the 


autocorrelation function of € (t) 
We assume that this condition holds for the system 
described by the differential equation 
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If we consider, instead of the nonlinear differential equation (7), the equivalent system of two first-order equa- 
tions x = ys y = Fly, x, © (t)), we may then, by employing the usual procedure [4, 5], find that the two-dimen- 
sional probability density w(x, x, t) for x and x must satisfy the following two-dimensional Fokker-Planck equa- 


tion: 
0 
0 OF 
Ws = —2 “2 2 {(F + \ K{ Fe] de) ws) 4 


0 
+ 5rr (Ws \ KIF, Fldc}. 


—oo 





(8) 


Here, in computing F and the integrand correlation functions, denoted by the symbol K, the arguments x and x 
are considered fixed. 

In seeking the solution of (8), it is necessary to add to it the initial and boundary conditions, which 
are determined by the physics of the concrete problem. Certainly, the probability density w, (x, x, t) must satisfy 
the normalization condition. By giving some distribution , wy (X, x, t)= Wap (X, xX), at some initial moment, we then 
find, by solving (8) , the subsequent deformation of this distribution for t > t. If the random functions x(t) and X(t) 
can assume all possible values between plus and minus infinity then, as the boundary conditions, one should use the 
convergence of the probability density to zero as x, X -» + 0, Les, W2(+ ©, + @, t)= 0. If the random functions 

x(t) and X(t) assume values which are limited to a certain region, then it is necessary to seek the solution of (8) 
only in this region. In the problem we are considering , (6), we shall use, as the boundary condition, the condition 
that the probability density be periodic in ¢. 

If we limit consideration to a stationary solution of (8), when the probability density wg does not depend ex- 
plicitly on time, i.e., W, = 0, we then obtain the following Fokker-Planck equation for the stationary state: 


Bo | eur, ae} S f(r [LS -eJe)a} 



























As applied to (6), it is now necessary to set 
F =aA,—ap—aAsing + &(t, 9, @2=@, (10) 
We now compute the individual quantities which enter into (9). It can be shown [6] that 






E(t, 9) = ¥+a¥+4%% (E, cosp —E sing) = 0. 
Therefore, 





F = aA, —ap—aAsing, 











K|F, F:)=|F —F| |F:—F.] =&@ DEE +4 9 = 





= ($+a$+ ‘= (E, cos p—E; sing))[ $e +a b, + SO (Ex cosp— E;,8in ¢)]= 


=$ be +204 da" b 44+ (S)or%r@, 


where it is assumed that the correlation function of the noise n(t) has the form 


Ry (t) = 07 (t) COS Wot. 
From (10) we also have 
See Fo. KIS) Pleo 
o9 on va ee" Es ‘|= 


By substituting the individual quantities in (9), we obtain the following definitive form of the Fokker-Planck 
equation which describes the stationary operating mode in the PLO system under consideration: 


K =e | 






a A . 2 Ow. 
A ears {(ado — ap — aA sin ¢) w,) — pt = 0. (11) 













Ix° 


11) 


Here, 


—___ 


Ki = | (9 $4 20d de +a%h d+ (F)'ot (x) de. (12) 


—0o 


Second-order partial differential Eq. (11) is linear and of parabolic type. Its general solution depends on two 
arbitrary constants, which are defined by the condition that the probability density be normalized and be periodic in 
¢, Lee., 


—co 


2m oo 
| | wi, dede = 1, wa(g, @ + 2m) =u (>, 9). (13) 


An exact solution of (11) for (13) exists only in the particular case when the initial mistuning of the generators 
equals zero: Ag = 0. We easily convince ourselves that the solution of the sings 


Sats 4+ — = (ag +- aA sin 9) we} —9% = 0, (14) 


satisfying (13) has the form 


. = x i aq? aA 
wa (p, 9) = VY InKy aay eXP s oF + —— cos e| ; (15) 
2n] o( K ) . 8 
2 
where I)(z) is a Bessel function of zero order and imaginary argument. 
From (15) we find the one-dimensional probability density for the frequency difference ¢: 








w (9) = \ W, (?, 9) dp = V =x. exp (— Se) : (16) 
From this we easily obtain 
0 
+r a 4 a le sae 
9=0, 0 =gt=—K, ( {RC bbe+ po bet ne (az) 7} de. (17) 


In many practical cases one can neglect the generators natural fluctuations in comparison with the external noise. 
Then, the basic role in (17) will be played by the last term in the integrand. 

Thus, if the initial mistuning of the generators at the moment when the PLO fs switched in is zero, then the 
average frequencies of the fundamental and locked generators remain equal to each other during PLO operation, and 
the dispersion of he frequency difference o *%, as is obvious from (17), increases with increasing intensity of the 
external noise o” and the generators’ natural fluctuations, and decreases with increasing time constant of the RC 
filter. These results are comprehensible physically. 

In the sequel we shall consider approximately the cases when, in (11), the initial mistuning of the generators is 
different from zero (Ay # 0). 

Approximate Method of Solving the Fokker-Planck Equation 

It is difficult to find a solution in compact form of (11) which satisfies (13). Below, we present an approximate 
method of solving (11). 

We shall seek a solution in the form of a series: 





ws (, 9) = Dy w™ (9) La (@)s (18) 
n=0 


where the functions L,, ( ¢) are to be defined, and 


) = exp (— %* 
w(@) = exp (— 3p). (19) 
By direct differentiation, we easily convince ourselves that the following recursion formulas are valid for the 
derivatives of the function w(¢) : 


gu (G) = — 4 Kw" G) — pw) (m= 0.4.2... (20) 
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By substituting the expression for w,(%, ¢) from (18) in (11) , then using (20) and collecting similar terms, we 
obtain 


= K, ) w+” @) L, (@) — a (bo — Asing) 5) w™*” ) Ln () = 


=a >) nw (9) La (9) — >) nw” (9) Ln (9) 
n=0 


n=0 


or 


J wi) Ln) — Fy (40 — Asing) 3} w"*” @) La (@) = 


n=0 
2 or “a = — — 21 
= Dnw™ ) Ln (@) — Fe Dy rw"? G) Enq). se. 
n=0 n=0 
By equating coefficients of like orders of the derivatives w°") (9), we obtain the following system of dif- 
ferential equations for the definition of the L, (¢): 
L; (9) = 0, 
’ 3 3 , 
Lo(9)— (Ao — Asin 9) Log) = F-Li (9) —2 ZL) 
’ a? 2 , 
Ly (9) — (bo — Asin) Ly (9) = 2 FL (9) —3 Z-La (9), 
, 2 ; 2 , 
La (9) — 3; (So — Asin 9) Ly (9) = 3 &- Ls (9) — 4 - Li), 
a re ee a a ee we Or eh ae ee a ee ek ee et ee a (22) 


In general, 


Ln—1 (9) — - (Bo — Asin ¢) Lua (9) =m GLn(g) — (n+ 1) Lng (@). 


We may suggest the following order for solving (22). The equations are solved successively , starting with the 
first , whereby, in solving each successive equation, we ignore the last terms in the right member, containing deriva- 
tives of higher-order functions Ly (¢)- 

For an approximate estimate of the value of the mean frequency difference of the generators, we shall henceforth 
limit ourselves to computing the first two functions, Ly(¢) and Ly(¢). From the first equation of (22) we have that 
Ly = c = const. By neglecting the term containing Lj (¢) in the second equation, we obtain the following differential 
equation for defining Le (¢): 


Ly () — %- (Ao — Asin g) Lo (¢) = $c. 


The solution of this equation, as we easily convince ourselves, can be written in the form 
? 


3 
L, (¢) = cx exp (D,9 + Dcos ¢) \ exp (— Dyy — Dcos 4) dy, 
Ve 
where @p is some arbitrary constant and 
oe aA, 
a 


Thus, in our approximation, we obtain the following solution of (11); 


Ws (, 9) = w (g) Lo (¢) + w’ (9) Lig) = 


, D=—. (23) 


' (24) 
a ao? . 4 
= ¢ &- exp (~ ar) {—9 + aexp (D,p + Dcos¢) \ exp (— Dyy — Dcos*) dy}. 


oe 
To satisfy the condition for periodicity (13), it is necessary that the following equation hold : 


o+2n 
exp (Dop + 2D, +- Dcos ¢) \ exp (— D,y — Dcos4) dy = 


Pe 
? 


= exp (Doe + Dcos #) | exp (— Dox — Dcos 4) dx. 


Ve 
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23) 


24) 





From this equation, the following relationship derives: 


o otin 
xp (2nD, 
{ exp (— Dyy — Dos) dy = 7— ae ae \ exp (— D,y— D cos) dy, (25) 
Pa + 
which permits the arbitrary constant @ to be eliminated. 
To determine the second arbitrary constant c , we use the normalization condition 








2m co _— 2n 
° . 2" 2xD 
\ Ws (p, p) dedp = ca x. i oe ea \ exp (Dop + D cos 9) dp x 
ofan ; (26) 


x exp (— Doy — Dos 1) dy = 1. 


It can be shown [1, 7] that 
on ot+an 


J= \ exp (Dop + Dcos 9) dp \ exp (— Dox — Dos x) dy = 
0 ? 
= 4n*exp(— Dp) | lin, (D) |, 


where Ij» (z) is a Bessel function of imaginary order and imaginary argument. Therefore, we obtain from (26) 








— ‘Ky 1—exp (2Dpo) 
C= FaiV tea exp(eDa) | /t0.(D)[* (27) 
By using (25) and (27) , we can write the approximate expression for the probability density in its definitive 
form: 
: fT 1—exp (2nDp) Boy’ ag? 
ws (? 9) = V 2naK, 4n?® =: tmDe) | Tin, (D) [* exp (— 2K x 
exp (2nDo) tg (28) 





x{-¢+e, exp (Dyp-+Dcos¢) | exp (— Det — Dos y) dy}. 


— exp (2xDpo) 


It is easily verified that this formula goes over into the exact formula (15) for zero initial mistuning (Dp = 0). 
This attests to the fact that the suggested approximate method of solving the Fokker-Planck equation gives logically- 
noncontradictory results. 

By integrating (28) over g and # between the proper limits, we obtain the one-dimensional probability densities 
for, respectively , and gy; 





=——— ’ (29) 
Cai V a ° 2x 1—exp(2xDo) » __ a9 
w (9) -{ anKy -39V aK, 4n?® anaes | 4in(P)| +} exp ( K; ). 
e+2n 


w (>) = — exp (rDp) | Jip, (D) |"? exp(Dop + D cos @) \ exp (— Doy— Deos x) dx. 


The last formula for the probability density of g coincides in form with that obtained in[1]. The difference 
consists in this, that here the coefficients Dy and D depend on the low-frequency filter's passband[cf. (23)}. This 
fact also gives us some confidence that, in taking only the first two functions in (22) into account, We obtain results 
which are in qualitative agreement with those obtained in [1]. 

By using (29), we find the average difference of the generators’ frequencies and the dispersion of this difference; 


(30) 


Ky; exp (2xDp) —1 


iS BSC oe sh xD, 
= °° exp (nDpo) 


De Ip, (D) |, 0. at Ky. (31) 








| Tip, (D) % = Ay 


213 





For Dy = 0 , these formulas coincide with (17). 
Particular Examples 


Let the externai noise n(t) be obtained as the result 
of the action of a sufficiently wide-band fluctuation on a 
narrow band oscillatory loop withresonant frequency wp 
and, consequently, let the noise have correlation function 

k(t) = ote—F!*! cog wot. 

We assume that we may ignore the natural instability 
of the frequencies of the fundamental and locked genera- 
tors in comparison with the noise n(t) Then, 


K,= (aA +) \ Of ds = = (as +): 


oo 





By substituting this value of K, in (23) and (31) , we 
BAo /A;\* B /A,\* 
r= (S). =a (FZ). 


a 
—A = FHAAs, o: = 7 A? (=) ° 
In the given concrete case, Dy and D do not depend 
ona. With this, as is clear from (32) and (31), the 
dispersion of the frequencies’ difference depends on a, 
but the average discrepancy of the frequencies does not 
depend ona. This last result is seemingly a consequence 
of the approximate solution of (11), Recently, A. G. 
Zhuravlev showed that, although the term containing 
L2( ¢) in (18) for the probability density does not of 
itself give a correction to the mean frequency difference, 
its being taken into account does result in the function 
Lo (¢) assuming a different form: 


a@ exp (2xDpo) " 
Io (9) = C RT exp (2%Dy) °*P 
o+ar 


get 


(32) 





« (Dog + Deos 9) (a — Acos7) X 


(33) 


x exp (— Doy — Dos 7) dy. 
We now make some numerical estimates. It is 


of the Fokker-Planck equation. Then, A * sAq = 100 cps, 
Dp = 1.5 (Ay/0)',D =3(Ay/0) and oy * 700@/Ay), 
By using the graphs of I. G. Akopyan, given in[1], we get 


as A,\* 
W — Wo = 35cps, 0. = 36 cps for (2) x 0.7, 


Ay 


5 a 
w — @ = 10cps, o. =18cps for (2) = 2.0 


From what has been considered above, we may draw 
the following conclusions : 

1, From (32) and the quantitative estimates adduced, 
it is clear that, when the conditions hold for the applica- 
bility of the Fokker-Planck equation (a<< 6 ), the mean 
frequency difference and its mean square decrease as the 
ratio A;/o increases. For sufficiently small signal -to- 
noise ratios, the average frequency difference may turn out 
to be significantly less than the mean square value, and 
it is necessary to take it into account. 

2. When stability condition (34) holds, a decrease 
of the time constantRC provides a shorter duration of the 
PLO's transient response, but the dispersion of the frequency 
difference is thereby increased. Therefore, the 
requirements of short establishment time and high noise 
immunity are, in a very definite sense, contradictory. 

In conclusion, I wish to thank R. L. Stratonovich for 
his interest in the present work. 
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Unax _ + uArAy g 2A. 
We take U4, ¥ 0-5U max » Lee. AgAs /2 Ay. 
Whence p 2/ A, and Aw sAy. 
The parameters of the PLO scheme are so chosen 1. 
that stable operation is guaranteed. As follows from [9], 
forAy/A < 0.6, stability is guaranteed if the following 


inequality holds; T_ 2h 8. 
RCA 7> TT: (34) 


Let Ay  lv,s=O.1lkg/v, a = 20 sec™!,8=300sec1, 9. 
Ae/ 4 = 0.5. With these parameters, Ineq. (34) holds and, 
simultaneously, the conditions hold for the applicability 
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ANALYSIS OF A TWO-POSITION COMPENSATING CONTROL SYSTEM 
UNDER CONSTANT PROLONGED DISTURBANCES 


A. I, Cherepanov 
(Pert’) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 310-316, March, 1960 


Original article submitted October 28, 1959 


It is shown that the simplified graphico-analytic method developed in [1] for computing the processes in a two- 
position regulator may be employed for analyzing the processes in two-position compensating control systems [2] 


under constant prolonged disturbances. 


In [2] of this author there was described a new 
method for improving the quality of the processes in a 
two-position regulator by the use of a special relay-type 
correcting device. The basic advantage of this method 
in comparison with prior ones [3, 4, 5, 8] is that, with 
essentially reduced amplitude of controlled quantity 
oscillations, of temperature for example, it provides the 
minimum increase in the required number of switchings 
of the object's heater in and out of the system, computed 
on the basis of a nonsymmetric mode of operation. 

It should be mentioned that, at the present, of all 
the types of two-position compensating controllers, those 
with thermoelectric correcting (compensating) devices 
have been most widely used (6, 7]. 

However, electromechanical correcting devices have 
certain advantages: e.g., such devices can be used both 
in regulators which employ some electrical means of 
measuring the controlled parameter and in many industrial 
two-position regulators of dilatometric, bimetallic, or 
magnometric types which use some nonelectric means 
of measuring the controlled quantity. Moreover, the use 
of an electromechanical correcting device permits the 
employment of the regulator's measuring portion for 
measuring or recording the true temperature in the object. 
When a thermoelectric device is used, it is necessary to 
have an additional intrument for this purpose, since the 
action of the device distorts the output signal of the 
measuring portion. 

Simple formulas were obtained in [2] for calculating 
the correcting device's adjustment parameters and the 
basic parameters of the compensating control process in 
the steady state. It was noted in[2] that the character of 
the possible disturbances in the system should be taken 
into account in choosing the correcting device’s parame- 
ters. 

In the present paper we provide an analysis of the 
processes with disturbances in a two-position control 
system with an electromechanical correcting device. The 
possibility is demonstrated of the existence in such a 
system of complex periodic motions, somewhat similar to 
those considered in [5, 8]. Recommendations are given for 
the choice of the adjustment parameters of the correcting 


device which will provide the very best quality processes 
given the disturbances. 

As in[2] , the analysis is made by means of a 
simplified graphico-analytic method of calculating the 
two-position controller's processes [1], thanks to which the 
formulas obtained are simple, and may be used without 
difficulty for engineering computations. 

By the adjustment parameters of the correcting device 
we shall understand the following quantities: tet, is the 
duration of the make state of contact K-1, tye is the dur- 
ation of the break state of contact K~1, thre = tytg tmak 
is the total duration of the make state of contact K~1 
during one disk rotation, trak = &ty+ ts — tyte, tho Is the 
time for one complete disk rotation, tro = tots - 

1. Pulse Repetition Period as a Function 
of the Load 

The distinguishing feature of a two-position control 
system with an electromechanical correcting device fs 
that the repetition period of the pulses of heat admission 
to the object, for invariant adjustment parameters, depends 
on the load. This is achieved due to the use in the device 
of two controlling disks, one of which switches in the flow 
of energy to the object at a definite previously specified 
time while the other disk returns the first disk to its 
initial position in the time following a positive deviation 
of the parameter from its given value* . With this, if 

the magnitudeof the energy efflux changes within certain 
limits, for example, if it increases, then the new switching 
in of the relay's motor would occur somewhat earlier 
than for the previous value of the efflux. The pulses 
switching in the influx would have the same duration, but 
their repetition frequency would increase, and the 
correction process would not be ‘disrupted. 

Figure 1, a shows the graph of the process in a two- 
position controller with compensation, where the influx 
multiplicity is ny = 8, and the degree of correction is 
k = 2,86. In Fig. 1 and in the sequel, we adopt the 
nomemclature; 4 ® is the deviation of the controlled 
parameter (temperature), A 9’ is the deviation of the 
object's controlled temperature with compensation 

* For a more detailed description of the correcting 
device's operation, cf. [2]. 
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(correction), A9,,- is the temperature deviation of the 
sensitive element with compensation, A 9'r4) 1s the 
positive amplitude of the deviation with compensation, 
A §m(-) is the negative deviation amplitude, Q;, is the 
influx of heat, Q.¢ is the efflux of heat, Tj, is the 
duration of switching in for control without correction, 
Tn 1s the duration of switching in when the corrective 
device is used, Toyz is the duration of switching out, T 
is the period of oscillation, t) is the moment of switching 
in, ty is the moment of switching the influx off when 
correction is used, t, is the moment of the new making of 
controlling contact K~-1, t, is the moment of disk cessation, 
At is the system's lag time, At" is the residual lag time 
when correction is used. 

Figure 1,b shows the control process in the same 
system when the heat efflux is increased twofold, L.e., 
for n, = 4. As is clear from the graph, the correction 
process is not disrupted. The duration of switching in of 
the heater, T';,, , remains as before. The duration of 
switching out T*,,,,, 1s decreased, and becomes equal to 
Toute ° The positive amplitude of the deviation is also 
decreased. The negative deviation amplitude AFn- 
increases, as with ordinary two-position controllers. ihe 
period of influx switching in is shortened to T',;. The 
total amplitude of the oscillations, A 9/4) + A Sny-) 
is somewhat decreased. 


It may be seen from a comparison of Fig. 1, a and b 
that the degree of correction k became larger in the 
second case; k = 6,2. 

We now find the dependence of the variation in the 
degree of correction k on the magnitude of the influx 
multiplicity n for a constant (given) duration of switching 
in Ty, - As is known [2], 
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At 
k == AW . (1) 

The quantity At’ can be expressed in terms of the 
duration of switching in T'j, and the influx multiplicity 
n by use of the formula [2] 





Qin ef 
Tin = At’ ~-——~— + (At — At’) ——_. 
in Qin as Qef ( ) Qn inte Q ef (2) 
From (2) we obtain 
ay 1 
At =T in — At —,. (3) 
By substituting (3) in (1), we find that 
At 
k =x \ eon i a 4) 
Tin — Ate 


From (4) one can calculate the degree of correction 
for various values of influx multiplicity n and switching 
in duration T’;,, « 

Equation (4) shows that, for some value ofn , the 
denominator vanishes and k becomes infinite. 

This value ofn will be, in some sense, a critical 
one; we denote it byn,.t. Obviously, 

n,,=1+ a 
in 

When the influx multiplicity has the value n,,, the 
temperature in the object during the given time T";,, 
rises exactly to the given limit. For values of n less 
than n.,, the temperature in the object during the 
switching in time T*;,, does not attain the given value, 
and the normal course of the correction process is disrupted 





(5) 





T In[2) no, has a different connotation. 








(1) 


city 


(2) 


(3) 


(4) 


*tion 
hing 


al 


(5) 








For the case shown on Fig, 1,a, the calculated value 


of Nor = 3 On Fig. 1,b, the dashed line represents the 


temperature variations in the object for this value n = 3+ 
The calculated value of n,, is well borne out by the graph. 


Thus, in the system under consideration, the repetition 
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2. Complex Periodic Motions 
As was shown in the previous section, the correction 

process can proceed normally for disturbances up to some 

value ny With this, as the disturbance increases (the 
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period of the influx switching-in pulses depends on the 


load, and the correction process flows normally only for 


influx multiplicities greater than the critical value. 


For values of influx multiplicity n< n., 
course of the correction process {is disrupted , and there 


, the normal 


arise complex periodic temperature oscillations in the 


system, examples of which are considered in the following 


section. 


























value of n decreases) there will be an increase of the 
negative amplitude and a decrease of the positive am- 
plitude. The given condition (n > n,,) for the guarantee 
of normal correction process is necessary but not sufficient. 
As is clear from Figs. 2, 3,and 4, it is necessary, for the 
normal flow of the correction process, that the following 
additional conditions be met: 







t>Tin. to <7. (6) 
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If the correcting device's parameters, selected so as 
to permit processing of the greatest possible disturbances 
beso chosen that inequalities (6) are transgressed for values 
of n higher than nor, then complex periodic motions, 
similar to those considered in[5, 8], will arise in the 
system. 

Examples of such processes are shown in Fig. 2 and 3. 
Figure 2,ashows the process with correction for n = 8 up to 
line cd. At the moment which corresponds to line cd, 
and instantaneous twofold increase of efflux (n = 4) occurs. 
As is clear from the graph, after several oscillations the 
process is reestablished with the new period T, With 
this, the correcting device completes two revolutions 
during a period. It can be assumed that the form of the 
postdis turbance -reestablished oscillations does not depend 
on the moment when the disturbance {s applied. This 
is confirmed graphically on Fig. 2, b. 

Figure 3 shows the control process with correction up 
to line ab forn=8. At the moment corresponding to 
ab, a disturbance is applied (a fourfold increase in heat 
efflux) and the process continues with n= 2. The new 
period of complex oscillations Ty is, in this case, some- 
what less than in the previous case. 

Thus, complex oscillations occur in the system if one 
of the following conditions is not met: n> Mor, t2> Typ» 
to < T’- 

3. Optimal Parameters 

We now seek the optimal adjustment parameters of 
the correcting device for a given accuracy of control 
(degree of correction k), starting from the condition that 





the correction process not be disrupted for the very greatest 


disturbance, i.e., from the condition of minimal Dor * 

By Nom We shall understand the critical value of 
influx multiplicity below which the controlled parameter 
(temperature) can not be maintained at a given level. 


This occurs due to an excess of efflux of heat over influx. 


For ordinary two-position control, n.,, = 1 when Q,, = Qi 


For regulation with correction , disruption of control 
occurs for a value of Nem > 1, since the influx of energy 
to the object is switched offindependently at some time 
three the case when the temperature in the object 
remains below the given value, the correcting device 
works continuously. ‘During one rotation, the amount 
of heat applied to the object equals Qintmak Am equal 
quantity of heat will flow out with the condition 

Qin lmak= (tbre+ 4nal Gef> 
where Qer is the critical value of heat efflux. If 
Qer> Qos, then the temperature does not rise to the 


given level. 
From whence we find Dom 48 


5 Q in tmakt lbre 

alt le ee 
Qef 

It is clear that tak = tro ~ thre» We finally obtain 
t bre 

o ‘bre 


wit 2e.. 


'mak mak 


Nom= 1 + i (7) 








tro 
j=. 
bre 


Neom= Ay » where 
We have (5) for n,; Figure 4 shows the control 
process for the critical value of influx multiplicity 
N = Mc = 3. For this critical case, as is clear from the 
graph , (6) can hold only for 
lL, = Tin 


’ ’ 8 
tio = 1" = Tin ®) 


Thus, (8) gives the optimal adjustment parameters 
for the correcting device. 

We now find the value of ner for the case of op- 
timal adjustment . By substituting (8) in (7), and by 
taking into account that t,, = T'jn + t,, aid that 
thre =At, we find that 


Tin + At 


ol 
in 


(9) 








"em ’ 
i.e., in the given case, Nar, = Mer + 

The correcting device's adjustment parameters 
corresponding to (8) should be taken in those cases when 
very protracted disturbances are not observed in practice 
and n is always greater than n.,- We note that complex 
motions cannot arise in the system under these conditions. 

If very protracted disturbances are possible in the 
control system (n< Nor), then the device's parameters 
should be chosen from the condition that the system be 
capable of processing these disturbances. 

As is clear from (8), this can be achieved by 
increasing the time of disk rotation t,, or by decreasing 
the time of the contact's break state t},,.- In practice, 
for thermal objects (for example, of the dessicator 
hood typé), one might recommend that the rotation time 
be chosen within the limits of (0.6 to 0.7)T',since with 
this are guaranteed a rapid rise in temperature at the 
initial moment of switching in, and good processing of 
the disturbance (for example, upon opening the door 
of a dessicator hood). 


CONCLUSIONS 


1, The simplified graphico-analytic theory of 
calculating processes in two-position controllers may be 
used for the analysis of processes with disturbances in 
two-position controller systems with time correction. 

2. The distinguishing feature of the two-position 
control system with an electromagnetic correcting device 
is that the pulse repetition period (the mark-space ratio) 
depends on the load. 

3. Under definite conditions, complex periodic 
motions of the controlled quantity (temperature) arise in 
the system. 

4. With optimal adjustment parameters of the 
correcting device, normal course of the correction process 
is guaranteed for the very greatest disturbances. 
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ON THE ERRORS OF AN INTERPOLATING DEVICE 


FOR A GIVEN CIRCULAR MOTION 
V. V. Karibskii 


(Moscow ) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 317-322, March, 1960 


Original article submitted July 7, 1959 


An algorithm for a digital interpolating device for a given circular motion is described. A theoretical investi- 
gation of the interpolation errors is given. A brief description of the device's circuit is provided. 


Today, in various domains of technology, the neces- 
sity arises of so controlling an object that it moves along 
a given trajectory. Frequently this is implemented by 
means of sampled-data servosystems. Each system 
controls the motion corresponding to one coordinate and, 
at each input pulse, there results a translation by unit 
length in the scalechosen, this length being called the 
value of the pulse. The pulses for each coordinate axis 
must be distributed in accordance with a given trajectory 
of the object's motion. Distribution of the pulses is 
ordinarily implemented by means of special digital 
computers. 

In the present paper we consider the design 
principles, and estimate the errors, of a device of this 
type which provides motion along a circle. 

The equation for a circle centered at the origin has 
the form 

x? + y*? = RP, (1) 
where x andy are the current coordinate values and 
R is the circle's radius. 

We assume that the motion begins at point A with 
coordinates x = 0, y = R (Fig. 1), and must proceed 
clockwise about the circle. 

We now find the value of x” for which the coordinate 
y is decreased by unity, i.e., becomes equal to y=1, for 
which we substitute the new value of the coordinate y 
mets 2— R2 _4y 

xz? = R°—(y — 1)’. 

By taking into account that y = R, we get xy" = 2R-L 

Analogously, we can find the values of xs”, | ee 
x for which there are corresponding changes of the 
coordinate y by 2, 3,..., k units. Knowing these 
values, we can find the corresponding differences, which 
will have the form 


zi = 2R—1, 
28 — 29 = 2R—3, (2) 
Te — thy =2R—(2—1) (k =4, 2, 3,...). 
Based on the formulas obtained, we can suggest the 
following method of constructing acircle. First, we must 
vary the coordinate x until the value of xe equals 2R-1 
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and, as soon as this happens, we must change y by one 
unit, after which we change the coordinate x until the 


difference Xe"=x;" becomes equal to 2R=3, at which time 
we again change y by one unit, etc. It is obvious that, 
in this case, the approximation error is nowhere larger 
than unity. However, due to the discreteness of the 
system's operation, such a method is not implementable, 
since translations of less than the value of a pulse are 
impossible. Consequently, R and the quantity x}? are 
integers. However, the’quantity x}; may be expressed 
by a noninteger at points where y, is expressed by an 
integer, and conversely. We shall, therefore, in moving 
from point A, vary coordinate x by unity and then 
compute the value of x*, As soon as the coordinate x 
reaches the integral value x;_,, for which the condition 
xi-m = 2R-1 holds, we change coordinate y by unity 
and determine whether the condition holds that 
A, = Zj3_-m—(2R — 1) > 2R—3. 
If so, then a unit increment of coordinate x at the given 
point of the circle corresponds to an increment of co- 
ordinate y greater than, or equal to, two units. There- 
fore, coordinate y should be changed by one unit, the 
quantity 4, = zj-m—(2R — 1) —(2R — 3) com- 
puted, and the condition A, =(2R-5) verified etc. Such 
operations are carried out until the condition 
Ay < 2R-[2(k + 1)—1] begins to hold. As soon as 
this happens, we again change coordinate x by 


unity and then compute the increment of x? obtained, 
with account being taken of the residue, i.e., in this 


2 
case X4¢1-m—* tom +A, where X44 2-m7%jem = 1h 
If now the following condition holds: 


Tt1—m — Ti-m + Ay > 2R—([2(k+1)—14}, (9 


then coordinate y should be changed by one unit, and the 
process just described should be repeated. The sequence 
of actions described above and (2") comprise the algorithm 
in accordance with which the scheme described at the end 
of this paper operates. In this case, naturally, coordinate 
y will be changed by one unit somewhat later than the 
time when x reaches a value corresponding to a one-unit 
change of y , rather than simultaneously. As a result, the 
circle will be approximated by a broken line with an 
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error. We now estimate this error. For this, we prove that, 


with the method chosen, varying the magnitude of x by 
unity at any point x, in the interval er x< R/ ¥3 cannot 
give rise to a pore. of the quantity x ie i greater 
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Fig. 1. 


than what is necessary to change y by one unit. 

A. With motion along an approximating curve in 
accordance with the method we have adopted, let the 
moving point reach point K (Fig. 1) where, at this 
moment, coordinate x is changed by x; units and 
coordinate y by k units. Then, a subsequent increment 
of coordinate x by one unit engenders the increment 


Thay — SF = 27; f- 1. (3) 

For coordinate yj to be changed by one unit, with 

the condition that A,=0 at point K, it fs necessary 
that the following equation hold; 

Tay — Xi = 2R — [(2(k + 1)— 1]. (4) 


It follows from this that for all relationships of x 
and y which satisfy the condition 


2R — (2k + 1) > 22; + 4, (5) 
an increment of one unit in coordinate x will give rise 
to an increment of coordinate y equal to or less than 


unity if, at these points, it can be assumed that A; = 0 
After transformation, (5)assumes the form 


R—k>x,+1 (6) 
or, taking into account that y; = R-k, 
Yi > Ti+1- (7) 


For x equal to y;, a unit increment of coordinate x 
engenders , for A, = 0 , 4 unit increment of coordinate 
y . which corresponds to an angle of 45 degrees, or 
x =R/¥2. In other words, for motion from point A to 
point K, each unit change of coordinate x engenders a 
smaller variation in y . Only, when , as the result of 


successive variations of x by unity, coordinate x becomes 


equal to coordinate y at the preceding point does the 
coordinate y also change by unity. 

B. We now determine the maximum possible de- 
viation of the approximating curve along the axis 
of ordinates from the given circle in the interval 


0< x < R/ V%, with account being taken of the Ay 
Let the motion begin at point A (Fig, 1) and be directed 
clockwise, with the coordinate x being given a unit 
increment. By virtue of what has already been proved, 
each unit increment of coordinate x requires a variation 
of coordinate y by a quantity less than unity. 

However, in accordance with the method chosen for 
constructing the circle, this change does not occur until 
he resulting increment of x? satisfles the condition 
x*> 2R-1. This means that, for such a value of x, the 
deviation along the axis of ordinates of the approximating 
curve from the given one is 6, = 1. Let such a point, 
where the condition x*> 2R-1 holds, be the point with 
abscissa xj.,,- Then, at the point x;__,_, where 
%-m7*{-m-1 = 1, by < 1. But by virtue of the fact 
that on this segment a unit increment of coordinate x 
cannot engender an increment of coordinate y by more 
than unity, the transition from point x4_-,-1 to point 
X4-r C40 increase 6, by at most unity. Thus, the 
resulting error along the axis of ordinates can lie within 
the limits 0=6,< 2. It should be kept in mind that, 
after transition to the point with coordinates (xj~m, 
Yj-m)» there occurs a change of y by unity, and 6 
becomes less than unity. 

By an analogous argument we can successively check 
each value of x in the given interval, coming to the same 
conclusion each time. 

We now consider the very worst case. Let the moving 
point be at point K with coordinates (x; , yj) such that the 
next increment of coordinate x by unity engenders the 
incrementation of y by unity for Ay = 0% be, 

71,17", = W-(2k+1) = 2yy-1, where k is the number 

of units by which coordinate y was changed at the moment 
of transition to point K. We assume that, at the moment 
of transition to point K, the total accumulated increment 
is x" =x4-1 , where x;-1 is the abscissa of the point at which 
the last change of y occurred, or that the magnitude of 
the residual A, . if we hit point K after the latest change 
by unity of coordinate y rather than x , had the maximum 
possible value A; .,,,* In accordance with the algorithm 
for the operation of the scheme, Ay may Must, at point K 
be at least one unit less than 2yj-1, Le., must satisfy the 
equation Ay max = 2¥j72- 

In making the transition from point K to point L, the 
residual Ay, max increases by 2y;~1, as a result of which 
coordinate y changes by unity and executes a transition 
to point L’. However, the residual at point L’, Ay, 4 = 
= Ax max: turns out to be larger than the quantity, 
2yj-1)-1, necessary for changing the ordinate of point L’ 
by unity. As a result of this, a transition occurs to point 
M, and the following residual is formed: 


Ante = 2y¥1 —2—2(y —1)+1=1. 


y 


A third change of coordinate y by unity can only 
occur for the condition 2(y;-2)-1 = 1, Le., for yy = 3. 

By taking (7) into account, we can write x41 Sy, = 3 
It follows from this that variation of coordinate y at 
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point L by three units in succession is only possible for 
R= 42 xj, $372. For all values of R>3V2, 
change of coordinate y at point L can occur only by two 
units. For large R (R s 100), the magnitude of 

2 (y,~2)-1 » 1 and, by neglecting unity in comparison 
with 2 (yj-2)-1, we can assume that in the worst case the 
change of coordinate y is by exactly two units, and point 
M lies on the circle. Thus, the maximum error, equal 
to two, will be at point L. 

If the order ofapplying incrementation is somewhat 
varied, this error can be reduced to unity. For this, we 
shall first carry out all the computations and then, when 
it becomes clear that a unit incrementation in y must 
occur, we execute it simultaneously with the incrementa- 
tion in x which was giving rise to i. In the actual scheme, 
this leads to motion at an angle of 45 degrees from point 
K to point L" (Fig. 1). As a result , the errer in coordinate 
y » equal to by: will everywhere be less than unity. 

The expression for x in terms of y which is obtained from 
the equation of the circle coincides with the expression 
for y in terms of x. Therefore, all the argument already 
given for x fs also valid ary within the limits of varia- 
tion of angle a from 0 to 45°, i.e., for variation of a 
within the limits R = x = R/V, and the approximating 
curve is symmetric about the line passing through the 
first quadrant at a 45° angle. At the limits of the interval 
0< x < R/ V2, the radial error at any point x will equal 


bn = V (VR?—2? + dy)? +2°— RK. (8) 
The error 5,, equals the difference between the true 


value y corresponding to a given x and the next larger 
integer if y itself is not an integer. 
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2 r From pulse 
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Fig. 2. 


Figure 2 shows the simplified block schematic of 
the interpolating computer which operates in accordance 
with the principles described above. The circuit consists 
of counter 3, to which quantity x is fed, adder 1 which 
stores the current value of the difference x 2x} + Ay, 
and adder 2 which stores the current value of 2R-[2(k+1)-1}. 
The dashed lines on Fig. 2 mark off the adders" least 
significant places. If the motion is clockwise, starting 
from point A (Fig. 1) then, before operations commence, 
the counter and adder 1 are set to zero and the number 
2R-1 is placed in adder 2, Then, from the pulse generator, 
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one pulse is applied to the x counter, storing a"1" in the 
counter and in adder 1. From the contents of the adder 

1 the number 2R-1 {s subtracted and, if the difference 
formed in adder 1 is negative, the number 2R-1 fs added tp 
it, and one pulse is applied along the x axis. The sign 
of the difference is staticized in the sign bit of adder 1, 
After this, another pulse is applied to the counter, and 
the number 3 is added into adder 1. As a result, the 
number 2 is formed in the counter, and the number 4 (2°) 
is formed in adder 1. Formation in adder 1 of the square 
of the number is implemented automatically, thanks 

to a special scheme for connecting the counter with 
adder 1, suggested by B. L. Ermilov[1}. Then the arith- 
metic acts just described are repeated. If the difference 
turns out to be positive, then addition does not occur, one 
pulse each is applied simultaneously along the x and y 
axes, and the number 2 is subtracted from the contents of 
adder 2, i.e., in adder 2 there remains the number 2R-3, 
This number is subtracted from the contents of adder 1 
and, if the difference is positive, then another pulse fs 
applied along the y axis, and the number 2 fs again 
subtracted from the contents of adder 2, etc., until the 
difference becomes negative. As soon as the difference 
does become negative, addition again ensues, and the 
next pulse is applied from the pulse generator to the 
counter, etc. 

If motion begins from some arbitrary point with 
coordinates [x%. Ye]. then the initial establishment of 
values must be the following : in the counter,the number 
Xe is established, in adder 1,the number xp=R? + Yq, and, 
in adder 2, the number 2y, - 1, where yp, is the value of 
y corresponding to the value x, and rounded off to the 
next higher integer. 

For processing the complete circle, the initial estab- 
lished values and the direction of motion along the axes 
are changed in making the transition from one quadrant 
to the next. 

Thus, by means of addition and subtraction, compar- 
isons are made in the scheme proposed of the quantities 





xis -=" j+ 4k with the quantities 2y,-1, where x; 
and y, are the coordinates of any point of the approx- 
imating broken line. In a comparison the indication 
(sign of the difference ) is developed which is used for 
controlling the distribution of the generator pulses along 
the x andy axes in correspondence with the given law 
of motion (circular) With such a computational algorithm 
for circular motion there occurs a variation of the path 
velocity by a factor of V2 which, in the majority of 
cases, is admissible. Moreover, by using velocity feedback, 
one can so control the frequency of the pulse generator 
that the velocity can be kept constant and equal to a 
given value. (For a more detailed description of the 
scheme, cf. [2].) 


The scheme suggested by B. L. Ermilov [3] operates 
by such a mathematical algorithm, but with the use of 
another computational algorithm, the mathematical 
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algorithm not being described in general form. The 
scheme of B. L. Ermilov consists of two counters and one 
adder of the accumulative type. Pulses from the generator 
are applied to the x axis and to the input of one of the 
counters, forming the quantity x in it; with this, the 
quantity x? is formed in the adder. The adder’s overflow 
pulses are applied along the y axis. By means of the 
second counter, each overflow pulse introduces the 
quantity 2Ma(2y;-1) into the adder, where n is the 
number of places in the adder, adding this number to the 
number remaining in the adder after transmission of the 
overflow unit. If the overflow pulse does not give rise 

to a second overflow pulse then, from the first counter, 
the number 2x-1 continues to be added with the arrival 
of each generator pulse until overflow again occurs, 
whereby since x = 1, 2, 3,..., x" is formed. 

Since, on the segments of the circle where the 
coordinate x is close to R, there are emitted Y2R-1 
successive overflow pulses with the frequency of the 
scheme's elements’ operation, there occurs a sharp change 
in the path velocity. Indeed, the generator's frequency 
must be approximately Y2R=1 times less than the fre- 
quency of the overflow pulses so that V2R-1 overflow 
pulses can be disseminated between any two generator 
pulses. As a result, motion at the beginning of the 
processing of a quadrant of the circle (from 90 to 45°) 
will occur with a velocity smaller by a factor of V2R-1 
than that at the end. Such a sharp variation of path 


velocity in cases when the device is designed to guide 

a cutting tool is inadmissible. It should be mentioned that 
control of the path velocity by variation of the generator 
frequency is not possible in the scheme given, since the 
frequency of the overflow pulses does not depend on the 
generator frequency. 


SUMMARY 


1. For constructing a circle, we use the derived 
formula 73,, — xi + Ay>2R — [2(k + 1) — 4]. 

2. In constructing a circle in accordance with the 
computational algorithm suggested, there occurs no 
accumulation of errors, and one obtains an absolute radial 
approximation error 5p which fs less than one unit in the 
chosen scale. 

3. The scheme suggested makes it possible to control 
path velocity by varying the pulse generator's frequency. 
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The results are given of an experimental investigation of rubber-fabric membrane characteristics, these 


membranes being used in pneumoautomation devices. 


1. Aim and Methodology of the Investiga- 
tion 

Modern pneumoautomation devices, designed for the 
regulation and control of various continuous productive 
processes, are constructed on the principle of force 
compensation. 








As one of their basic physical elements, these devices 
usually contain either bellows or rubber-fabric membranes. 


The pneumatic instruments of unified aggregate systems 
(UAS), produced domestically, as well as the instruments 
of the majority of foreign companies —Moore, Sunvic, 
Taylor, Bristol , Brown, Electroflo, DRD, etc. are cons- 
tructed with membranes. The membranes delimit the 
individual pneumatic chambers inside the instruments, 
and it is on the membranes that forces from the com- 
pressed air in these chambers are balanced. 

Such indicators of quality of operation of modern 
pneumatic devices as their accuracy and sensitivity, as 
well as their static and dynamic characteristics, depend 
to a large extent on the properties and characteristics of 
their membranes. 

There are many works in the literature on the inves- 
tigation of metallic membranes, but there are only a few 
works [1-4] devoted to the investigation of rubber-fabric 
membranes. In connection with this, and because 
precisely these rubber-fabric membranes have attained 
their greatest dissemination in pneumatic devices for the 
control and regulation of productive processes, the 
necessity for a more careful investigation of these 
membranes arose. 








Fig. 1. Forms of the membranes investigated. 


The membranes investigated had plane and conic 
profiles (Fig. 1, a and b), prepared from brand 1520 
caprone, from glass cloth and TD cloth covered with 
brand 44-9 rubber, and also included corrugated mem- 
branes of AM-93 rubber, used in UAS instruments (Fig. 

1, c) 

The experimental investigation was to determine the 
following characteristics ; 

1) the dependence of the membrane's effective area 
on its travel , F, = f(S); 

2) the dependence of the membrane's effective area 
on the pressure drop across it, F_ = f (Ap); 

3) the membrane’s hysteres 

In addition, we determined the effect of the membrane! 
geometry, of the dimensions of its rigid center, and of 
the material's quality on the course of the characteristics 
y* f(S) and . = f (Ap), investigated the stability of 
membrane operation,and determined the possible skewing 
of the membrane which might arise during its operation. 

The experiments were carried out on stands specially 
developed for this purpose, in which atmospheric pressure 
acted on one side of the membrane and, on the other 
side, compressed air or a vacuum, the magnitudes of which 
were controlled. 

The forces arising from the pressure drops acting on 
the membranes were equalized by weights attached to the 
membranes’ rigid centers. The membrane characteristics 
were taken off for the following weight sizes ; 0.540, 
0.728, 4.62, 8.85, 13.09, and 17.32 kg. A cathetometer 
measured the membrane displacements. 
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For measuring skewing, pressure was introduced into 
a chamber over the membrane, with a long pointer 
attached to the membrane's rigid center on the load side. 
Membrane skewing was determined by the position of 
the end of the pointer fixed to a scale under the pointer. 
The scale was a uniform lattice of concentric circles 
and radial lines. Angles were read off from the radial 
lines, while the displacement of the pointer's end from 
the center was read off from the concentric circles. 

The magnitude of the membranes’ effective area 
was computed from the formula 

FP, — a ’ 

where F, is the magnitude of the effective area in cm?, 
G is the size of the weight placed at the membrane's 
rigid center, in kg, and Ap fs the pressure drop across the 
membrane, in kg /cm’, 

For the testing of the plane and conical membranes 
we used rigid centers of twosizes: d = 20 mm and 
d = 45 mm; the diameter D of the seals of these mem- 
branes was the same for all, and equal to 60 mm. A 
corrugated membrane (Fig. 1, c) was tested with a rigid 
center d = 45 mm and a value of D equal to 62 mm. 


EXPERIMENTAL RESULTS 


Stability of Membrane Characteristics 
and Membrane Skewing 

The characteristic P = f(S) , for new membranes, was 
only stabilized after the membrane had operated under 
pressure for several hours. The characteristics of one and 
the same membrane, in an unchanged set-up, may be 
repeated only in the case when they are taken off, one 
after another, without any interruption in membrane 
operation. Membrane characteristics taken off after an 
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interruption in operation are initially unstable, but then, 
after some time, are stabilized and begin to coincide with 
each other; but they almost never coincide exactly with 

the characteristics taken off before the interruption. The 
difference between these effective area characteristics 

can reach 10 to 12%. 

With a readjustment of the membrane, the character- 
istics taken off before and after the readjustment never 
coincide. As the pressure drop acting on the membrane 
becomes larger, the characteristics become more stable. 

In investigating membrane skewing, it was observed 
that they are unstable, but depend basically on the 
membrane’s positioning and, particularly, on the presence 
of nonconcentricity between the membrane’s rigid center 
and the seal's circle. 

Figure 2 shows part of the scale plus the graphs of 
the motion of the end of the pointer for two different 
positionings of one and the same membrane. The solid 
line shows the course of the pointer as the pressure in the 
chamber over the membrane was increased (forward travel), 
and the dashed lines show it for decreasing pressure 
(reverse travel). It was established that the skewing angle 
between the plane of the flapper, fixed on the membrane’s 
rigid center, and the plane of the nozzle's face end, 
which this flapper must cover, can oscillate, for different 
membranes, from 0.5° to 1.5°. 

Investigation of Plane Membranes 

A comparison of the characteristics F, = f (S) , taken 
off for membranes prepared of caprone, glass cloth and 
TD cloth, led to the following conclusions; 

a) the membranes of TD cloth possessed the greatest 
hysteresis; for them, the relative magnitude of the 
maximum hysteresis per effective area AF, can reach 
18% , while for membranes of glass cloth, having min- 
imum hysteresis, it was about 7 to 8 % (computing AF, 
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Fig. 2. Graphs characterizing the skewing of a membrane's rigid center. 
























as the maximum absolute value of the hysteresis divided 
by the maximum value of effective area F. max for 
forward travel of the membrane) 


b) for one and the same pressure drop on the mem- 
brane (i.e., for ome and the same value of F.) , the least 
travel was observed in glass cloth membranes; the travels 
of membranes of caprone and TD cloth were almost 
identical, and were larger than the travels of glass 
cloth membranes. 
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Fig. 3. Characteristics of plane membranes 
for d = 45 mm, D = 60 mm, G = 0.540 kg. 
Solid lines are for forward travel, dashed 
lines for reverse travel; @ is caprone, A is 
TD cloth, O fs glass cloth. 


A comparison of plane membrane characteristics 
(Figs. 3 and 4) with large (d = 45 mm) and small (d = 20 
mm) rigid centers leads to the following conclusions: 

a) for large rigid centers, the characteristics are 
nonlinear, while for small centers they are linear for 
almost the entire length of travel; 





b) with small centers, the hysteresis per area equals 
1.6 to 2%, i.e., fs much less than for membranes with 
large centers, in which it may oscillate between 7 and 
18 %4 


c) when the rigid center is decreased, the membrane 
travel is increased, but not by the same amount; in 
caprone membranes this increase is about 53 % of the 
travel for d = 45 mm and, for glass cloth membranes, 
about 23 %; 


d) as the diameter of the rigid center is varied from 
45 to 20 mm, or as its area varies from 16.2 to 3.2 cm’, 
the maximum effective area F, (with S = 0) is varied, for 
all the membranes, approximately from 21 to 12 cm’*, 


Figure 5 shows the characteristics y,. = f (S) for 
plane caprone membranes with a large rigid center 
(d = 45 mm), taken off for weights of different 
magnitudes. 


By considering the characteristics obtained, one can 
establish the following regularities in their behavior; 

1. The larger the weight, i.e., the greater the initial 
value of the drop necessary for the membrane to begin 
to move, the greater is the steepness of the characteristic 
Fe = f (S) as the pressure drop thereafter increases, f.e., 
the less does the effective area change over the course 
of the membrane's travel. 

This is explained by the circumstance that, with a 
large weight, there fs a large pressure drop across the 
membrane even before it begins to move, under the 
action of which pressure drop the membrane’s free surface 
is deformed, forming a semispherical crimp. With a 
further increase in pressure drop, a motion begins which 
is analogous to the motion of a crimped (corrugated) 
membrane. 

The formation of crimps in plane membranes 
occurring under the action of a pressure drop (thanks 
to which the variation of effective area as a function of 
travel decreases) increases the possibility of using plane 
membranes. Thus, for example, in the regulators of the 
English Sunvic company, instead of using doubled 
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Fig. 4. Characteristics of plane membranes for d = 20 mm, 
D = 60 mm, G = 0.540 kg. x is caprone, A is TD cloth and 
O is glass cloth. 
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Fig. 5. Plane caprone membrane characteristics, taken off for different weights. 


d =45 mm, D = 60 mm. 


corrugated membranes with atmospheric pressure between 
them, as is done in the UAS blocks, doubled plane 
membranes are used, with the pressure between them 
being obtained from the supply channel. Since the supply 
pressure is necessarily greater than the pressure in any of 
the regulator chambers, crimps are always formed under 
the action of this pressure. 


2 The greater the size of the weight, the less the 
hystersis per area AR, : if AF, = 3.6 % for a weight 
G = 0.728 kg, then AF, = 0.2% for G = 4.62 kg and, 
for G = 8.85 kg, it vanishes completely. 

The decrease in hysteresis observed for increase in 
weight is explained by this, that the internal frictional 
force which engenders the phenomenon of hysteresis 
becomes vanishingly small in comparison with the forces 
generated in the membrane from the pressure drop on it 
which is necessary to equalize the given weight. 


3, Figure 6 shows the characteristic F, = f (Ap) for 
a caprone membrane with d = 45 mm and D = 60 mm, 
constructed for different values of travel S® . It fs clear 
from the graphs that the effective area varies only for 
small values of pressure drop, after which, with increasing 
pressure drop, it remains virtually constant. 

This means that if the displacement of the mem- 
brane fs severely limited, Le., if it is virtually clamped 
into the position when it is initially established for values 
of travel of S = 0.4 to 0.8 mm then, as the pressure drop 
across the membrane varies within the limits of 0.1 to 
1 kg/cm’, the effective area will vary very insignifi- 
cantly. The explanation of this fact is very important, 
since it opens up the possibility of such an establishment 


of initial membrane position in devices for which 
deviations of effective membrane area should be 
minimal. 


Investigation of Membranes with Conic 





Profiles 


In the investigation of membranes with conic 
profiles, it was established that, with small weighs 
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Fig. 6. Plane caprone membrane characteris- 
tics. d = 45 mm, D = 60 mm. 





* The magnitude of $ is figured from the neutral 
position, in which the plane of the rigid center coincides 
with the plane of the membrane's seal. An increase in 
the travel S corresponds to a displacement of the 
membrane upward from the neutral position. 














(G = 0.540 g), it was necessary to take the membrane'’s 
natural rigidity into account in computing the effective 
area, since calculations with the formula F, = GAp give 
incorrect results. Indeed, in calculating with this formula 
for small values of the travel S, the magnitude of the 


for different weights and for different diameters of the 
rigid center, and with respect to hysteresis, turned out 
to be valid as well for membranes with conical profiles, 


Investigation of Crimped UAS Membranes 
For d = 45 mm and D = 62 mm, the characteristics 
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Fig. 7. Crimp formation in conical 
membranes under the action of a 


weight. 


effective area was obtained as equal to 30 to 35 cm? , 
while the magnitude of the piston area, whose diameter 
would equal the diameter of the membrane’s seal, 
equalled only 28.8 cm?, 

This is explained by the fact that, in the process of 
taking off the characteristics, the membrane was so es- 
tablished in its initial position that the rigid center's plane 
coincided with the seal's plane. With this , the membrane 
was deformed, and crimps were formed (Fig. 7). As a 
result, a force Q arose, resulting from the membrane's 
inherent rigidity, which counteracted the weight G applied 
to the membrane. 

An experiment was specially carried out to determine 
the inherent rigidity of conical and corrugated membranes, 
in which the free membrane was gradually loaded, and 
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Fig. 8. Characteristics of inherent mem- 
brane rigidity. 


the translation S of its rigid center under the action of the 
applied weight Q was recorded. The corresponding graphs 
are shown on Fig, 8. Curves C and D on this figure 
correspond to different samples of glass cloth conical 
membranes. 


It is clear from these graphs that, for small weights, 
the force Q becomes commensurable with the magnitude 


of the weight G , so that the computation of the effective 
area in this case can be carried out with the formula 
F, = (G-Q)/A p. 

All the fundamental regularities observed in the 
testing of plane membranes, both with respect to the 
behavior of the characteristics F, = f(S) and F, =f(Ap) 
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of these membranes are shown in Fig. 9, from which it is 
clear that the hysteresis and the variation of effective 
area as a function of travel are, for small weights, very 
large. 

However, for large weights (starting with G = 4.62 
kg), i-e., for large pressure drops, the amount of 
hysteresis sharply decreases, and there is a simultaneous 
change in the character of the behavior of characteristic 
F, = £(S). These characteristics, starting with a travel 
of about S = 0.45 mm, become practically linear and, 
what is particularly important, steep. 

The magnitude of a UAS membrane's effective area 
strongly depends on the position of the rigid center with 
respect to the plane of the membrane’s seal. This is 
particularly apparent for small weights. Figure (10) 
shows the characteristics of a UAS membrane for d = 45 
mm and D = 62 mm, taken off for weights of G = 0.540 
kg and G = 17.32 kg, and for various distances h of the 
membrane's rigid center from the seal plane. Charac- 
teristics A correspond to the neutral position, when 
h = 0; characteristics B correspond to h = -2 mm, Le., 
to a position of the membrane below the neutral position, 
and characteristics C correspond to the position when 
h = 2 mm, i.e., the membrane is above the neutral 
position. 

For UAS membranes we investigated the hysteresis 
as a function of the travel (Fig, 11) This confirmed 
the conclusion, given in [3], that the hysteresis 
decreases with decreasing travel. 








In Fig. 8, given earlier, the characteristics Q = f (S) 
were given, from which the inherent membrane rigidity 
can be determined. 


Characteristic A corresponds to the case when the 
UAS membrane is loaded on the side of crimp concavity, = 
and characteristic B to the case when the load is on the 
convex side. 


From a comparison of these characteristics it is 
clear that the membrane possesses different rigidity, 
depending on which side of the crimp the force is acting, 
specifically, the rigidity is greater in the case when the 
force acts on the crimp's concave side. 


SUMMARY 
From the experimental results, we may draw the 
following conclusions: 
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Fig. 9. Characteristics of UAS membranes, taken off for various weights. D = 62 
mm, d = 45 mm. 
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Fig. 10. UAS membrane characteristics, taken off for various positions of 
the UAS membrane with respect to the plane of the seal. 
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Fig. 11. Variation of the UAS membrane's hysteresis as a function 


of the travel. 


1. The characteristics of all the membranes can be 
stable only in the case when the membrane is always found 
under pressure. With large pressure drops, the 
characteristics are more stable than with small ones. 
Interruptions in membrane operation lead to a confusion 
of its characteristics which are then, after protracted 
operation, again stabilized , but frequently not to the 
same value which held prior to interruption of operation. 


Such membrane behavior can have a deleterious effect 
on the operation of the instrument containing the mem- 
brane and used for the automation of some productive 
process; for example, the adjustments may be disrupted. 


2. The hysteresis of plane meinoranes is less than 
that of corrugated ones. The magnitude of the hysteresis 
depends on the magnitude of the pressure drop across the 
membrane; the greater the drop, the less the hysteresis ; 
with large drops, hysteresis generally vanishes in plane 
membranes. Hysteresis is decreased with decreasing 
membrane rigid center. The greater the membrane travel, 
the greater its hysteresis. Among the plane membranes, 
those of TD cloth have the greatest hysteresis. 


3. The variation in effective membrane area along 
its travel, i.e., the behavior of characteristic F, = f (S), 
depends on the magnitude of the pressure drop across the 
membrane. For a given membrane travel, the greater the 
drop, the less the variation in the membrane's effective 


area. 






4. The effective membrane area varies as a function 
of the drop; the less the drop, the greater this variation. 
However, starting with a drop of about 0.15 kg/cm’, 
further increases in it lead to virtually no changes in 
effective area, f.e., it may be considered constant. This 
last statement is valid both for plane and for conical and 
corrugated membranes. 
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From whence the following important conclusion 
may be drawn ; if the membrane’s travel fs limited, then 
its effective area will be constant for all practical 


purposes for variations of the pretense drop on it within 
the limits of 0.15 to 1 kg/cm’. 
It is precisely in such conditions that membranes 
operate in instruments designed on the principle of force 
compensation. In these instruments, the pressure drop 
across the individual membranes might vary within the 
limits of 0.2 to 1 kg/ cm’, and the static translation 
of the membranes equals about 0.05 mm. The effective 
membrane area, during the instrument's operating 
processes, either remains completely unchanged or changes 
only very slightly, thanks to which an increased accuracy 
of such devices is obtained. 


5. The variations in effective area, both with travel 
and with pressure drop, are least of all in UAS membrane 
with semispherical crimps. 


6. The magnitude of the membrane’s effective area 
strongly depends on its initial establishment; this 
dependency is the more strongly pronounced, the less 
the pressure drop on the membrane. When the membrane’ 
rigid center is displaced downward with respect to the 
plane of the seal ( h< 0) the effective area increases in 
comparison with the case when the membrane fs in the 
neutral position (h = 0) and, when the displacement is 
upward (h > 0), it decreases. 




















In connection with this, the accuracy of the initial 
placement of the membranes in UAS instruments plays 
an important role, since by it one determines the 
deviation of the effective area of one membrane from 
another. In isodrome regulators, for example, this is 
inadmissible, since it might lead to the appearance of 
residual nonuniformities. 
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1. Experiment showed that the formula usually given 
in the literature for calculating effective membrane area, 


F, = zy, (D* + Dd + d?), 


which does not take into account the variation in area 
as a function of travel and pressure drop, can be used 
only for determining the value of F, for zero or very 
small travels of plane membranes, and only in the case 
when the membrane is established in the strictly neutral 
position (h = 0). 


8. Calculating membrane characteristics from 
experimental data by the formula F, = G/A p can be 
done only for large weights and, consequently, for large 
pressure drops. Calculations for small pressure drops 
give incorrect results, due to the fact that the 
membrane’s inherent rigidity has not been taken into 
account. 


The magnitude of the force arising from this rigidity 
is not a constant for different membranes, and may 
reach a magnitude of 300,. If this force is commensurable 
with the magnitude of the weight, it is then necessary 
to take it into account in the calculations, 





9 In certain devices of pneumoauto: atior as, for 
example, lead blocks, where there are not very stringent 
requirements on the invariance of membrane effective 
area, doubled corrugated membranes can be replaced 


by doubled plane membranes with supply pressure be-- 
tween them. 
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OPTIMAL FREQUENCY CONTROL 
OF ASYNCHRONOUS ELECTRIC MOTORS 
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(Leningrad) ; 
Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 333-339, March, 1960 
Original article submitted November 25, 1959 


The question is investigated as to which law should be followed in controlling voltage and frequency in 
asynchronous electric motors so as to provide optimal transient responses. The optimal control laws are derived 
for various requirements on the character of the transient response for the case when the impedance moment de- 
pends linearly on the rotational velocity *. 











For those stands where it is possible arbitrarily to vary ? ee 
the frequency of the voltage supplying the executive i}--- — 7 ’ 
asynchronous electric motors, the problem arises as to . 


which law the frequency should be varied by so as to 
provide optimal motor operation in the transient states. 
Electric drives can be divided into two basic groups as a 
function of the requirements imposed on the character of &>—- —-——%-—— -- 4 — ~ 
their transient responses. 

Electric drives of the first group must provide max- 
imum translation of the executive mechanism in a given | 
time for given values of initial and final rotational 

















velocities. Electric drives of the second group must 0 / Zin 
provide maximum change of rotational velocity of the Fig. 1. 
executive mechanism in a given time. The relative units will be : for motor rotational 


As examples of the first group of electric drives, one — torque, M = Mmo/ Mmax: for impedance moment, 
may take traction motors while, as examples of the second y = M,/ Mmax. for time t = t/T for rotor rotational 


group, one can adduce electric starting motors. speed, v = n/n., for rotor current Soquancy. w= frot/ for: 
Optimal control laws, as we shall see below, will be for current i = 1/1, and, for stator frequency, f = f,/fmay 
different for the different groups of electric drives. We make the following assumptions: 
We shall start from the basic equation of electric 1. We shall replace the true motor magnetization 
drive dynamics: curve a by the piecewise-linear curve b , consisting of 
two segments (Fig, 1). 
Mpo —M,=J = 9 (1) 2. We shall not take into account the electromag- 
netic transient response (the free component in the rotor 
current). 


where Mino is the motor's rotational torque, M, is the 


In relative units, the dynamics equation (1) takes 
impedance torque, depending on rotational velocity, and 


the following form: 


J is the moment of inertia about the motor shaft of the de (2) 
motor and drive. M—p=—z- 
We now introduce relative units. As the unit rota- The rotational torque M and the heat dissipation in 


tional torque we take the maximum (tilting) torque Mpax,» the rotor Q will be expressed by the formulas 
corresponding to critical slipping; as the unit of rotational 


speed we take n,, the synchronous velocity for the supply M = er , (3) 
line's maximum frequency; as the time unit we take Ty, D 

the mechanical time constant, numerically equal to the = _ te (4) 
time for accelerating the drive zero up to the synchronous : 1+ w? 

velocity for a constant rotational torque equal to the 

maximum; as the unit of current frequency in the rotor, * The basic propositions of this paper were presented 
we take fo,, the critical frequency corresponding to on Sept. 29, 1958 at the All-Union Conference on the 
tilting slipping; for the unit of current we take I,,, the Theory and Applications of Discrete Automatic Systems 
rotor current for critical slipping. (Section on Optimal System Theory). 
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where ® is the motor's magnetic flux in relative units ; 


@M _ mo 
Omax 





The yield of electric drives of the first group will 
be proportional to the integral 


t 


Aa = \ ide, (5) 


0 


and that of the second group to the integral 


Ay = ( dy = (va. (6) 
Ys 0 


In seeking an optimal control law, we should take 
into account the limitations imposed by the properties of 
the electric motors. The most important of these is the 
limitation of heating : the temperature of the winding 
must not exceed a limiting admissible value. Since the 
thermal transient responses in a motor flow much more 
slowly than the mechanical transient responses (heating 
time constants are measured in minutes [2]), one can 
ignore the heat emission during the start-up and braking 
times. With heat emission ignored, the heating limita- 
tion on the quantity of heat which can be generated in 
the rotor in a given time, i.e.. 

t 

2 42 
\ naa dt <Qy. (7) 
0 

With complete use of motor heating, (7) goes over 
into an equality. Steel saturation introduces a limita- 
tion on magnetic flux. With the assumptions we have 
made, this limitation, in relative units, is written in 


the form 0<O<1, (8) 


Finally, the motor's rotational speed is limited by 
the supply line's maximum frequency 


¥< Vmax (9) 
and, for =0, 


Ymax = 1. 


Optimal Control Law for Electric Drives 
of the First Group 

The optimal control problem for electric drives of 
the first group is formulated mathematically as follows ; 
to find w = w (T) and o*= er) which maximize (5) 
for a given value of (4), for given relational equations 
(2) and (3) and given limitations (8) and (9). This is the 
general Lagrange problem of variational calculus, com- 
plicated by the fact that only one-sided variation is 
admissible here. Consequently [1], the functions sought, 
Ww (T) and eT), must : 

1) either be extremal, i.e., satisfy the Euler equations 








oF  ¢0F _, oF _@ OF 
Ov 426dtayvOtC~ do dt dw’ 
OF d OF 
=0, 3or — ar 307 = 9 


for the auxiliary function F ; 
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M2, 

where Ag and \ are Lagrange multipliers, with A, 
constant and A a function of Tf ; 

2) or be boundary curves, f.e., satisfy (8) or (9) 
with the equality signs; 

3) or, finally, may consist of extremal segments 
plus segments of boundary curves. 

We turn now to the search for extremals. 


The first Euler equation leads to the expression 
- - _ & 11 
dh -— ho = a (11) 


from which, for p = Ht, + kv , we will Have 
ho 


A= yy — e"* . (12) 
and, for the most important case when pt = [lp = const, 
a ay C, — Kot. (13) 


The second Euler equation, after some transforma- 


tion, is written in the = 
@ 


ts 7* (14) 
We note particularly that (14) does not contain © 
and is valid for any magnetic flux (in addition to the 
trivial case @ = 0). 
To find the control law for magnetic flux, we use 
the third Euler equation. 
By differentiation of F with respect to owe get 
OF _ 2 2 


=. 15 
0m ~~ 1+ 0% a rere (*) 





By substituting the value of A from (14) in (15), 
we convince ourselves that 0 F/@ @" can not vanish and, 
consequently, B (1) which maximizes integral (5) can 
not be extremal, but is a boundary curve, Le., @=1. 

The optimal voltage control law is a control such 
that, with account taken of the voltage drop in the 
stator winding, the magnetic flux will be maintained 
constant. 

The optimal laws of current frequency variation in 
the rotor will be ; 

for # = const, 


2 


= = C.— Kot} (16) 


for fp = Hg + kv 
20 
i—o? 

For constant magnetic flux, the law of rotor fre- 
quency variation also determines the optimal law of 
variation of the motor's rotational torque as a function of 
time. By knowing the motor's torque and its mechanical 
time constant, we determine the rotor's speed of rotation 
and, thereby, the optimal law for supply frequency 
variation. 

This law will be valid until the rotor’s speed of 
rotation reaches its limiting value, bounded by (9) 
After reaching this limiting speed, one must transfer 


uN t 
=> Ca". (17) 
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Since the extremals must, at the junction point, coincide 
with the boundary curves tangent to them [1], then the 
equalitites M =p and dv/dr = 0 must hold for the 
transition points. 


Optimal Control Law for Electric Drives 
of the Second Group 

The problem of optimal control for electric drives 
of the second group is formulated mathematically as follows: 
to find w(t) and e*(r) which maximize (6) for a given value 
of (4), for given relational equations (2) and(3) and 
given limitation (8). Thus, the optimal control problem 
for electric drives of the second group is analogous to 
that for those of the first group, with the sole difference 


that the auxiliary function will now have the form 


29% ’ P 2a? 
Fi = Ta +2(y +e—7>u)- (18) 








The first Euler equation is now written in the form 
> _ ah 
Tae Oe ° (19) 
The second Euler equation will be 
1 
i * 
From the third Euler equation, with account taken of 
(8), we get, just as for electric drives of the first group, 
that the magnetic flux, for optimal control, must be 
maintained constant (®=1). For w = tg + kv we 
obtain, from (19) and (14), the expression for the optimal 
control law for the current frequency in the rotor: 
pas = Cee, (20) 
It follows from this that the rotational torque must 
increase as the drive starts up. For k = 0, Le., for 
Bt = const, it follows from (20) that w = const. Thus, 
we arrive at the important conclusion that, for constant 
impedance torque, the optimal control law for electric 
drives of the second group amounts to maintaining 
constant current frequency in the rotor, and consequently, 
constancy of the rotational torque. With this, the drives 
rotational speed will vary linearly and, consequently, 
the supply frequency must also vary as a linear function 
of time [6]. 


On the basis or the duality principle of variational 
calculus [1), all the optimal laws derived by us from the 


wnt 


assumptions that the quantity of heat was given, and 
that it was required to provide the greatest translation 
(or change of speed), will also remain valid for the case 
when the translation and the time are given, and it is 
required to provide minimum dissipation. The same laws 
also remain valid when the time of translation and the 
heat dissipation are given, and it is required to provide 
the minimum time of translation. 

The optimal control laws for asynchronous electric 
motors have great similarity to the optimal laws for dc 
motors [3-5]. If we compare the dependence of rotational 
torque on time for optimal control, then the optimal 
control law for de motors turns out to be the limiting 
case of the optimal law for asynchronous motors with the 
condition that the latter's maximum torque tends to 
infinity. In an optimally controlled asynchronous motor, 
the frequency in the rotor is always less than critical, and 
the rotational torque fs less than maximal. 

It should also be mentioned that the maximum 
torque of an optimally controlled asynchronous motor is 
larger than the maximum torque of a motor with a constant 
voltage on its terminal, since for critical slipping, due to 
the voltage drop in the stator winding, the motor's 
magnetic flux is less than the nominal value. With 
optimal control, when the magnetic flux is maintained 
constant, the maximum torque is 1.5 to 2 times greater 
than what is stated in the catalogs. 

Figure 2 shows an example of optimal control for 
a motor in which critical slipping equals 0.2, the 
maximum torque multiplicity (for complete magnetic 
flux) equals 5, the mechanical time constant is Ty = 1,82 
seconds, the impedance torque is constant and equal to 
the nominal value, and the electric drive falls in the 
first group. 

Determination of the Coefficients in the 





Extremal Equations 





We had found the extremal equations in a general 
form. For practical use of the optimal control laws 
found, it is necessary to know,in each concrete case, 
the coefficients in the extremal equations as functions 
of the given electric-drive parameters. 

We consider electric drives of the first group, with 
constant impedance torques. Extremal (16) is known, 
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from control by extremals (16) or (17) to control by 
boundary curves , Le., it is necessary to maintain 

Vv = Vmax and M = and, consequently, w = const. 

and we now derive useful auxiliary formulas. By carrying 
out a change of variables in (2), (3),and (4) on the 

basis of (16), and then integrating, we get 


2 1 1 
ial yal ~ Be =) 
2 rh) at 
Q= * Fe ~arc tara | . (22) 
We recall the self-evident formula 
v= f— o, (23) 


where f is the stator frequency in relative units and 
Sor is the critical slipping. 

By using these formulas one can, on the basis of the 
given parameters of the electric drive and the 
characteristics of the transient response, find the coef- 
ficients in the extremal equations. For the example 
given, let it be required that the electric motor with 
critical slipping s,, and impedance torque [p» start the 
drive , within time Ty, from initial speed vy =0 to 
speed UV, , corresponding to stator frequency f= 1. The 
dissipation during starting time Ty must not exceed Q, . 
The electric drive falls in the first group and, conse- 
quently, the optimal control problem {s to provide 
maximum translation in time T,. We know extremal 
(16). We determine the coefficients C, and A, which 
appear in (16), We start by determining rotor frequency 
Ws, corresponding to the end of the start-up period. 
From the equation 

203 
1+ w; 





= Po 


we find w, , after which, from (23), we find v, = 1=s,,W» 
Now , (21) and (22) give us a system of two equations 
for determining the two unknowns w, and Ag: 


= Eo 
“= [al 


@3 





= Yolo, 


Q = clits arc caro] , 


Knowing wy, we determine the coefficient 
Cy = 2wy/(1-wy") in (16). 


SUMMARY 


1, For the frequency control of asynchronous motors, 
there exists an optimal control law of frequency asa 
function of time which provides maximum yield of the 
electric drive for given limitations on heating. 

2. The optimal control law depends on the condi- 
tions imposed on the electric drive, and on the law of 
variation of the impedance torque. For electric drives 
of the first group, the optimal control law will be given 
by (17) and, for those in the second group, by (20). The 
voltage control law is the same for both groups, and 
amounts to maintaining the constancy of the motor's 
magnetic flux. 
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ON A DRIVE SCHEME WITH A GIVEN EQUATION OF MOTION* 


I, A. Boguslavskii 


(Moscow) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 340-343, March, 1960 


Original article submitted May 15, 1959 


A method is presented for the construction of a drive scheme with a given equation of motion wherein opera- 


tional amplifiers are not employed. 


In designing an automatic control system controlled 
by disturbances (cf.[1]) and also in various modelling 
schemes, it is frequently necessary to have a drive with 
a given equation of motion. This means that the drive's 
output coordinate x(t) (angular velocity, or angular 
displacement of the drive's output axis) must be related 
to the input coordinate f(t) (for example, with the 
arbitrarily varying angle of rotation of a potentiometer 
slider) by a given differential equation of the form 


n m ° 
d*x (t) 1 d*; (t) 
2 ae = ee (1) 
k=0 k=0 





The problem just formulated is easily solved under 
the following conditions ; 

1) if the drive is shunted by a total feedback path 
for angle or for angular velocity such that the servosystem 
thus obtained follows, with high static and dynamic ac- 
curacy, the output coordinate x(t) by some input voltage 
“ u(t) and, consequently, 


x(t) = ku (t), (2) 


where k is a scale factor; 

2) if, using operational amplifiers, one constructs an 
electronic scheme in which the output voltage u(t)** is 
related to the input f(t) by the equation obtained by 
replacing x(t) by ku(t) in (1). 

However, in certain conditions of industrial usage, 
an electronic scheme may turn out to be inapplicable 
due to its complex voltage supply, zero drift of the 
operational amplifiers, etc. It is therefore advantageous 
to consider the possibility of replacing an.electronic 
operational circuit by one composed of passive R and C 
elements (resistors and capacitors) which is free from the 
disadvantages mentioned above. However, a circuit of 
passive elements can be used directly only in solving 
certain special types of differential equations. For 
example, even for an equation which describes a simple 
oscillatory link, an RC circuit, as is well known, does not 
permit an accurate solution. In addition, the direct use 
of a circuit of passive elements (even in those cases 
where this is theoretically possible) is inconvenient in 
that a change in the coefficients of (1) requires, as a 
rule, a change in the parameters of all the passive 
elements. We shall show that a circuit with passive 
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elements can be freed of the disadvantages cited if a 
voltage proportional to x(t) is applied to certain of its 
links. In this case, the drive, together with the basic 
direct negative feedback , is shunted by several indirect 
feedback paths (Fig, 1), and thus participates directly in 
the solution of (1). The servo amplifier which directly 
controls drive operation does not fall into the category 
of operational amplifiers due, for example, to the 
completely different requirements on gain stability. 
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Fig. 1. 1 is the circuit of passive elements, 


2) is the drive. 





Fig. 2. 


For simplicity, we present the proposed principle 
for drive design in the case when (1) has the form 


(p® + ap* + a,p + ao) x (t) = 
= (bep® +- bp + bo) i (t), 


where the a; and bj are given constant coefficients and 
p= d/dt 

©The material in this paper is contained in the author's 
paper at the Second All-Union Conference on Automatic 
Control Theory. 


(3) 


* *Note that "u(t)" is used for both input and output voltage 
[translator's note]. 
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In the general case, the circuit design is carried out 
analogously. 

We consider the electric circuit (Fig. 2) in which the 
R, and Cy are constant resistors and capacitors, chosen 
from technical considerations, and the ej are certain 
voltage sources with zero internal impedance. In 
analyzing the circuit, for example, by the method of 
loop currents, we find that u , the voltage taken off 
capacitor Cj, satisfies the equation 


(p® +- App? + Ayp+ Ao)u = (B,p? + = (4) 
= Byp + Bo) e, + (Dip + Do) ee + Coes, 
where the Aq , B;, Dy, and Cy are constants, easily ex- 


pressed in terms of the circuit parameters, 
We require the holding of (2) and of the conditions 


et fhe (t) + hy f (5) 
¢: = Uhre (t) + kye/ (0), (6) 
¢s = + [hss (t) + kf ()), (1) 


where the k,, and kg; are certain positive or negative 
constants, defined successively by the relationships 


Ay — Bokxy ed ao, (8) 

A, — Byksxy —- Dykxe == aj, (9) 

Ay — Boku — Dokxe — Cokxs3 = 4, (10) 
Brky, = be, (11) 
Byky, +- Dyky, — bi, (12) 
Bok, + Doky, + Cots, = 50. (13) 


By a direct substitution of ey, €,, and e, in (4), we 
easily convince ourselves that, if (5)-(13) hold, x(t) and 
f(t) satisfy the original (3) and, consequently, the 
construction just described solves the problem posed. As 
Was previously stated, (2) holds as the result of construct- 


ing a high-quality system for foliowing x(t) by the voltage 
u(t). If x(t) is the angular velocity of the output shaft's 
rotation then, to realize (5), (6), and (7), with the shaft, 
in addition to the tachometer generator of the direct 
feedback path, there must be connected three dc 
tachometer generators whose output voltages are divided 
by special potentiometers in accordance with the coef- 
ficients k,; and then introduced into the proper places 

in the circuit (Fig. 3). In addition, to the axis rotating to 
the angle f(t) there must be attached the sliders of three 
potentiometers with insulated voltage supplies and out- 
put voltage dividers operating in correspondence with 
coefficients ks - The practical values of Ry, Rz,and R, 

in Fig. 3 are not less than the order of the large divider 
impedances shunted by the internal impedances of the 
tachometer generators. Therefore, Aj, Bj, and C; can 

be determined without the internal impedances of the 
voltage supplies being taken into account. 

If x(t) is the angular rotation of the output shaft, then 
to it, in addition to the potentiometer of the direct feed- 
back path, there are connected the sliders of three 
potentiometers with insulated voltage supplies and with 
the corresponding dividers. We emphasize that, for any 
change in the coefficients of (3), the parameters of the 
passive circuit (the quantities Ry and Cj ) are not changed. 
When there are new values of the k,4 and the kg , only 
the positions of the divider sliders are subject to change. 
In an analogous manner, by the use of a passive circuit 
composed of n series-connected RC networks, n + 1 
mutually insulated voltage supplies proportional to x (t), 

n insulated voltage sources proportional to f(t), and the 
corresponding number of dividers, one can construct a 
drive with equation of motion (1) forn > m. By having 
several circuits of passive elements and several drives 
with the corresponding number of insulated voltage sources 
proportional to their output coordinates,one can, by 
applying portions of these voltages to various points of 

the circuit, construct a system of drives with given 
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interrelationships of the equations of motion. The method 
described also applies, in an uncomplicated way, to the 
case when coefficients aj and b; of (1) are given functions 
of time, differentiable a definite numper of times. For 
example, in (3) let a = a(t), ay = a;(t), and a, = a(t). 





Then, coefficients k,; in (5)-(7) should be considered 
as unknown functions of time. As a simple verification, 
we show that x(t) will satisfy (1) for variable a, if the 
functions k,4 (t) are successively determined from the 
following relationships ; 


A, — Bokx, = a, (t), 
dk 


A, Lome Bykxy _ 2B, 


~ — Dykxg = 4, (t), 





x 
dt 
dk 





dk d*k . 
A, — Bokx; — B, > _ B= — Dokxs —D, 7 _ Cokxs = Ag (t). 


If we , by using the proper elements, move the 
sliders of the dividers in accordance with the laws k,,(t), 
kx,(t), and kxs (t) , we can implement the realization of 
a drive with equation of motion (3) for variable a9(t), 
a;(t) and a,(t). 

The static accuracy of the solution of (1) by the drive 
is determined by the stability of the coefficient k in (2), 
the stability of the values of the Rj and C, and by the 
accuracy of positioning the dividers’ sliders for the for- 
mation of coefficients k,; and kg. The dynamic 
accuracy of the solution of (1) by the drive will be the 
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higher, the wider is the frequency passband of the drive's 
servosystem, providing for the holding of (2), as compared 
with the bandwidth of the fundamental frequencies of the 
solution of Eq. (1) 
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ON THE NOISE IMMUNITY OF THE DISCRETE-INCREMENTAL METHOD 


OF INFORMATION TRANSMISSION 
V. M. Baikovskii 


(Moscow ) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 344-358, March , 1960 


Original article submitted August 28, 1959 


The potential noise immunity (noise stability) of discrete-incremental modulation (DIM) with fluctuation noise 
is determined. Comparisons are made of the noise immunity of pulse-code modulation (PCM), deltamodulation(DM) 
and DIM. It is shown that, in a number of cases, the noise immunity of DIM is higher than of PCM and DM. 


{NTRODUCTION 

In remote action telemetry systems, the communica- 
tions channel is the most essential element, and deter- 
mines the entire construction and, often, the accuracy 
of the system. This is true of telemetering apparatus in 
energy systems and gas and oil pipelines where it is 
necessary to transmit a series of parameters over distances 
of up to hundreds of kilometers. This is also charac- 
teristic of radiotelemetric systems, widely used in aviation 
for controlling various moving objects (including those 
used in space explorations), etc. 

In designing such telemetry systems, the form of 
signal and the type of modulation are selected, not only 
from considerations of the best construction of the trans- 
forming devices, but also as a function of the channel's 
noise level and other characteristics. 

Recently, many designers have been turning to 
discrete methods of transmitting data which, distinguished 
by their high equipment accuracy, permit high noise 
immunity to be attained. These methods acquire partic- 
ular value from the broad development of digital 
technology. 

Certain forms of discrete modulation had already 
been frequently considered as applied to communications 
problems. Pulse-code modulation is the first, and most 
highly developed, of the digital methods of transmission. 
There are many different PCM systems, differing only 
by the implementation of the individual elements of 
the apparatus ; transformers, accumulators, etc. But 
there is one thing common to all these systems; the 
complete instantaneous value of a parameter (with 
account being taken of quantization) is transmitted, in 
coded form, along the communication channel. 

There later appeared discrete modulation methods 
in which the signal in the channel corresponds only to 
the increment of the parameter,for example, DM and 
delta-PCM. 

Recently, still another method of this class for 
transmitting information was suggested—digital-incremen- 
tal modulation ("incremental coding*)[1, 2, 3]. 
Incremental coding arose in connection with the desire 


to increase the efficiency of transmitting data along a 
channel by exploiting the statistical properties of the in- 
formation to be transmitted. If the assumption in PCM 

is that each successive parameter value fs not related 

to the previous one, in DIM the design is based on the 
assumption that there is a comparatively peaked (at least, 
nonflat) correlation function for the information. 

The DIM operating principle is clearly elucidated 
in Fig. 1. The basic difference between this modulation 
method and those of PCM and DM is that here a signal 
is sent through the channel only at the moments when 
the measured quantity passes through a quantization 
level. If it stays within the limits of one quantization 
step, then no signal appears. Thus, transmission occurs 
irregularly; the more slowly the, information changes, 
the more rarely are pulses sent through the channel. 


Measurement limit 
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Fig. 1. a is the parameter to be transmitted 
b 1s the error signal in the communications 
channel,and c is the voltage at the output of 
the receiver's memory block. 
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In this paper we consider the noise immunity of DIM 
as compared with that of PCM and DM. As criteria for 
judging noise immunity of discrete telemetry systems, 
we adopt here the mathematical expectation of the error 
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(the mean error) and the dispersion of the error (the 
squared mean-square error). We have obviously made 
the assumption that the error's probability distribution 
function is sufficiently completely characterized by its 
first two moments. 


Potential DIM Noise Immunity for One 
Particular Case 

To estimate the effect of fluctuation noise on a discrete 
system for information transmission, it is convenient to 
use V. A. Kotel'nikov's theory of potential noise im- 
munity. It consists essentially of the following : 

Let the set of signals { A(t)} be used for transmission. 
When signal Aj;(t) is sent, then, due to the action of 
noise, the following total quantity is applied to the 
receiver: 








X (t) = Ay(t) + W (0), (1) 


where W(t) is a normal fluctuation oscillation. Depend- 
ing on the value of X(t), the receiver identifies it with 
(or assimilates it to) one of the possible signals Aj(t) of 

{ A(t)} with some probability of error. The potential 
noise immunity defines the probability of signal distortion 
as the result of the action of noise, for reception by an 
ideal receiver. 

An ideal receiver is characterized as one for which 
this probability of incorrect reproduction will be the least 
for a given form of signal (given energy) and a definite 
specific energy (per unit bandwidth) of the fluctuation 
noise. In other words, the ideal receiver reproduces the 
information, corresponding to that signal, which maximizes 
the expression 

P(At) Pg, (X). (2) 


Here, p(A;) is the probability of transmitting the 
ith signal and pa,(X) is the probability that the receiver 
reproduces the sum in (1) as the signal A, (t). 





and 
= 


A 
he, = & + xq in BA 1) 


P(Ag) ’ 
Thus, the potential noise immunity is defined by the 
two quantities :p(A;)/ p(A,) —the ratio of the a priori 
probabilities of the two signals—and the quantity a, which 
characterizes the ratio of the mean energy difference of 
these signals to the specific energy of the noise, i.e., 


Pp (A,) 


Qe = a+ -—lIn . 
” T 2a" p(Ai) 











+? 
2 
a = sa \(A@—4,OPat. ©) 
y 


2 


Here, T is the signal duration and o is the specific 
energy of the noise (intensity), which is defined so that 
co 
W(t) = (eo? (@) dw. (6) 
0 
Since the quantity to be transmitted in telemetry is 
always bounded by definite limits, then, on the average, 
over a sufficiently long interval of its variation, one 
can consider the appearances of pulses and pauses for 
DM and PCM to be equiprobable. 
Consequently, 
eee ee 
P(A,) = p(Aa) = in 
Ais = Koy 7 a. 
As a result, we obtain 
p(A,/As) = p(Ag/A,) = V (2). (8) 
Discrete-incremental modulation differs from the 
foregoing in that the signal can assume three values. In 
the channel appear the values +1 or ~1 when the meas- 
ured quantity goes through aquantization level and the 
value 0, ifthe quantity does not go beyond the bounds of 


In systems with DM or with PCM (with pulse-by-pulse © quantization step. Thus, there are the signals 


reception of the code combinations), there are only two 
possible discrete values of the signal in the channel, 
corresponding to a pulse and to a pause. This case is 
considered by Kotel'nikov [4] and in many other works 
[5, 6]. We therefore give only the necessary results here, 
omitting the detailed computations. For the case of only 
two signal values Ay and A», the probability of an error 
in transmission is expressed by the formula 


Pais = P(A) P(A2/A1) + P(A2) P(A1/A2) - (3) 


The quantities p(A;/ A») and p(A,/ Ay) characterize the 
probabilities that one signal transfers to the other as the 
result of the action of fluctuation noise. In correspon- 
dence with [4], they equal 


P(A2/Ay) = V (1), p(Ai/Az) = V (a4), (4) 





where ‘ co zt 
ad 24. 
V (x) ' Vin \< dz, 
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A,(t) = +1, Ag(t) = 0 amd A(t) = -1, with probabilities 
PC Ay), P(Ag) and p(As). With this, and for the same rea- 
sons as for(7), we can assume that p( Ay) = p(As) and, 
moreover, that p(A,) > p(A,;) = p(As), since only in this 
case is DIM meaningful. 


(+1) A,* PU) °A, (+1) 
(0) A, * A; (0) 
(-1) Ay? °As(-1) 


Fig. 2. Scheme of dis- 
crete channel (DIM). 


The scheme for a channel with discrete-incremental 
modulation is shown in Fig. 2. According to this scheme, 
the probability of error for DIM may be found from the 


formula 
3 3 


Per =D) Dy P(Ai) p(Ay/Ai) 


i=1j=1 


(i=*7). (9) 





or 
of 
de 


WU 





The values of the probabilities P(A A;) depend on 
which pulse signs are used for transmission. 


In its general form, the problem is difficult to solve 
even with three signals. It is impossible to find an ex- 
pression for the probability of distortion in a simple and 
graphic form, convenient for computations and for 
comparisons with the analogous expressions for other 
systems. We therefore consider one simple special case 
of signal transmission by DIM, and we then show that the 
results obtained can be used for an approximate estimate 
of noise immunity in a number of other cases. 


We assume that to signals A;(t) and A,t) there cor- 
respond identical video pulses of different polarity. Let 
their transmission be implemented by the so-called 
synchronous method, i.e., the pulses can be applied to 
the channel only at regular time intervals, or multiples 
thereof. Therefore, the intervals between pulses are 
multiples of some quantity Tin chosen, just as in DM, 
from the condition that there be reproduced the portion of 
the information with the greatest steepness (the overshoot). 
In this case, the analysis of noise immunity (cf. the 
Appendix) allows us to construct , for the ternary channel 
characterized by the scheme of Fig, 2), the matrix of 
transition probabilities ; 


A Ay As 
A, |1—V(a’) V(a’')[1—V(a")] V(a’)V a”) 
A; | V (a”) 1 — 2V (a”) V (a”) 
As | V(a’)V (a”) V(a’){1—V(e")] 1—V(a’) 


(10) 


To determine the mathematical expectation and 
the dispersion of the error, it is necessary to find the “value” 
or the weight of each individual error 4; The “value” 
of an error for discrete-incremental modulation is easily 
determined if we turn back to Fig. 1. 


When, under the action of noise, signal A, goes over 
to Ay , or Ag goes over to A», the step function voltage 
in the receiver will be changed by quantity A, ie., for 
(A, + Ay) and (A, + A,), € =+1 and, analogously, for 
(Ap Ag) and (Ay +A), E =—1. If (Ay Ay) oF 
(Ay As), then € = + 2 


Since it was accepted that signals Ay and A, are equi- 
probable, the system of signals is symmetric, and the 
average of the error equals Zero ; 


3 3 
MUer = > D P (Ai) (Aj/Ai) Fj = 0 


i=1j=1 


(i = /). 


(11) 
According to probability theory, 


DE = M (§ — ME)’. 
Thus, taking into account that ME = 0 , we find that 


DU oy = A*Di = A?MF. (12) 


By using (9) , we find that 
DU er = A?{2p( Az) V (@")4-2p (Ay) V (@’) [ —-V (a) + 


+8p(Ay)V (a) V (@’"). 


We now introduce the concept of the probability p,, 
of a transition of the quantity to be transmitted through a 
quantization level. Obviously, this quantity depends on 
the character of the variation of the parameter to be 
transmitted. The steeper the correlation function, the 
smaller is the value of p,- From what was said above, 
it follows that p, = 2p(Ay). By keeping in mind that 
P(A2) = 1=2p(A,), and by denoting p(A,)/ p(Ay) = K,we 
obtain the final expression for the dispersion of the error 
for DIM : 


DU, = A*pa{ KV (a") + V (a’) -|- 3V (a’) V (@”)].44) 


(13) 


It follows from (14) that the noise immunity of this 
type of modulation is connected with the statistical 
properties of the ensemble of possible communications. 
An analysis of this dependence was given in [2]. 


On the Approximate Estimate of Noise 
Immunity of DIM in the General Case 

It is obvious that, in addition to the one considered 
other methods of DIM transmission can be used. However, 
of all the possible methods,the only interesting ones from 
the practical point of view are those which are charac- 
terized by small values of the probabilities p(A;/ A;) 
and p(A;/ Ay) - 

The determination of an exact analytical expression 
for noise immunity for each case in even this limited 
class poses mathematical difficulties. 








In [4] there was given a method for the approximate 
estimation of noise immunity for the transmission of 
several discrete signals. It is based on the theorem that 
the probability P of at least one of m possible outcomes 
Bj lies within the limits 


P(Bimax <P < Yi p(Bi). 


k=l! 


(15) 


Here, the lower limit corresponds to the complete 
interrelationship of all the outcomes, i.e., to the case 
when one outcome implies the occurrence of all the 
others. To obtain the upper limit, it is necessary that the 
outcomes be incompatible. 


As applied to DIM, this inequality takes the form 


% 3 3 
p- p (Ai) p” (Aj/Ai) max < Pais S >> > p(Aj) p (A; Aj), 
i= imi j=l 
(16) 
where each p tAj/ A;)is determined, in accordance with 


(4), for each pair of signals without taking account of the 
others. 
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We considered several methods of transmission; for 
example, sending,for the signals A, (t) and A,(t), radio 
pulses with orthogonal frequencies, unipolar orthogonal 
video pulse, etc. Analysis showed that, for identical p 
pulse energy, the bounds on the error given by (16) for 
many methods of DIM transmission virtually coincided 
with those obtained in the case considered above. With 
this, the upper and lower bounds are so close that the 
region of error has the form of a very narrow ribbon 
(Fig. 3). By keeping this in mind, we can assume that 
the value of the error found for synchronous transmission 
of different polarity pulses characterizes, with adequate 
accuracy, the noise immunity of a number of other 
forms of DIM transmission. Therefore, in making our 
subsequent comparisons of noise immunity for the three 
forms of modulation, we shall use (14) to define the 
error with DIM (the dashed line on Fig. 3). 


U 2 3 @ 











IAG 
































k=10- 
0° 
wo 
P iis 
Fig. 3. 


It was previously noted that the noise immunity of 
DIM is determined on the basis of the probabilities of 
suppression or of formation of false signals A,(t) and 
A,(t) (since A, (t) = 0). 

If we therefore replace the “ternary” ideal receiver 
by two, one of which receives A;(t) and A,(t) while the 
other receives A,(t) and A,(t) , then the noise immunity 
will differ very little from the potential value. Com- 
putations show that the error arising with this is 
negligibly small. By taking this into account, we have 
the possibility, in making approximate estimates of 
noise immunity for DIM in many cases, to consider two 
identical binary channels separately. It is obvious that 
such an approach simplifies the selection of the most 
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efficient methods of transmitting and receiving, and 
eases system design. 


Comparison of Systems’ Noise Immunity 





For the comparisons, it is necessary to obtain ex-" 
pressions for the dispersion and mathematical expectation 
of the error for DM and PCM. With deltamodulation, the 
transition of one signal to the other in the channel (sup- 
pression of a pulse or the emergence of a false one) 
gives rise to a variation in the step function voltage in the 
receiver by the amount 2A. The error voltage with this 
can be expressed in the form Uer = 2A€ , where € fs a 
random variable which assumes the values + 1 and ~1. 
Since these values are equiprobable , the mathematical 
expectation of the error equals zero; 


MU,, = 24ME=0 (17) 
The dispersion of the error may be determined as 
DU, = 4A*D§; (18) 


since ME = 0, then DE = ME”. 
For deltamodulation, 


Me = 1 pA) +14 p(Av Ar). 9) 
By taking (8) into account, we obtain 
DE = V (a). (20) 
Thus, for DM 
DU, = 4A°V (a). (21) 


In the case of code modulation, the value of the 
error depends on which pulse in the code group will be 
distorted. The value of each individual pulse with 
ordinal number k in a code group of n pulses is found 
from the formula 


Ny, = 2*—-1A., (22) 
Therefore, the error voltage is determined as 


n 
Uer = >) 2*-1AE, (23) 
k=1 
Here the random variable, when noise acts, takes 
the values + 1 and —1 and, otherwise, the value zero. 
Exactly as in the previous case, we find that 
MU,, = 0 and 


n n 
DU ., = D >) 2k1AF = >) (2*)2A2DE, (24) 
k=1 k=1 
After some transformation, we find the formula for 
the dispersion of the error in the case of PCM: 
Y aealanned'. 


DU er = 3 





A*V (a). (25) 


A) The case of a given average transmitter power. 
We now consider the conditions under which the systems’ 
noise immunity should be compared. 

1, We assume that there {fs no filtering of the signal 
in the receiver , but simply its recording by a sufficiently 
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fast-acting reproducing device. We therefore carry out 
the comparison by using (14), (21), and (25). 

2. We assume that all the systems are designed for 
transmitting information with identical upper-frequency 
bounds Fy . 

3. We assume that the quality of transmission is 
identical for identical numbers of levels of partitioning, 
i.e., that the discreteness errors A are equal. 

4. We assume that the mean transmitter power P 
and noise level o are the same for all cases. 

5. We assume that the bandwidth and, consequently 
the duration of the pulses in the channel, is identical for 
all the systems. 

For simplicity, we shall assume that the pulses are 
rectangular. With this, the signal energy entering into 
Formula (6) will equal Ppt (here,P,, is the transmitter’s 
pulse power). 

Then , for the quantity « we obtain the expression 


1 Pp 1 By 
ed ed ee 


for code and deltamodulation, and the expression 


1 / eer 
an=ty/prnt)/ 7, (27) 


for discrete-incremental modulation. 

Here, the quantity T,,, is the average duration of the 
time interval between neighboring pulses. It depends on 
the statistics of the information to be transmitted and, 
obviously, is related to p, - Finding this relationship 
analytically in the general case entails mathematical 
difficulties, so that the magnitude of T,,, 1s determined 
from the corresponding distribution, taken off 
experimentally from recordings of information typical 
of the given telemetry system. 

The implementation of this investigation for several 
parameters subject to telemetering in industry, in 
energetics in particular [2], showed that the distribution 
of the intervals between moments when the parameter 
changed levels is close to a Poisson distribution. Indeed, 
if it be assumed that the transition probability Py, on the 
time segment At = Tmin depends on the length of this 
segment, i.e, Py = Pp At, and is unrelated to the behavior 
of the process on the remaining segments, then this 
distribution will be described by the expression 


H (t) = poe At, 





(28) 
and the magnitude of T, can be found from the 
formula on 

i 
Tay == ( t poe? dt aa ss. (29) 
0 
Since we have agreed that p, = PoTmin, then 
L 
Tay me (30) 
Pn 


In order to express the repetition period in terms of 
the quantity Fy for code modulation, we use the well- 
known relationship 


fi = 2Fy (31) 


The pulse repetition frequency, in the case of delta- 
modulation, is chosen, aswas mentioned earlier, from 
considerations of reproducing the maximum slope in the 
information to be transmitted | dS(t)/dt|,,4,* Just as in 
(7) and (8) , we use the simplified expression 


|A|= | ae fi = Oil. (32) 





Thus, the repetition period for DM, or Tmin for 
DIM, will equal 


A 


TipM = Tmin = ‘OnF 
u 


(33) 

On the basis of what has been said earlier, we can 
write the expressions defining the quantity a for all three 
forms of modulation; 


a A (34) 
7DM~ “5 VY 2F,V = 
tf a i 
* pem= ~~ |/ OF +A’ 
u 2 igs 
- (35) 


where IgJ (1+A)/A] = n is an integer; 


DIM = iV vV du + 2) , (36) 


Figure 4 shows the curves which characterize the 
potential noise immunity of the three forms of modula- 
tion considered, for two different values of A. It is clear 
from the figure that, for a given mean transmitter power, 
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the noise immunity of a system with discrete-incremental 
modulation is always better than in a system with delta- 
modulation, and only sometimes worse than the noise 
immunity of a system with PCM. With this, pulse-code 
modulation turns out to have higher noise immunity than 
DIM for large numbers of levels to be transmitted and 
high transmitter power. For small signal-to-noise ratios, 
conversely, the dispersion of the error is less when DIM 

is used. 

B) The case of given peak pulse power . We now 
compare the modulation forms we have been considering 
when the transmitter's peak power is given, retaining all 
the other conditions assumed under paragraph A. 

In this case, (26) and (27) are changed. If we denote 
by q the mark-space ratio of the pulses, then these 
expressions will have the form: 





a = Aj PP (37) 
PCM © oF 9 ] 1+A , 
“ul 8A 





{ io 
om. aig (38) 
As for the value of O&pj)4, we must consider two 
cases; 
a) when the duration of the pulses and the interval 
T are rigidly related to the quantity q ; 
) when the following conditions are imposed on the 
pulse duration; 


tT “#Y and t< Tin 
In the first case, the magnitude ofr will be the 
same as for DM. However, it should be taken into account 
that, if transmission is implemented by heteropolar video 
pulses then, for DIM, their amplitude must be half that 
for DM and PCM for identical dynamic ranges of trans- 
mission. Therefore, the formula for apypy now takes the 


form 
1 rv. 
mee i set we 
IM = Vw; ro (39) 


Now, if radio pulses with different carrier frequencies 
are sent through the channel, then the magnitude of api, 
is defined by Formula (38) , but the system's bandwidth 
is doubled. 

Figure 5 shows the comparison of the system's noise 
immunity for given peak power. It is clear from the 
figure that, under these conditions, DIM is inferior to the 
other forms of modulation except in the region of low 
signal-to-noise ratios. 

However, in practice (for example, in radiotelemetry) 
it is sometimes required that the conditions of "b” hold. 
In this case, limitations are imposed simultaneously on 
the transmitter's peak power and on its average power, 
but the pulse duration is not strictly given (except that 
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Fig. 5. 


T <= Tmin® With this, the magnitude of apy can be 
found from the formula 





_1,/ Pp A K+2) 
pim > s |/ oF Gx lg 





a (40) 


Thus, by increasing the duration of the pulses, one 
can improve the noise immunity of transmission and, 
simultaneously, reduce the required channel bandwidth. 
Broadening of the pulses is possible up to the point when 
the quantity (K + 2) /2q becomes equal to unity, after 
which the magnitude of apy Must remain constant 
The improvement of the noise immunity of DIM as a 
function of K for various admissible mark-space ratios is 
shown in Fig. 6. 
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SUMMARY 


1. The theory of potential noise immunity allows 
one to estimate the telemetering error with fluctuation 
noise for the discrete-incremental method of informa- 
tion transmission. An accurate and, at the same time, 
sufficiently simple expression for the probability of 
distortion with DIM could be found for only one special 
case synchronous transmission of heteropolar pulses. 
However, this expression can be used for an approximate 
estimate of the error in a number of other cases. 

2. If the signal-to-noise ratio is not very small and 
the probabilities p(A, / A,) and p(As/ Ay) are much less 
than the remaining transition probabilities, then a 
communications channel with DIM may be considered 
as twoindividual binary channels, and transmission and 
reception methods can be chosen in correspondence with 
this. 

3. Noise immunity of DIM is immediately related 
to the statistics of the messages to be transmitted; the 
greater the slope of the correlation function, the smaller 
the probability of error 

4. The comparison of the noise immunity of the 
discrete systems* showed that, for given average power, 
DIM is always better than DM, and is inferior to PCM 
only for high signal-to-noise ratios and large numbers of 
quantization levels. 

5. For given peak power and pulse duration, DIM 
has a practical noise immunity which is comparable to 
that of DM and significantly worse than for PCM. 

For a given mark-space ratio, a gain can be 
achieved only at the cost of widening the pulses. 


APPENDIX 

Computation of the Transition Probability Matrix 
for DIM with Synchronous Transmission of Heteropolar 
Pulses 

We determine the noise immunity for the given 
particular case of transmission with three discrete signals. 
For this we return to the basic propositions of the theory 
of potential noise immunity and to the definition of the 
ideal receiver which was given previously. We recall that, 
by Kotel'nikov's theorem, a signal which is limited to 
the time T and whose energy is lumped into some 
frequency band F may be completely characterized by 
a finite number 2FT of its values. These values are 
considered as the coordinates of the signal vector in a 
multidimensional space. The length of the vector is 
determined by the signal's average power, and equals 


a ae 
Ya=)V 40. 


r=] 











(A.1) 


Here,the a, are the signal values at 2FT points 
uniformly distributed in time. Analogously, the fluctua- 
tion noise is represented by a random vector of length 

he 2FT 


= am, (A.2) 
1 


aT & 


where 6 is a normal random variable for which 





p(z<0<2+d2)= 6 (A.3) 


dz 
V2n 

The terminus of the total vector, corresponding to 
the oscillation X(t) can, with some probability, fall at 
various points of space. Therefore, the problem of the 
receiver is the correct identification of the received 
signal X(t) with that one of the signals which was trans- 
mitted during time T. Thus, the receiver's design 
principle is characterized by some division or other of the 
entire signal space into regions of correct reception of 
each of the signals A, (t). 

In particular, for an ideal receiver, the configuration 
of the regions is so determined that it must reproduce the 
signal with the maximum value of the probability 
described by (2). With transmission of A,(t), the probabi- 
lity that signal X(t) falls in region dv of the 2FT-dimen- 
sional space equals 


dv 


Pa, (X) = (= 


T 
x exp(—-4 XO —4, OP). 








4 x 





(A.4) 


Therefore, the ideal signal, in receiving oscillation 
X(t), will reproduce that signal for which the quantity 





T 
oa (X(t) — A; PP — In p(A;) (A.5) 


will be minimal 

Thus, when signal A;(t) ofn possible signals is sent, 
the ideal receiver correctly reproduces the information 
corresponding to this signal if n-1 inequalities of the 
following form hold ; 


T [X (t)— A, (t)}* —o* In p(A) < 








(A.6) 


If at least one inequality does not hold, then distor- 
tion arises. Analogous conditions are valid for all n signals. 


There are the signals A(t) = + 1, A,(t) = 0 and 
A;(t)=-1. Signal A,(t) is found at the origin of coordi- 
nates, and the radius vectors corresponding to signals A,(t) 
and A;(t) are equal in modulus and oppositely directed. 
The boundaries between the signal regions are two parallel 
planes, perpendicular to these radius vectors and sym 
metrically located with respect to the origin of coordinates 
(Fig. 7). 


<T (X@—A, OF —o* in p(A)). 


* We recall that the comparisons made were valid for 
messages in which the distribution of the intervals be- 
tween pulses in the channel was close to a Poisson dis- 
tribution. 
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Re gion of l 
signal 





It is easily seen that, for distortion of signal A,(t), it 
suffices to violate just the one inequality 





T [X (t) — Ax (t)]? — o* In p(Ai) < 





<T [X (t) — Aa (t)}*? — 0? In p(A3), (A.7) 


since it is only possible to pass from the region of signal 
A,(t) to that of A;(t) by going through the region of 
signal A,(t). 

We consider the inequality inverse to that of (A.7). 
By using (1), we may write 


T [W (t) + Ai () — Ai (®)]}? — 0? In p(Ay) > 








>T (W(t) + Ai (®) — Aa (é)]® — 0? In p (Aq) (A-8) 


or 





TWO ST WO:+27TW iA —Aro) + 





+ T [Ai ()— 420)? — o? In Sent (A.9) 


By then taking into account that, in our case, 
A(t) = 0, we use (2.60) of [1]: 
WOAM=7aV4OH%  — (A. 10) 


Here, 9 is the normal random variable of (A.3). 
After some transformation we obtain 


0>0 V2T V A*() 14-7 A2(t) —o? n bi. (A-11) 


The probability that this inequality holds is V( @"), 
so that 
Paty (40 =V (2’)- (A-12) 


It is obvious here that 
V Tari) 0 p(Aa) 
a’ = — —_ In (A.13) 
oV2 ¥2V T A? (t) P (Ai) 
V TA) g LAs) 


— oo ae 





By substituting a = 





K, 


we find . 
a’ = 1 In kK. (A.14) 
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Distortion of signal A,(t) corresponds to the holding 
of the inequalities 


T [X (t)— Aa (t)}? — o* In p(Ay) > 








>T [X (t) — Ax (t)}® — o? In p(Aj), 





T [X (t)— Aa (t)}? — o? In p(Ay) > 





mT [X (t)— As (t)]® — 0? In p(As).  (A.15) 


It is easily shown that the probability of each ine- 
quality holding is the same for all, and equals 


p (As3/A2) = p(Ay/A2) = V (a"), (A.16) 


where { 
a”=a+>5-InK, (A-17) 


and the probability 
Psi, (As) = 2V (a”). (A-18) 


Due to the symmetry of the system of three signals, 
we have 


(As) = p 


A)). 
ais | 1) 


"Ais 
We can thus obtain the expression for the probabi- 
lity of error with discrete-incremental modulation in 
this particular case. 


Since 
3 
Pp = A. ph 2 Ae19 
er 2 P( dp. (A) (A-19) 
then 
P,. = 2p(A;) V (a’) + 2p(A2)V(a"). (A-20) 


er 


Equation (A.20) characterizes the probability of 
distortion as a whole. For us to compute the errors, it is 
necessary to know the transition probabilities of each 
signal to the two others separately. For this, we consider 
in more detail the transmission of signals A,(t) and As(t). 
Let signal A;(t) be transmitted. The probability that the 
ideal receiver reproduces precisely signal A,(t) at its 
output equals the probability of the simultaneous holding 
of the inequalities 





T 
p (A,) exp (— aIX@—Ai OP) < 





T 
<p(4:) exp(—3IXO—AOP), 





T 
p (A) exp(—3 IX —A wh) < 





T 
<p (As) exp (. st LX (t)—As OF) _ (Ae21) 


The probability that the first of these hold was 
already found earlier (A.12). 





$i 


hi 





We consider the second inequality. We rewrite it in 
the form 





T [X (t) — Ag (t)}* — o? In p (Ag) < 
T [x (t) — Ag ())? — o? In p (As). 





(A.22) 


After multiplying out the parentheses and cancelling, 
we have 


on Wan iP ain P (As) 
2T W (t) Ay (t) + T A? (t) — 0? In= (Ay) = 


«2T W (t) [Ay (t) — As (t)] + T [Ax (t) — As (t)]?. (Ae 23) 








By taking into account that |A,(t)| = |A,(t)| and 
p(Ay) = PCAs), we can write 


2T W (t) Ay (t) < 2T W (t) [Ai (t) — As (t)] + 





+3T A(t) +o? In K. (A.24) 


We now use Formula (A.10). Then, 


V2T V 4 (t) 6, <2 Y2T V 420) 815 1 








+ 3T A? (t) +07 1n K. (A.25) 


For the case considered, the quantities 6, and 6», 
are equal, since they are the projections of one and the 
same vector on one and the same line. With this in mind, 
we obtain the inequality 


0<2V2T V A? (t) 6, + 37 A? (t) + 0% In K.(A.26) 
The probability that this holds equals 
1—V(a”), (A.27) 


where 


1 
a” =3a+ 55 In K. (A.28) 


The probability that the two inequalities in (A.21) 
hold simultaneously equals the product of the partial 
probabilities, Le., 

P(A3/A1)=V(a’)[1—V(a")}. (A-29) 


The probability of reproducing signal A,(t) instead 
of A;(t) is defined by the inequalities 





T 
p(A) exp(— ZZ O—AOP) < 





ee 
<p (As) exp (— IX —As OP). 





y 
P (Az) exp (-3 [X (t) —AOP) < 





T 
Pp (As) exp ¢ ga [LX (t) —As OF) »  (A-30) 


Since the regions of signals A,(t) and A;(t) do not 
have common boundaries, the first inequality should be 
replaced by the inequality 





4 
p(Aexp(— 3 XO—A OP) < 





\ 


T 
<p(Asexp(—S(XQ@—AWP). (4.31) 


After transformations analogous to those given above, 
we find? that 


p (As/A;) = V (a’) V (a). (A.32) 


Since p(A,/ A) and p(A;/ As) have the same value, 
it is then easy, from the data obtained , to construct the 
matrix of transition probabilities (10). 
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T Generally speaking, this can be obtained from simpler 
considerations, by taking into account (A.29) and that 


Pris (Ay) = p(Ae/Ay) + p(As/Ai). 
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FINITE AUTOMATA. II* 


M. A. Aizerman, L. A, Gusev, L. I. Rozonoér, 


I, M. Smirnova, and A, A, Tal' 


(Moscow ) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 359-368, March, 1960 


In the second part of this work we continue the exposition of the authors’ point of view on finite automata, the 
problems in the theory of finite automata,and the interrelationships of this theory with the ordinary theory of switching 
(in particular, relay contact) circuits. The general problem is posed of constructing automata operating in one time 
scale (i.e., with one partitioning of the continuous time axis into discrete points) from elements which operate in a 
different time scale. As examples of this type of system, we consider abstract neurons (neurons of the type proposed 
by McCulloch and Pitts) and neuron nets. Relay contact and, more generally, switching circuits are then simply one 


possible technical realization of abstract neuron nets. 


6. Realizations of Automata Operating in 
a Given Time Scale from Automata Operat 
ing in a Different One 

In the first part of this work it was assumed that both 
an automaton as a whole and all of the elements of 
which the automaton was assembled operate in one and 
the same time scale (the term “sequentiality” was used 
previously) As applied to the last remarks of section 5, 
related to technical realization, this meant, for example, 
that all the delay elements simultaneously received sig- 
nals denoting the appearance of the next cycle (i.e., 
next moment of discrete time) from one and the same 
cycle transducer. 

Here we shall be interested in the more general 
problem of constructing automata which operate in one 
time scale from automata operating in a different one. 

We consider automaton A, defined by the equation 

x (p) = F [x(p— 1), p(p — 1)1, (@) 
for which there is given partitioning of the time axis into 
cycles t= 0,1,2,.., p,-.-. , defining the time scale T, 
and we construct a tape L of the states of this automaton 
which arise for some initial automaton state K(0) = Ky 
and for some sequence of inputs p = p(t) We transform 
this tape in the following way. From the sequence t = 0, 
1,2,.++,P,+++, we choose some increasing sequence 
of values of time ty, ty, tz, + +, ty, + + +, and we pick out 
from the original tape those columns p, K which corres- 
pond to the moments of time t, ty,...t; ,- +. The 
sequence of these columns, written in correspondence with 
the sequence of moments of time &, ty,-.+,tj,- ++, 
forms some new tape L® . Stillone more transformation 
may be carried out on tape L® by recoding all the sym- 
bols written in row"k*T . The recoding amounts to 
this, that all symbols from the {« } alphabet which 
appear on tape L® are divided into groups , and each 
group is put into correspondence with some new symbol 
(which , in particular, may again be a symbol from the 
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{« } alphabet). The symbols within any given group are 
not distinguished from each other, and one and the same 
new symbol corresponds to each of them. By replacing 
the symbols in column “K" of tape L* by the correspond- 
ing new symbols, we obtain tape L’. 

It may certainly happen that there are contradictory 
triads on tapes L® and L’, although there were none on 
tape L. It is important to note that even if there are 
contradictory triads on tape L", it is possible that, thanks 
to the recoding, there will be none on tape L'. If there 
are contradictory triads on tape L", then this tape cannot 
be the state tape of any finite automaton. But if L* 
contains no contradictory triads, then it could have been 
obtained as the result of the operation of some finite 
automaton, operating in time scale T', defined by the 
sequence of moments of time t, ty, - ++, tj, + +s 

We now assume that, for each possible state tape L 
of automaton A, it is possible to stipulate a sequence of 
moments of time t, tj,.--, tj, +--+ and a method of 
recoding + such that the transformed tape L' would not 
contain contradictory triads. With this, let there be no 
contradictory triads even among the set of all the triads 
contained in the entire collection of transformed tapes. 
Then this set of transformed tapes can be considered as the 


* For the first portion of this work, cf. Avtomatika { 
Telemekhnika No, 2 (1960). 

T The symbols in row "p” will not be recoded although, 
to be sure, this is possible. 

+ For different tapes it is admissible to choose dif- 
ferent sequences of time and different methods of coding 
(for example, the sequence of moments ty, ty,. ++, \» 

« « - may be determined by the sequence of inputs p(t) 
for which the given tape L was obtained; for an example 
of this type, see section 7 below). As a particular case, 
the sequence t, ty,.. +, 4, +--+. and the method of 
coding may be chosen independently of the appearance 
of any given tape. 
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set of state tapes of one and the same new finite auto- 
maton A‘, operating in the new time scale T'**® . The 
alphabet {x'} of states of automaton A’ fs defined by 
the accepted method of recoding and contains, as fs 
readily understandable, a number of symbols less than, or 
equal to, the number of symbols in alphabet {«}. The 
alphabet {p'} of inputs of automation A' is composed 
of all or some of the symbols of alphabet {p} . By num- 
bering the moments of time t = &, ty, -. - by the series 
of natural numbers 0, 1,...,i,... and by introducing 
the new time variable t’= 0, 1, 2,..., p",...+, wecan 
present the operation of the new automaton A‘ by an 
equation of the form 


x! (p’) = F’ [x (p' — 1), p' (p’ — 1)I, 
where p’=i1ifp=t- 


The problem of constructing a given automaton A‘, 
operating in time scale T', from elements operating in 
time scale T reduces to the realization, based on these 
elements, of an automaton A which, although operating 
itself in time scale T, may nonetheless be used, with 
some coding system and with the proper choice of a 
“new time", as the automaton A", operating in time 
scale T's Automaton A‘ may be represented as consisting 
of automaton A and a special “reading” device. This 
device records, for each state tape L of automaton A, 
the moments of time &, ty,..., Users, reads the symbols 
for the state and inputs of automaton A at these moments 
of time, and puts out the notation for these symbols in the 
adopted coding systemtt. In contradistinction to the 
cycle transducer which was discussed at the end of section 
5, the reading-coding device just described in no way 
acts on the elements of automaton A. 

Automata thus constructed, operating in time scale 
T', must themselves serve as elements for the construction 
of nets (section 5), mapping automata operating in the 
new time scale. 


The transition from automaton A to the new automaton 
A' was implemented by transforming discrete time t to 
t' and alphabet {x } to { K'}. As was stated earlier, 
because of the possibility of contradictory triads making 
an appearance, not each such transformation translates a 
finite automaton into another finite automaton. The 
problem arises as to what sort of transformation of discrete 
time, in conjunction with a recoding of symbols, does not 
take the device under consideration out of the class of 
finite automata. Closely related to this question is the 
question as to which class of device the aforementioned 
teading-coding device belongs to, and whether it is itself 
a finite automaton. 


We note in conclusion that the contents of sections 
5 and 6 reduce to the determination of different structures 
which are maps of one and the same automaton. What 
we are saying , essentially, is that different structures, 
thanks to recoding of states and transformations of time, 
may be mappings of equivalent automata. 


7. Two Methods of Constructing Automata 
by Transforming Time Scales 

Let there be a set of elements, consisting of arbitrary 
instantaneously acting binary logical transformers and of 
one-cycle delay element symbols from the alphabet 
{0,1}. We assume for simplicity that the cycles are 
defined by a partition of the continuous time axis into 
identical intervals of duration T . It is required to 
construct, from this set, an automaton forming a given 
automaton A! 


x (p’) = F [x(p’ — 1), p(p’ — 1)) (2) 


with alphabets {*} and {:p}, operating in a time scale 
defined by a partitioning vf the time axis into identical 
intervals of length 17 (/ O-= is an integer), i.e., into an 
i-fold decelerated time scale. 








It was shown in section 5 that a net, mapping any 
given automaton of any time scale T', can be cons-~ 
tructed from a set containing any logical transformers 
and from delay elements of the same time scale T’. We 
are now given logical transformers and delay elements in 
time scale T. The problem thus reduces to the cons- 
truction, from delay elements in time scale T, of delay 
elements in time scale T". 


u O Z; O a z z z; 


Fig. 1. 





** A time scale fs defined, not by giving a concrete 
partitioning of the time axis into equal or unequal 
intervals, but by giving a method which permits such 
a partitioning to be made in each concrete case. Thus, 
for example, the statement that each change of the 
input state means the beginning of a new cycle (Le., 
defines a new discrete moment of time) thereby also 
determines the time scale, since for each concrete giving 
of the sequence of inputs it permits one to carry out a 
definite partitioning of the time axis into cycles. 
Earlier, before there was any discussion of trans- 


formations of time scales, this remark would have been 
meaningless. 


TT Since, for different tapes of automaton A (Le., for 
different inputs and states of the automaton), different 
sequences of moments of time ty, ty,.. «ty +--+ and 
different coding methods are permitted, then this 
device must, generally speaking, obtain signals as to the 
inputs and states of automaton A at each moment of 
time t= 0, 1, 2,44, Drees 


This device is intended, essentially, for the represen- 
tation of events. Specifically, it represents some event 
characterizing the operation of automaton A (for example, 
each change of state of A's input , the appearance of 
some state in automaton A itself, etc.) by putting out 
a reading signal at the moment when the event transpires. 
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For this purpose , we join / delay elements, operating 
in time scale T, in series so as to form a delay line 
(Fig. 1), which is described by the recursion formulas 


x, (p) = 2, (p— 1), 
Zo(p)= 23(p— 1), 


a, (p) = u(p— 1). 


By eliminating from these relationships all the x; 
except xy, we obtain, for the delay line as a whole, 


2, (p) = u(p—?). 

The delay line's output at time s/ is uniquely 
determined by the input at time (s-1)/ . Therefore, if we 
set t' = s for t = sl , then such a delay line may be consid- 
ered as a delay element operating in time scale T’. 

We now construct an automaton by the assembly metha 
(section 5) but now, instead of delay elements in time 
scale T, we connect, in each instance, a delay line as 
just described. The automaton thus constructed will 
operate in time scale T. But if we read the input symbols 
of this automaton, and the output of each delay line, 
only every IT seconds, then we shall obtain the required 
automaton, operating in the time scale we require. 

The reading-coding device of which we spoke 
earlier can, in the given case, be a time counter , which 
gives a reading signal for the automaton inputs and the 
delay lines’ outputs every! T seconds; here, coding is 
implemented automatically since only the delay lines’ 
outputs are read, by virtue of which those states of 
automaton A which are not differentiated by the states of 
the terminal elements of all the delay lines are therewith 
lumped into one group. 

We consider me particular case of the problem of 
realizing an automaton, operating in a given time scale, 
by means of an automaton operating in another time 
scale. This case is of interest in its own right, since it 
plays a basic role in the understanding of certain important 
technological methods of realizing finite automata.in 
general, and switching circuits in particular. 

We consider some automaton A, given by (1). At 
time t), let the input assume the state denoted by symboi 
p’, and let the input state thereafter remain unchanged 
until time ty. During the entire time t) = t= ty, the 
automaton operates as an autonomous automaton, and its 
changes of state occur in correspondence with the equation 


x(p) = F [x(p— 1), p*), (3) 


We now assume that, for each p , there exists a 
numbers such that; 

a) this autonomous automaton is decremented, Le., 
its graph fs such that, fron any initial state, it can reach 
one of several possible equilibrium states [i.e., a state 
which, by virtue of (3), is transformed into itself] within 
a finite numberh (h=s) of cycles, and 

b) the following change of the automaton's input 


by< pk. 
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occurs only after the automaton has reached an equilibrium 
state, Le., ty~te = s. 

At time t = ty, input state p° ts changed to p*, which 
then remains unchanged until time t= t,. There is then 
formed a new autonomous automaton 

x(p) = F [x(p — 4), o'], h< Pty, 
about which we make the same two assumptions. Once 
again, after a time less than t,-ty, the automaton reaches 
an equilibrium state, which may differ from the previous 
one. At time t = ty, input p! is changed to p’, etc. 

We assume that the automaton's inputs and states are 
read only at the moments to, ty, ty, - , when 
the input states change. It is not necessary here to recode 
the states. 

Tape L", written by the corresponding “reading” 
device, 


cor Gees 























does not contain contradictory triads, as is easily seen. 


Therefore, tape L' can be considered, for any input 
sequence p°, Pp, p’, - ++, as the result of the operation 
of the new automaton A’, operating in time scale T’ 
which is defined by the moments at which the inputs 
change. Serving as the states of automaton A‘ are the 
equilibrium states of the autonomous automata formed 
from automaton A for fixed (sustained) inputs. 


It is immediately evident that the operation of 
automaton A‘ is completely independent of the time 
scale of automaton A. It is important only that automaton 
A operate sufficiently rapidly, i.e., that equilibrium is 
attained in it before a new cycle commences in automaton 


8. Abstract Neurons and Neuron Nets 

The name “abstract neuron” is given to an element 
with a finite number of input lines and one output line, 
the states of which can be denoted by symbols from the 
alphabet {0, 1}; the symbol at the output is uniquely 
defined by the symbols at the inputst seconds previously, 
where T is a given number# + 





Thus, a neuron is a combination of logical trans- 
formers and delay elements with lag time tT. The 
specific feature of neurons is that the logical transformers 
and delay elements cannot be distinguished from each 
other. 


+ The term “neuron” for these abstract elements fs 
used for historical reasons, and is only slightly related to 
the physfological meaning of the term. 








on 





A neuron may be considered as an elementary 
finite automaton in which only two states are possible 
whereby, with a fixed input, an equilibrium state is 
reached in one cycle (the time between cycles is T 
seconds). 


Neurons are connected to each other, forming 
neuron nets. For the formation of a net, the states of 
several or all input lines of the neurons from some set 
are immediately identified with the output states of 
some neurons from the same set. In contradistinction 
to the nets treated in section 5, neuron nets do not 
contain symbol transformers, and the original elements 
are not any automata, but automata of a very special 
form, namely, neurons. 


The logical functions which can be realized by one 
neuron can be most diverse. But in neurons of a given 
type only a completely concrete set of logical functions 
may be realized. By setting up a net from neurons of a 
given type, this net having one output line and several 
input lines, one can extend the set of realizable logical 
functions, in the sense that the state at the output of such 
a net will be a logical function of the input states, but 
with a delay , not of T, but of lt seconds, where / is an 
integer. We now assume that, for any logical function 
whatever, one can, from neurons of a given type, cons- 
truct a net which realizes this function in not more than 
tr seconds. The least value of € which satisfies this 
condition will be called the index of the neuron. The 
number € depends on the logical capabilities of the 
neurons from which the net was set up. *** 


It follows from what has been presented that, from 
neurons with a finite index € , one can construct any 
logical transformer but, in contradistinction to the 
transformers studied in previous sections, these will not 
be instantaneously acting: up to € T seconds will be 
required to implement a logical transformation. 

A neuron net is a particular case of a net in the 
sense of section 5, so that any neuron net is therefore 
a finite automaton. The time scale of this finite 
automaton is determined by the partitioning of the 
continuous time axis into equal intervals of duration T. 
Obviously , if one is considering the operation of an 
automaton in such a time scale, then neurons with index 
¢> 1 do not constitute a complete set. Indeed, in such 
a time scale it is impossible to construct from neurons 
even the simplest automaton. 


x(p) = L[x(p— 1), u(p— 4), 


where L is a logical function for the realization of which 
more than one cycle is required (if €> 1 then, by defi- 
nition, such a function L exists). 

Now, if the time scale is defined by a partitioning 
of the time axis into identical intervals equal or greater 
than ¢7, then any automaton may be constructed from 
neurons with index ¢ . In fact, any logical function is 


realized within ¢ + seconds, as well as (by the use of feed- 
back paths) a system of equations of the form 


x4 (p’) = Li [x1 (p’ — 1), «209 %n (p’ — 1); 
u,(p’—1),...,Us(p’—1)] (i= t4,..., 2), 


where p'=s if p=s¢, i.e., an abstract structure which 
is the mapping, as was shown in section 5, of any finite 
automaton of time scale-p', 

By using automata operating in such a “uniform” time 
scale ¢ T, we can construct, from neurons, automata 
operating in other time scales. On the basis of the 
arguments given in sections 5-7, we can now give at 
least two methods for realizing such automata. 
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First method. From neurons with finite index ¢ 
we can construct delay elements for the time defined by 
the externally supplied special cycle signals u,- Figure 2 
shows an example of a circuit of such elements, with 
input u and output x . The transformers L shown on the 
circuit realize the following logical function with a 
delay! r (where 1 =< €): 


z= (a&c) V (b&c). 
In correspondence with this, the equations which 
describe the delay elements have the form 


a(p) =ly(p—) &u(p—)IV [z(p—) &u(p—D)I, 
y (p) =(u(p—) &u(p—)IV ty(p—) &u(pP—YI, 


Correct operation of the scheme shown can be 
guaranteed in the case when the new cycle signal u; 
lastsfor / old cycles and, moreover, if input u does not 
change before signal u, is applied. Figure 3 shows the 
changes in time of allthe net coordinates for certain 
given changes inu andu,- Figure 4 shows an example 
of a McCulloch-Pitts neuron net which realizes a delay 
element for a time defined by synchronizing input u, - 

From such delay elements and from logical trans- 
formers which, as was shown above, are constructed from 
the same neurons, one may construct any finite automaton, 
operating at a tempo defined by synchronizing input u, - 





*** It is certainly possible that ¢ =o , or that there 
exist logical functions which, in general, cannot be 
realized by neurons of the given type. 

The McCulloch-Pitts neuron has index 2, since no 
more than 27 seconds is required to realize any complete 
normal form in a net constructed of McCulloch-Pitts 
neurons. 
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Second method . If it is necessary to construct an 
automaton whose time scale is defined by the changes in 
the input states, then one may use, in neuron circuits, the 
method of constructing automata from the equilibrium 
states of the corresponding autonomous automata (cf. 
section 7). Certainly, the first method is the more 
general; it is also valid in the given case, when the 
signals u,are engendered by each change of input states 
in the automaton. 





9. Relay Contact Circuits as Finite Auto- 
mata 

In the consideration of relay contact circuits there 
arise many important technological and scientific pro- 
blems related to the peculiar features of this branch of 


technology. The problems alluded to will be briefly 
stated at the end of the present section. We, however, 


will not be interested here in these problems, however 
important they may be for the construction of actual 


relay contact circuits. 
Here, we consider relay contact circuits only from 


the point of view of their special features which permit 
their assimilation to finite automata and which permit 
the use of the methods of the theory of finite automata 
in the construction of relay contact circuits. Therefore, 
we shall henceforth in this section consider an idealized 
model of a relay circuit. 

By a relay we understand a device with a finite 
number of input lines and one output TTT, and which 
satisfies the following two requirements : 








a) each of the input lines and the output may have 
only two states, denoted by zero or unity; 

b) the output state is some logical function of the 
states of the input lines with some time lag fT . 
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For example, in ordinary electromechanical relays 
with one winding, each relay has one input (the presence 
or absence of current in the winding) and realizes a lo- 
gical function of one variable of the type " yes” (nor- 
mally open relay contacts) or “no” (normally closed 
relay contacts) The lag of T seconds is due to the relay’ 
operating time. 

In using relays with several windings, each relay has 
several inputs, and the character of the logical functions 
realized (with a delay of T seconds) in the relays is 
determined by the special characteristics of the windings 

Individual relays are combined in a relay circuit by 
means of electrical circuits (contact circuits) These 
circuits define all or several inputs to relays as logical 
functions of the states of relay contacts and of contacts 
whose states depend on external stimuli. These logical 
functions are realized virtually instantaneously in practice 
Therefore, the contact circuits play the role of the instan- 
taneously acting logical transformers. 

From what has been presented it follows that a relay 
is an ordinary neuron, the value of whose index ¢ depend 
on which logic can be realized by means of one relay. 
A relay circuit differs from a neuron net only in its 
capability of being used as an instantaneously acting 
transformer. This difference is inessential, since the 
essential fact for a neuron net is that its time scale, as 
an automaton, is predetermined by the time T— the 
neurons’ lag. 

Therefore, relay circuits are finite automata of the 
“neuron net" type, and everything presented in sections 
6 and 7 applies to relay circuits as well. In particular, 





TTT Or two outputs, one of which is the negation of the 
other. 
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the two methods given in section 7 for constructing 
finite automata which operate in time scales different 
from the neurons’ natural time scale, also define two 
possible ways of constructing relay circuits which map an 
arbitrary finite automaton whose time scale differs from 
the natural time scale of the relay circuit, as defined by 
the relays’ time lag. Both these methods are used for 
synthesizing relay circuits. The first of these (with the 
construction of delay elements) is used most frequently 
in the construction of the so-called cyclical circuits. 
However, it is clear from section 7 that this method is 

a general one, valid for the relay realization of any 
finite automaton. As an example, Fig. 5 shows the 
schematic realization of a delay element for the time 
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Fig. 5. 


defined by synchronizing signal u_, constructed from 
ordinary electromechanical relays. The second method, 
using the equilibrium states of the corresponding auto- 
nomous automata, is actually used in almost all 
sequential circuits. 

Although relay circuits are one of the possible me- 
thods of technically realizing finite automata, specific 
analytical and synthesizing techniques are applied to 
them. This specificity is due, in particular, to the fact 
that the lag times in different relays of one circuit are 
not identical. Therefore, the time axis cannot be 
partitioned beforehand into equal cycles, thanks to which 
the analysis and synthesis of relay circuits requires that 
one obtain additional information which is not forth- 
coming from finite automata theory. Moreover, the 
differences in relay operating times might lead to relay 
competitions. 

Also specific are the minimization problems which 
arise in the synthesis of relay contact circuits (for 
example, minimization of the number of contacts, 
windings, etc Thus, relay circuits pose their own 
problems which do not arise in the general theory of 
finite automata. To our way of thinking, the aim of 
telay circuit theory consists of studying such types of 
specific problems and of general methods for the tech- 
nical realization of finite automata by relay circuits of 
some type or other. 

Since finite automata theory began to be crystal- 
lized as an independent scientific discipline only very 
recently, it is natural that part of the general questions 
which are important for the understanding of finite 
automata in general were developed within the purview 
of relay circuit theory. Due to this, there is ordinarily 


no distinction made between the general questions bearing 
on finite automata theory and the questions specific to 
relay realizations. 

The development of the theory of logical devices 
naturally led to the introduction into consideration of 
abstract “finite automata” and to the splitting off of 
finite automata theory as an independent scientific 
discipline. 


10. Concluding Remarks 





In designing concrete automatic devices, it is 
necessary to distinguish several independent steps. 

Design begins with the analysis and idealization of 
the operations which the automatic device must execute. 
If this idealization leads to a formulation of the operating 
conditions in terms of discrete time and a finite set of 
variables, each of which may assume only a finite 
number of values, then the possibility arises of seeking 
the automatic device to be designed in the class of 
finite automata. With this possibility there is closely 
related a problem subsumed under the general problem 
of representing events, and the eluciaation of what may 
and may not be implemented by finite automata. 

If this question is answered in the affirmative, there 
arises the second problem, namely, by which equation of 
type (1) may one describe the finite automaton which 
satisfies the conditions posed , i.e., what are the tapes 
of this automaton and of what triads do they consist. 
This problem relates to the problem of the abstract 
synthesis of a finite automaton from its given operating 
conditions. 

Even within boundaries of abstract consideration, one 
may cite a number of ways of realizing automata. One 
and the same automaton may be made, for example, 
from elements implementing functional transformations 
of two-valued or of multi-valued logic, with the use of 
special delay devices for externally given cycles or without 
their use if the equilibrium states in abstract neuron nets 
are employed, etc. The choice among the theoretically 
different ways of realizing an automation is related to the 
choice of the automaton's structure, to its abstract 
assembly. 

This entire group of questions comprises the contents 
of the general theory of finite automata. 


Each technical realization of a finite automaton 
gives rise to its own scientific problems, requires its 
own methods of analysis and synthesis, the taking into 
account of the technical special features of the actual 
apparatus, the economic and other requirements on the 
design or the circuit. 

The application of finite automata theory ends as 
soon as one has chosen the structure or the net, L.e., as 
soon as one has written the equations of type (6) of the 
first part of the present work . 

For the problem of the technical realization, the 
equations of type (6), or any of their equivalents, are the 
initial data, and the aim here {s the engineering choice of 
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design or circuit which is optimal in some sense (tech- 


nical or economical) 


Certainly, the boundaries between the problems of 


finite automata theory and the problems of technical 


realization are not well-defined. 
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Even in the choice of an abstract structure, one 
has to take into account the requirements and the 
special features of the technical realization. But this 
“infiltration” of technical requirements into the abstract 
considerations does not impede the detachment of the 
general theory from the very important, but specific, 
problems. 








A VISUAL-MATRIX METHOD 
FOR MINIMIZING BOOLEAN FUNCTIONS 


A, D, Zakrevskii 
(Tomsk) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 369-373, March, 1960 


Original article submitted April 25, 1959 


A semiintuitive method is presented for the minimization of Boolean functions which is based on the use of 
certain properties of visually perceived information. The method {fs applicable to the minimization of arbitrary 
Boolean functions with up to 10 to 12 variables, and is particularly effective for the minimization of incompletely 


defined functions. 

One of the most important problems which must be 
solved in the synthesis of switching circuits— the 
minimization of the number of circuit elements ~ is ordi- 
narily solved by transforming the algebraic form of 
Boolean functions until one has arrived at formulas which 
are superpositions of minimal numbers of elementary 
functions. 

The difficulties which arise in minimizing arbitrary 
Boolean functions increase catastrophically fast as the 
number n of arguments increases, so that a number of 
methods become inapplicable in practice even for n>6. 

The increase in the number of arguments of the 
function to be minimized may be compensated for either 
by restricting the class of functions to be considered or by 
increasing the speed of solution of the minimization 
problem by using modern computing technology. The 
realization of these possibilities is the main subject of a 
large number of published works [1-17]. 

A common characteristic feature of the contents of 
a large part of these works is the algorithmization of the 
process of minimizing Boolean functions. Significantly 
less attention has been given to the development of 
methods based on the use of certain properties of the 
human perception of information which are difficult to 
translate into algorithms. Despite this, just such methods 
are necessary for manual minimization of Boolean 
functions. 

The original information on Boolean functions must 
be presented in compact and readily perceivable form. 
Since the basic “memory device” for manual minimi- 
zation is the two-dimensional sheet of paper, the most 
compact form of representing information in this case 
will be a two-dimensional—matrix—form. Matrix forms 
for presenting Boolean functions were used by O. Plekhl' 
[10], A. Svoboda [18], M. Karnaugh [19] and others. 

In these methods, the elements of information 
concerning the Boolean functions are the functions’ values 
given by sets of elements in two-dimensional space, 
placed into a single-valued correspondence with the sets 
of values of the variables. The order of the mutual 
spatial arrangement of the elements, a matter of in- 
difference in an algorithmic approach, plays an es- 
sential role for the visual perception of the matrix by a 


human being. An advantageous choice of such a mutual 
spatial disposition of elements, i.e., such a method of 
coding the spatial elements by sets of values of the var- 
fables, is one for which the properties of the Boolean 
functions which are of interest to us in trying to carry out 
minimization, would be readily eye-catching. 

In the method suggested by A. Svoboda for mini- 
mizing Boolean functions, the potential possibilities for 
minimizing a Boolean function, as represented by a 
matrix, are expressed by the periodicity properties of that 
representation. With large numbers of variables, one has 
to do with the superposition of periodic representations 
with different periods which is difficult to preceive visually 
as a whole. In connection with this, A. Svoboda uses 
lattices (grids), corresponding to the individual variables, 
by means of which one successively {isolates the portions 
of the matrix which correspond to definite values of 
these variables. 

In the author's opinion, it is easler to perceive 
visually the symmetry properties of the representation, 
which gives the possibility of using the superiority of 
Gestalt perception of the original information, thus 
dispensing with cards and leading to the successive 
analysis of the individual portions of the representation. 

The visual matrix method of minimizing Boolean 
functions which the author proposes is sermialgorithmic, 
semiintuitive, so that its practical implementation is 
required in order to elucidate its superiorities over ocher 
methods. 

It should be remembered that, for a large number of 
arguments, due to the complexity of the problem, the 
minimization of a Boolean function in the general case fs 
considered to be successfully implemented if it has led 
to a sufficiently minimal algebraic form (i.e., strictly 
speaking , to a form which is sufficiently close to the 
absolutely minimal one). 

An aribitrary Boolean function f(x , X,, . - +x,) can 
be given by a matrix ( a4),) whose binary elements a;), 
are the function values, while the subscripts 1, k are it 
arguments. The pairs of values of i, k are put into a 
one-to-one correspondence with the sets of values of the 
arguments xy, Xp, + - «X,, Whereby the arguments are, for 
convenience, divided into two sets of as close to equal 
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size as possible. Values of the arguments which enter 

into one of these parts code the values of i , the values of 
the remaining arguments code the values of k . O. Plekhl’ 
and A. Svoboda coded the natural number i,k by a 
positional binary code with weight 2) for the jth place 
(bit). 

For coding the rows and columns of the matrix which 
represents the Boolean function to be minimized, we 
shall use a symmetric Gray code which, as will be shown, 
has certain superiorities over other codes. 

Unit ("1") values of a Boolean function will be 
represented by dots in the corresponding matrix elements, 
blank elements will correspond to zero function values, 
and blackened matrix elements will attest that the function 
is not defined on the sets of argument values correspond- 
ing to these elements. The set of all elements thus marked 
will be called the representation of the function. 


TABLE 1 








Example 1. Table 1 represents the Boolean function 
which takes the value *1" on the sets * 00000,00010, 
11011, 01110, 10100, 10110, 00111, is not defined on the 
sets 10010, 11010, 10001, 01101, 11111 and takes the value 
"0" on the remaining sets. One may consider that Table 
1 represents all 2° = 32 Boolean functions which are 
defined on all the sets but differ in their values on the 
sets 10010, 11010, 10001, 01101, and 11111. 

We note that when the number of variables does not 
exceed four, the method suggested here coincides with 
the method of M. Karnaugh[19] (cf. also [8, 20]), but 
loses its similarity as ne number of variables increases 
(in Karnaugh's method functions of six variables are 
represented in three-dimensional space while, for the 
representation of functions of five variables , Richards 
[8] suggests carrying out a diagonal decomposition of the 
elements in the tables representing functions of four 
variables, etc.) . 

Let y = fi(fy, fy, - +», f,,) be a Boolean function 
whose arguments are functions of the binary variables 
X4, Xg, + «+, Xp + In order to express y as a function of 
X4, Xgr ¢ + #/X,, it suffices to find, for each set of values of 
Xa Xgr Xge ee ey Xp the corresponding values of fy, f,,..-, 
fr, and then to determine the value of y corresponding 
to these latter values. 

This permits the easy implementation of operations 
on Boolean functions by means of their matrix forms, 
represented by tables on tracing paper. The sheets of 
tracing paper are laid on each other, the function values 
in the corresponding elements are compared and, in the 
corresponding element of the result matrix, the resulting 
value is noted. 
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The given method is not only applicable to manual 
computations; by using bit-wise logical operations, one 
can easily program it for computers. 


TABLE 2 





. 
o 
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Example 2, Table 2 gives the matrix representa- 
tions of the incompletely defined functions f,; and f, of six 
variables, as well as the results of certain operations on 
them. 

Note 1. The character ofthe representation matrices 
does not depend on the method of coding their elements 
if this method is the same for all the interacting matrices, 
so that the coding was not given here. _ 

Note 2. In obtaining the matrix for f, it was supposed 
that the inverse (negation) of an undefined value of a 
binary variable is also an undefined value. 

We begin the consideration of the method of 
minimizing Boolean functions with a comparatively 
simple case—obtaining the minimal disjunctive form 
(m.af.), while assuming that, in the general case, the 
function to be minimized is not defined on all the 
sets of argument values. 

In general, obtaining the m.d.f. makes it easier to 
find other minimal forms, although it may turn out that, 
in making the transition from the m.d.f. to another form, 
this latter will not be absolutely minimal. 

An algorithm which solves completely the problem 
posed was described by Yu. L. Zhuravlev [7], but it entails 
scanning variants, which leads to a verylengthy com- 
putational procedure, that hinders its practical employ- 
ment . 

Using the terminology of [7], we shall denote by an 
interval of the kth rank a sert of matrix elements cor- 
responding to the conjunction of k variables, some of 
which may enter with a negation sign. 

We now define the rules for isolating intervals. 


TABLE 3 











* By aset we mean a conjunction of argument values, 
for example the set 01000 corresponds to the values 
Xy = 0,xX, = 1,X, = x = Xs = 0. 
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1, Vertical lines in the table which separate 
elements corresponding to different values of xy we shall 
call axes of symmetry of the first rank; lines separating 
elements corresponding to different values of x, we shall 
call axes of symmetry of the 2nd rank, etc. We apply 
an analogous nomenclature to the table's horizontal 
lines (cf. Table 3). 

2. A set of elements which form a band of 
elements of width 2° , the central line of which is axis 
of symmetry of the rth rank,shall be called a zone of 
symmetry of this axis. (The hatched boxes on Table 3 
indicate zones for the first rank horizontal axis and 
second rank vertical axis.) 

3. Each element in the table for a Boolean function 
ofn variables is an nth-rank interval. 

A set of two kth- ‘tank intervals is an interval of the 
(k-1)th rank if the intervals are symmetric with respect 
to some axis of symmetry and are found completely 
within its zone of symmetry. 

We remark that, since the minimum width of a band 
of symmetry equals two, it follows that two series of 
elements forming a rectangle are always intervals; 
analogously, four series of elements forming a square are 
always intervals. 

The task of obtaining the minimal disjunctive form 
of a Boolean function reduces to finding such a decom- 
position of the function representative into intervals for 
which the sum of the ranks of all the intervals will be 
minimal. The decomposition is obtained by such a cover- 
ing by intervals of all the table elements corresponding 
to unit function values for which interval intersections 
are admitted, and covering by intervals of elements on 
which the function is undefined, but not admitting the 
covering of elements which correspond to zero function 
values. 

Example 3, To find the m.d.f. of the function 


f (a, b, c) = ab + be +- ae. 


TABLE 4 


A808 


£9 ab + be om f 


TABLE 5 


By successively marking with dots the elements o1 
the intervals which correspond to the conjunctions ab, be, 
and ac, we obtain the representation of the function 
(cf. Table 4, a). 

From the visual perception of the representation it 
becomes obvious that the advantageous decomposition 
of the representation is into the two intervals corresponding 
to the conjunctions ab and ac (Table 4, b). Thus easily 
do we obtain the m.d.f. of the original function f(a, b, c) = 
= ab+ ac. 

The special convenience of the visual method comes 
to the fore in the minimization of incompletely defined 
functions. 

Example 4. To find the m.d.f. of the function whose 
representation is shown in Table 5. 

The series of coverings of the representation by 
intervals is shown in Table 5, given only to elucidate the 
decomposition principle, since in practice one of the 
original matrices can be avoided by successively covering 
the representation by intervals and by blackening the 
elements which fall in intervals. With a little practice, 
one can carry out the decomposition of the representation 
mentally. 

The result of the decomposition gives the m.d.f. of 
the function to be minimized. 


f = 2y0rg + Lerere + Teste + 71% + Lorex. 


As this example shows, the visual estimate of the 
representation permits one to find the optimal variant 
of the incompletely defined function more rapidly 
than may be done with algebraic methods [6]. 

In certain cases, the algebraic form of the original 
Boolean function which contains a minimal number of 
variable symbols may be obtained in the form of a 
product of two disjunctive forms. 
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TABLE 6 
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Example 5. (Taken from R. K. Richards’ book [8], 
page 70.) To find the minimal algebraic form of the 
Boolean function f = ab + acd + be + bde+ade+ cde. The 
function's representation {s given as the intersection of 
two other representations, f, and f,, easily decomposed 
into intervals (Table 6). The final result, 

f = (a+ be + de) (b+ a+ cd), contains nine variable 
symbols.t 

The material herein presented allows the following 
conclusions to be drawn : 

1. The matrix representation of Boolean functions 
in which the matrix rows and columns are coded by the 
variables’ values in a Gray code has the convenience of 
immediate visual perception of certain properties of the 
function represented, the correct estimate of which is 
necessary in the finding of the optimal way to minimize. 


Visual analysis of the representation allows one to estimate 


the possibilities of minimization presented by the re- 


presentation by degrees of regularity of the representation's 


elements, the criterion for which is the symmetry of the 
elements’ mutual positioning. 

2. Being intuitive to a certain extent, the visual 
matrix method of minimizing Boolean functions was 
designed for using certain properties, difficult to render 
in algorithmic form, of visual perception, thus allowing, 
without recourse to computing technology, the imple- 
mentation of complicated transformations of the original 
function for its minimization. With a little practice, 
the minimization of Boolean functions of five or six 
variables may easily be done mentally (if the re- 
presentation of the function is in front of one’s eyes), 
while making notes on the representation elements 
during the minimization process allows one to increase 
the number of variables up to ten or twelve. 
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T This is one symbol fewer than in Richards’ result, 
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Gestalt perception of the representation which allows 
the connection between the representation elements to 
be taken into account in the search for the optimal 
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OF ELECTROMAGNETIC SYSTEMS 
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Questions of proportionality of electromagnetic devices are propounded. Geometric regularities, related to the 
conditions of rational use of electromagnetic energy, are used. As a result, there are found many optimal configura- 


tions from which, in correspondence with the technical and economic 


in each specific instance. 


In devices which use electromagnetic energy, the 
basic design element is the electromagnetic system, 
which is a ferromagnetic magnetic circuit with one or 
several windings. 

The electromagnetic system of an electromagnetic 
device is a function of the device's purpose. Thus, for 
example, there are electromagnetic systems for transfor 
mers, chokes, magnetic amplifiers, relays, and other 
electromagnetic mechanisms, among which can be 
included electrical machines. 

Despite the variety of the physical forms taken by 
electromagnetic systems, they have something in common, 
which speaks for the efficacy of a generalized consider- 
ation of the fundamental regularities which are valid 
for all electromagnetic devices. 

For this generalized consideration of electromag- 
netic systems, we use the fact that any electromagnetic 
system is some geometric figure whose absolute and 
relative dimensions are determined by the physical laws 
concerning the use of electromagnetic energy. 

We are justified in calling the set of questions as to 
the geometry and physics of electromagnetic systems the 
"geometry of electromagnetic systems .” 


Basic Postulates of the Geometry of 
Electromagnetic Systems 

Questions as to the geometry of electromagnetic 
systems in their generalized form have not yet been 
treated in the literature, so that it is necessary to first 
consider certain postulates and, for convenience of 
exposition, to introduce some terms from electromagnetic 
system theory. 

First of all, it is necessary to introduce the concept 
of the absolute magnitude of an electromagnetic system. 

The necessity for this arises in all cases when it fs 
required to compare different electromagnetic systems 
among themselves. Such comparisons can be made with 
sufficient simplicity if the electromagnetic systems to 
be compared are completely similar. In all other cases, 








requirements, the actual configuration is chosen 


this comparison gives rise to well-known difficulties. 

We shall henceforth take, as the absolute magnitude 
of an electromagnetic system, the sum of the volumes of 
its active materials, i.e., the sum of the volumes of its 
windings and magnetic circuit. With this, any change in 
the dimensions of the electromagnetic system without 
any concomitant change in the sum of the volumes of its 
active materials will be considered as a change of its 
relative dimensions with retention of its absolute mag- 
nitude. Conversely, we shall consider that a change in 
an electromagnetic system's absolute magnitude cor- 
responds to a change wherein all linear dimensions 
change while the conditions of geometric similarity are 
retained. Finally, if a change in the dimensions of an 
electromagnetic system is not accompanied by a 
retention of the conditions of geometric similarity and, 
simultaneously, the sum of the volumes does not remain 
a constant, we shall then consider that, in this case, both 
the absolute magnitude and the relative dimensions of 
the electromagnetic system have changed. 

As the relative linear dimensions of an electromag- 
netic system, we take its linear dimensions divided by a 
fixed one of them. 


Corresponding to the concepts of absolute magnitude 
and relative dimensions of an electromagnetic system are 
the concepts of the absolute and relative parameters of 
this system. As the relative (or reduced) parameters we 
take the absolute parameters divided by the absolute 
magnitude of the electromagnetic system. With this, if 
the parameters which interest us are proportional to the 
first powers of linear dimensions, then their reduction will 
be implemented by a division by the cube root of the 
sum of the first powers of the active materials’ volumes. 
For a quadratic dependency, the aforementioned sum 
will be raised to the second power , etc. 

By using the fact that the parameters of electro- 
magnetic devices can be reduced to an arbitrary number 
of turns of their windings, we shall consider their windings 
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to be reduced to the one-turn form. In this case, instead 
of thecross section and length of a circuit's winding, we 
shall consider its mean length and total cross section,with 
account taken of the corresponding space factor. 

Naturally, all parameters of electromagnetic systems 
are the corresponding functions of the geometric 
dimensions of these systems and of the properties of the 
active materials. For convenience of exposition, we shall 
henceforth, instead of such expressions as "a geometric 
function of the magnetic circuit's magnetic conductance,” 
use “the magnetic circuit's geometric conductance,” "the 
geometric impedance of the winding,” and also "geo- 
metric coefficient. . .” or “geometric constant. . .,” etc. 

The nonlinearity of ferromagnetic materials ordinar- 
ily introduces significant complication into the theory of 
elemectromagnetic devices, but it is not mandatory to 
take these nonlinearities into account in the theory of the 
geometry of electromagnetic systems. The explanation 
of this is that the geometry of electromagnetic systems 
is concerned with the regularities imposed on the relative 
(or reduced) parameters, in the aforementioned sense, of 
electromagnetic systems, and these reduced parameters 
do not depend on the properties of any material. As will 
become clear in the sequel, the materials’ properties, 
while having no infleunce on the choice of relative 
dimensions of electromagnetic systems, influence com- 
pletely its absolute dimensions. 

We should make particular mention of the questions 
in our generalized consideration of electromagnetic 
systems which deal with the classification of systems into 
those with closed or with open magnetic circuits. Thus, 
for example, a transformer has a closed, and an electro- 
magnetic relay an open, magnetic circuit. 

The author initially will consider all electromag- 
netic systems in the closed magnetic circuit mode (the 
mode of a short-circuited magnetic circuit) The tran- 
sition from closed to open magnetic circuits does not 
render it necessary to change the optimal geometric 
relationships obtained when considering closed mag- 
netic circuit modes. ‘The proof of this assumption in 
particular cases of electromagnetic system employment 
engenders no particular difficulty, and is in good agree- 
ment with experience . 

The consideration of proportionality of electromag- 
netic systems with closed magnetic circuits eliminates 
the necessity of taking magnetic flux leakage into 
account, which also simplifies the problem. 

The reduction of windings to one turn and the ignoring 
of flux leakage lead to an equivalence between current 
and ampere-turns, and also between magnetic flux and 
flux linkages. 

Energetic Bases of Electromagnetic 
Systems 

The use of electromagnetic energy in any electro- 
magnetic device is related to the property of the electro- 
magnetic system of accumulating energy obtained from 
a source of electric current. This energy is then partially 
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or completely utilized, and is partially converted to a 
source of supply or is lost in the form of various energy 
dissipations. 

As the characteristic of this property of electro- 
magnetic systems, we can use the quantity of accumulated 
energy divided by the power used in maintaining the 
static state of the energy stored in the system. To find 
this parameter we shall use regularities which exist for 
de circuits, but it should be kept in mind with this that 
the parameter in which we are interested does not depend 
on the nature of the current. 

As is well-known , when the windings of an electro- 
magnetic system are connected to a de source, a transient 
response ensues for which the following equation fs valid; 


uidt = i*rdt + idd. (1) 


Equation (1) shows that all the energy from the source 
consists of two parts, one of which is consumed in thermal 
losses in the windings, while the other serves to change 
the system's magnetic state. At the termination of the 
transient response, the current and the magnetic flux 
assume constant values. The energy which continues to 
be supplied from the source is now consumed solely by 
the thermal losses in the windings. These losses determine 
the power which goes to maintain the static state of 
electromagnetic energy accumulated in the electromag- 
netic system during the transient response. This power, 
defined in terms of the steady-state value of current i 
and the electric conductance gy Of the windings, equals 


je 
P=. (2) 
The quantity of electromagnetic energy stored in 
the electromagnetic system during the transient response 
is defined by the expression 
T 
W.= \idg, 


0 


(3) 


To find the differential magnetic flux which appears 
in Expression (3), we use the well-known dependency 
which exists between the magnetizing force of magnetic 
flux and magnetic conductance. With this, and also using 
(2), we obtain 
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where Gy and G, are the geometric conductances of , 
respectively, the windings and the closed magnetic cir- 
cuit, and k is a proportionality factor which takes into 
account the adopted units of measurement, the space 
factors of the winding cross section and the magnetic 
circuit, the specific conductance and permeability of 
the materials of the windings and the magnetic circuit, 
and other quantities relating the windings’ electric 
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conductance and the magnetic circuit's magnetic 
conductance with the geometric dimensions of the 
electromagnetic system. 

By dropping the factor k from (5) , we obtain a func- 
tion of the electromagnetic system's geometric 
dimensions which characterizes this system's property of 
accumulating electromagnetic energy. With this, we 


obtain 
SS 
=G.6. = MC 
7 m“Cc oe 





(6) 


Here, Syy and S,, are the cross sections of the winding 
and of the magnetic circuit, and L), and L,, are their mean 
lengths. 

As is clear from (6), to improve the properties of 
electromagnetic systems, it is necessary, in choosing 
their geometric dimensions, to attempt to increase the 
product of the cross sections of the winding and magnetic 
circuit, and to decrease the product of their lengths. This 
is the basic law of the theory of electromagnetic systems. 
However, (6) cannot be used for comparing different 
electromagnetic systems among themselves, since the 
quantity defined by this formula depends not only on the 





relative, but also on the absolute,dimensions of the electro- 
magnetic systems. To eliminate the dependence on 
absolute quantities, we use the relative value of the 
parameter of interest to us. We obtain with this 


p_ 1 Sudo 
vouente’ 7) 


where Q is the sum of the volumes of the active materials, 
i.e., the absolute magnitude of the electromagnetic system. 

In the theory of the geometry of electromagnetic 
systems, (7) is one of the fundamentals. In the sequel, 
this formula will be used for the analysis of the propor- 
tionalities of the geometric forms of electromagnetic 
systems. The parameter defined by (7) is called the 
"geometric constant” of the electromagnetic system. 

The introduction of the concept of the geometric 
constant into the analysis of proportionality of electro- 
magnetic systems allows geometric function (6) of any 
electromagnetic device to be presented in the form of two 
factors, one of which is determined by relative, the 
other by absolute, dimensions of the electromagnetic 
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Fig. 1. Certain possible geometric forms of electromagnetic systems. 












system, i.e., in the case of the geometric function given For this purpose, as an example, we consider the 
by (6), we shall have geometric form of an electromagnetic system as shown 
_ oshms in Fig. 2. 
‘< Ve. (8) As is clear, the given form has four dimensions, the 

absolute values of which are denoted by a, b, d, ande, 

To obtain relative dimensions, we relate all the 
remaining dimensions to one selected absolute dimension, 
For example, as such a fundamental dimension, we take 
the width a of the magnetic circuit's aperture . In this 
case, the relative dimensions of height, width of core 
section,and magnetic circuit thickness will equal, res- 
pectively, 






As is obvious, for a given magnitude of the electro- 
magnetic system's geometric function y , to decrease the 
expenditure of active materials, i.e., to decrease the 
quantity Q, it is necessary to try to maximize the geo- 
metric constant € . 

However, it should be taken into account with this 
that the prices of the materials of the windings and of 
the magnetic circuit are different, so that it is necessary 
to attempt to realize, not an absolute decrease in the —=C, aan, a k. 
sum of the volumes of the windings and the magnetic 
circuit, but a reduction of the total cost of the electro- 
magnetic system's active materials. 

Linkages of Geometric Dimensions of 
Electromagnetic Systems 

Electromagnetic systems, in their geometric aspect, Ly = 2a(m--k+- 2), L, = 2a(m +c +- 1). (10) 
are complicated figures which consist of linkages 
of simpler figures of windings and magnetic circuits. 

Linkage between windings and magnetic circuits may 
be simple or complex. In a simple linkage, the entire 
aperture of the magnetic circuit is filled by the winding's 
cross section, and the winding’s aperture is completely 
filled by the magnetic circuit. In a complex linkage, 

between winding and magnetic circuit there remain Q = 2a®[c(m + k + 2)-+- mk (m+ c+ 1)). ay 
certain spaces, filled by the surfaces of winding insulation 
or by other windings, or these spaces, for some reason or 
another, may remain unfilled. 

As examples, Fig. 1 shows a number of possible 










With the notation chosen the cross sections and 
lengths of the windings and magnetic circuit will equal 









Su = ac, Se = a*mk, (9) 















Knowing the cross sections and the lengths, we 
determine the volumes of the windings and magnetic 
circuits and then, by adding these volumes, we find the 
quantity which we have called the absolute magnitude 
of the electromagnetic system. For this, we get 









To determine the value of the parameter which we 
have taken as the fundamental one, we solve the last 
equation for a ; 




















geometric forms of simply-linked electromagnetic 3 0 (12) 
systems. Certainly, any simply~-linked geometric form, oinns V ; [c(m + k-+2)+mk(m+c + 1)|° 

with any gemmiaeite — ot en apertare . becomes We now have all the data necessary for determining, 
a figue with « complex age by (7), the geometric constant of the electromagnetic 






system under consideration. This constant equals 
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(13) 









3 
(m + k + 2) (m + ¢ + 1) ¥16* [ce (m + & + 2) + mk (m+ ¢ + 1)/? 





In our task here, we do not enter into an analysis of The problem of finding the optimal proportions of 
all the possible geometric forms of electromagnetic electromagnetic systems does not have a unique solution 
systems. We limit ourselves just to a consideration of a In each particular case it is necessary to take into accoull 
method of analysis which is applicable, with the proper the concrete conditions of use of electromagnetic energy 


specifications being made, to the analysis of any form and and the economic requirements. This bears first of all 
of any type of linkage. on the ratio of the volumes of the electromagnetic 


system's active materials. By bearing this fact in mind, 










































































































itil Ii ti we should seek, not the general case of the maximum of 

sy fa Hi IH the geometric constant of (13), but its particular cases, 

4 + f+ te Hit oa each of which corresponds to some previously given ratio 

| | | of the aforementioned volumes. With the notation we ha 

mn | % adopted, the ratio of the volume of the magnetic circult 

OS5d ya) ie an to that of the — equals 

mi e+ 4 

Fig. 2. Geometric form of a three-core o- = fe = 7" ‘ (14) 
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TABLE 1 
m 0.0 0.25 
I 3.0 3.0 

















By solving this last equation for one of the relative 
dimensions, and then substituting the value obtained in 
(13), we would have the capability, by computing the 
partial derivatives with respect to the relative dimensions, 
to determine the optimal values of these dimensions, 








which would correspond to the maximum geometric 
constant for the previously given constant ratio of the 
volumes of the active materials. However, this path to 
a solution leads to significant complexity, so that we 
shall solve this problem by finding the maximum of the 
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Fig. 3. Series of optimal geometric forms of three-core electromagnetic 
systems for various ratios of winding and magnetic circuit volumes. 












geometric constant of interest to us for a constant value 
of one of the relative dimensions, which is equivalent to 
the condition that the ratio of the volumes be constant, 
although the magnitude of this ratio itself is not given 
explicitly inthis case. As the constant quantity, we take 
the relative dimension m_ , since the ratio of the winding 
and magnetic circuit volumes depends most heavily on 
this dimension. 

By analyzing (13) for an extremum by the method 
of partial derivatives, we obtain the following optimal 
values of the relative dimensions k and c for a constant 
value of dimension m ; 


ne De |. as 
Pic a [1 tW 14am a |! 
24mk 


m + 1 16 
“no [1 + V1 + m +-k -- 2 + mk | si: 


Joint solution of (15) and (16) leads to cumbersome 
computation, so that we solve this problem by the method 
of successive approximations. The results are given in 
Table 1, where we also give the volume of the magnetic 
circuit as a percentage of the total volume of active 
material plus the corresponding maximum of the geome- 
tric constant for the electromagnetic system of the form 
under consideration. 

Thus, the imposition of the condition that the 
geometric constant be a maximum leadsto a completely 
definite regularity of variation of the optimal relative 
dimensions of the electromagnetic system as functions 
of the ratio of the volumes of its active materials. As a 
result, we obtain a continuous series of optimal variants, 
from which one should choose that variant which most 
completely answers to the technological and economic 
requirements of each given particular case of electro- 
magnetic system usage. 

















As an illustration of the regularities obtained, Fig. 3 
shows a series of forms of the type of electromagnetic 
system considered. All the sketches on the figure are 
drawn to scale and are in accordance with the condition 
that the geometric constant be a maximum. The order in 
which the sketches are placed on the figure takes into 
account the ratio of the volumes of active materials for 
a constant value of their sum. The limits of the series 
show the sketches with limiting ratios of their dimensions. 
On the top are shown the sketches corresponding to zero 
magnetic circuit volume and, on the bottom, those 
corresponding to zero winding volume. 


Knowing the optimal ratios of the dimensions, there 
is no difficulty in determining all the parameters of 
electromagnetic systems for any ratio of the volumes 
of the active materials. Thus, for example, the 


curves of Fig. 4 show, as functions of the ratio of the 
volumes of active materials, the relative geometric 
conductances of the windings and magnetic circuit as 
well as their product, i.e., the geometric constant of the 
electromagnetic system under consideration. 
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Fig. 4. Geometric conductance of wind- 
ings and magnetic circuit plus the geo- 
metric constant of a three-core electro- 
magnetic system as functions of the ratio 
of winding and magnetic circuit volumes. 
1 is the geometric constant, 2 is the mag- 
netic circuit conductance, and 3 js the 
winding conductance. 

It is interesting to note that the geometric ratios 
which correspond to the regions of the maximum of the 
maxima of the geometric constant (Figs. 3 and 4) are 
the ratios which have been most widely used in the 
practical application of the electromagnetic systems of 
the geometric form considered. This substantiates the 
correctness of the basic assumptions we have made for the 
theory of the geometry of electromagnetic systems. 




































































TABLE 2 

Sketch Form Sketch Form 

on figure} coefficient] on figure | coefficient 

x 16 x 108 

lea 3.40 le 6.70 

1. b 2.46 1. f 5.00 

le 3.62 1. g 8.00 

1d 4.85 lh 6.40 

















The method presented here for the analysis of 
optimal geometric proportions is applicable to any 
geometric form of electromagnetic systems. Thus, for 
example, the author carried out the corresponding analys 
of all the sketches given in Fig. 1. Omitting similarities 
we give only the value of the maximum of the maxima 
of the geometric constant for each of these sketches. The 
results appear in Table 2. The maximum of the maxima 
of the geometric constant may serve as the characteristic 
of the form of an electromagnetic system, for which reas® 
we call this parameter the “form coefficient” of an 
electromagnetic system. 














Economics of Electromagnetic Systems 





In choosing some variant from the series of optimal 
design variants of an electromagnetic system (cf. , for 
example, Fig. 3) for some particular case, one must take 
into account, not only technological, but also economic. 
requirements. 

Thus, economic considerations,while having no effect 
on the general regularities of the geometry of electro- 
magnetic systems, have an effect, in each particular 
case, on the system's geometric proportions. It should 
be kept in mind with this, that the basis of the economic 
influence is the difference in prices of the winding and 
of the magnetic circuit materials of the electromagnetic 
system. With identical prices per unit volume of these 
materials, the economic requirements would cease to be 
operative. 

The cost of an electromagnetic system, defined by 
the costs of its active materials, equals 

A = Qutu + Qeles (17) 
where Qyy and Q, are the winding and magnetic circuit 
volumes, and ty, and t, are the unit volume costs of, 
respectively, the windings and the magnetic circuit. 

The ratio of unit volume cost of winding to that of 
magnetic circuit equals 
: (18) 


. 
fo 


By introducing the coefficient 
seiiatiiemmal 
~ QutQ, Q’ 


?°~= 


(19) 


and by taking (18) into account, we transform Expression 
(17) to the form 


A = Qt, [n (to — 1) + 4]. 


By dividing (20) by Qt,, we obtain the relative 
cost of the electromagnetic system, Le., 


(20) 


Ay = Gy =R(to—1) +1. 


As fs clear, the relative cost when the prices of the 
active materials are equal (i.e., for t = 1) is unity, and 
does not depend on the ratio of the volumes of these 
materials. In this case, the optimum variant of the 
electromagnetic system to use is that one which corres- 


A 
f, (21) 


ponds to the maximum geometric constant. When these 
prices are not equal, a change in the ratio of the volumes 
of active materials changes both the geometric constant 
of the electromagnetic system and its cost. In this case, 
to choose the optimum variant of electromagnetic 
system to be used, one should aim at the minimum ratio 
of the relative cost to the geometric constant, {.e., 


PEER woe (22) 

Thus, for example, from the data of Table 1, if no 
account is taken of differences in active material costs, 
i.e., for t = 1, the optimal variant of the electromagnetic 
system to use is the one for which the winding volume is 
42 % of the total volume of active materials. If the price 
of the winding material is, for example, 20 times as high 
as the magnetic circuit material, then from the data of 
Table 1 and from (22) , we find that the optimal variant 
for this case is the one where the winding volume is 24% 
of the total volume of active material in the electro- 
magnetic system. 

What has already been said sufficiently adumbrates 
the principal interrelationships which exist between 
geometric proportions and the economics of electro- 
magnetic systems. 


Conclusion 

In this paper we touched on only the most general 
postulates of the theory of the geometry of electromag- 
netic systems. In the future, by basing ourselves on 
these regularities, we can find optimal relationships for 
standard electromagnetic devices; transformers, electro 
magnetic relays, magnetic amplifiers, induction relays, 
chokes, electrical machines, etc. 

As the final goal, the consideration of any configura- 
tion of the electromagnetic system of any electromagnetic 
device should allow one to obtain a continuous series of 
optimal variants of the configuration considered, the 
sequence of these variants being determined by the ratio 
of the volumes of windings and magnetic circuit. The 
study of such a series provides the capability of establishing 
the interdependencies of the relative dimensions and 
the technical-economic parameters which are of decisive 
value in choosing the physical design of the electromag- 
netic device from this series in each specific case. 
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The paper deals with relay winding overheating and heat transmission coefficient as functions of ambient air 
temperature. Empirical formulas are derived for determining the magnitude of relay winding overheating and relay 
heat transmission coefficient. 










1. Winding Overheating the absolute temperatures of the winding and of the am- 




















The excess of relay winding temperature over the bient air. Therefore, with identical magnitudes of power 
ambient medium's temperature is ordinarily determined dissipation, overheating of a relay winding is decreased 
under normal conditions for an air temperature of +20°C. as the ambient air temperature increases. 

However, relays are frequently used under conditions The working life of relay insulation is determined 
where there are large oscillations of the ambient temper- _ by the greatest magnitude of the working temperature. 
ature — from -60 to +60°C-— and special types of relays Consequently, it is necessary to know, as accurately as 
are used at temperatures up to +125°C and higher. possible, the relay winding overheating for the highest 

As is well known, with a raising of the winding ambient air temperature. Measurement of relay 
temperature the heet transmission coefficient (specific winding overheating with increased temperature is 
heat transmission) increases, since heat transmission by ordinarily carried out experimentally in a thermostat, 
convection depends on a fractional power of the overheat- for which a great deal of time is required. Moreover, 
ing temperature, while heat transmission by radiation is these measurements are usually distinguished by their low 
proportional to the difference of the fourth powers of accuracy, due to the significant relative errors in 
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Fig. 1. 1—@ = 21°C, 2 —@ = 44°C, 3—6, = 65°C, 4— 0 =84°C, 
. ™ 5 — G = 105°C, 6 — 0, = 141°C. 

















determining the winding overheating temperature for an 
increased temperature of the ambient air. 

Thus , for example, with ambient temperature of 
100°C and winding overheating temperature of 30°C, an 
error of 2%(2C) , in measuring the temperatures of 
winding and air gives a relative error of about 13% in 
determining the temperature of the relay winding over- 
heating. 

Relay winding overheating temperature can be found 
more accurately by analytical means. 

Figure 1 gives curves showing excess of relay winding 
temperature over ambient temperature as a function of 
the power applied for small-size-type RMU relay, taken 
off experimentally for various ambient temperatures up 
to+149P°C. 


The measurements were carried out in a 500-watt 
dessicator with automatic temperature regulation within 
the limits of +30 to +200°C with an error no greater than 
+ 3%. The working chamber of the dessicator had a 
diameter of 300 mm and a length of 290 mm. 

The dimensions of the relay types tested were: for 
type RMU, 24.7 x 38.5 x 41 mm, for type RES10, 11.1 x 
x 16.5 x 26 mm and, for type RKN, 25.6 x 56.5 x 97 mm. 


The relay weights were; for type RES10, 6 g, for 
type RMU, 70 g and, for type RKN, 290 g. 

The relay specimens had wire windings with brand 
PETKM heat-resistant insulation. The winding tem- 
peratures were measured by a copper-constantan 
thermocouple. To eliminate heat leakage at the relay 
supports, the relays were suspended by a fine wire in the 
center of the dessicator. 

The ventilation apertures of the dessicator were 
tightly covered by asbestos plugs, and were used for the 
outlet of the conductors and the thermocouple. 

The air temperature in the dessicator was also 
measured by a thermocouple,which was placed at the 
level of the specimen relay at a distance of about 30 mm 
from it. As protection against direct radiation, the 
thermocouple'’s junction on the relay'sside was shielded by 
a small fragment of white paper. For each value of 
supplied power, the relay was immersed for not less than 
an hour in the given ambient temperature. The actual 
temperature oscillations in the dessicator did not exceed 
2%, It is clear from Fig, 1 that, in the limits of over- 
heating temperature from 20 to 70°C, the curves differ very 
little from the straight lines drawn with dashes (the error 
does not exceed +1I°C). In Fig, 1, the dashed straight 
lines are shown only for the extreme curves, so as not to 
complicate the figure. 


In the aforementioned limits, the relay winding over- 
heating curves can be approximated by the formula 


}=—a+ bP, (1) 
where a fs the amount of overheating intercepted by 


the straight line on the axis of ordinates, and b is the 
tangent of the slope of this line. 


TABLE 1 





, Ambient . 
temperature a,° 
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b, °C/W 
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44 
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84 
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The values of the quantities a and b for relay type 
RMU for different values of ambient air temperature, 
determined graphically from Fig, 1, are given in Table 1. 
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From the data of this table, the dependencies of a 
and b on the ambient temperature were constructed, as 
shown in Fig. 2. These dependencies are straight lines, 
which can be expressed by the formulas 


a = Ay — CH, (2) 
b = bo = d§,, (3) 


where ap and by are the values of a and b for an ambient 
temperature of (°C , and c andd are the tangents of the 
slopes of these lines. 

By substituting in (1) the values of a and b from 
(2) and (3), we obtain a general expression for determining 
the relay winding overheating temperature as a function 
of the supplied power and the ambient temperature ; 


9 = dy — Cy + (bo — dq) P. (4) 


This formula gives sufficient accuracy only in the 
limits of overheating values from 20 to 70°C, which 
are the ones most frequently occurring . 

The largest errors are obtained at the initial portion 
of the curve, for overheating temperatures less than 20°C. 
The accurate determination of such overheating temper- 
atures ordinarily has no practical value. 

From Fig. 2 we find the values of quantities 4, c , 
byand d for relay type RMU: a = 50.7°C, c= 0.0175, 
by = 22.54°C / watt and d = 0.038, 


Consequently , for type RMU relay, winding over- Table 3 gives the values of the temperature over- 
heating ts expressed by the following formula: heating coefficient 6 for various values of power in 
$ = 5,07 —0,0175 6 + the winding, and the corresponding relay winding over- 
+ (22,54 — 0,038 6) P. heating temperatures of relay types RMU, RES10 and RRN 
By an analogous procedure, we found the values of the for these values of power and for ambient temperature 
quantities a, c , by,andd for miniaturized type RES10 +20°C, computed by means of (4) and (6). 
relay (without cover) and normal type RKN relay. Figure 3 shows the curves, constructed from the data 
The values: of these quantities are given in Table 2. of Table 3, of the temperature overheating coefficient 


(4a) 


TABLE 2 
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Relay type , be» ae 





RMU 0.038 


R 
REsio (in cover) 























RES10 (without covel , 





From the data given in Table 2, it is clear that, with 

increasing relay dimensions, the values of quantities b, 

and c increase, and the values of by and d decrease. 
Determination of the quantities a), c , by,and d 

for each type of relay is quite complicated and time- 

consuming, so that we find the relationship between the 

magnitudes of relay winding overheating for different 

ambient air temperatures and one and the same magnitude 

of power. 60 Ae 
Relay winding overheating for one and the same 


power P and different ambient temperatures 9 and Fig. 3. 1) Relay type RMU; 2) relay type 
O92 will equal RES10 without cover; 3) relay type RKN , 


9, = dy — c85; + (by — 951) Pand 4) relay type RES10 in cover. 
8 as a function of the windingoverheating temperature of 
9 = dy — Boa + (bo — don) P, these types of relays for 99 : +20°C. = aa 
purposes , the dashed line of Fig. 3 shows the approxi- 
mate curve for type RES10 relay (with cover). It follows 
from Fig. 3 that the curves for the coefficient 8 for rela 
$2 = 9, —(c¢ + Pd) (O92 — 9%) = 5 types RMU and RES10 (without cover) virtually coincide, 
91 [1 — B (Oo. — %;)], @) and the value of coefficient 8 for type RKN relay fs 
. about 10 % larger. 
Type rés10 relay, protected by a thin aluminum 
cover, has a coefficient 8 of about half the size, due to 


the poor thermal contact between the relay's magnetic 


from which we find the following expression for relay 
winding overheating for ambient temperature 99, : 


TABLE 3 





Relay type 
RES10 (without cover) 
P, w B $, °C 











0.00214 19.77 
0.002005 | 33.58 
0.001955 | 47.4 
0.00193 | 61,22 
0.001905 | 75.03 
0.001895 | 88.86 
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circuit and the cover (the relay’s pedestal is made of 
plastic). With air-tight relays having thick steel covers 
and fastened to a metallic base (pedestal), the magnitude 


of coefficient 8 is about the same as for relays unprotected 


by covers, since the winding of an air-tight relay 
generally has good thermal contact with the magnetic 
circuit and with the cover. As the relay winding over- 
heating temperature varies within the limits of 20 to 70°C 
{the limits of applicability of (4)], the size of coefficient 
6 varies comparatively little (by not more than + 10 %). 
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The magnitude of coefficient 8 depends on the initial 


temperature of the ambient air. Figure 4 gives the curves 
for the coefficient 8 as a function of winding overheating 
for various ambient air temperatures for type RMU relay. 
By means of these curves, the curves shown on Fig. 5 

of coefficient 6 as a function of ambient air temperature 
for various heating temperatures were constructed. The 
curves of Fig. 5, for 8 within the limits of 0 to 100°C, are 
practically straight lines, so that coefficient 6 as a 
function of initial ambient air temperature with average 
overheating temperature ( $ = 40°C) can, for small-size 
relays (without covers), be expressed by the formula 


B = (1.91 + 0,00368,,) 107%. (6) 


Consequently, with a constant magnitude of power, 
if one knows the winding overheating temperature of 
small-size relays (without covers) for the ambient tem- 
perature 99; , one can , with an accuracy sufficient for 
practical purposes, determine the overheating of these 
relays’ windings for another ambient temperature 992 by 
means of the formula 


d. =, [i —(1.94 ot 0.0036 601) 10-3 (92 val 901)).¢7) 


If the initial ambient temperature is 99; = 20°C, 
then 

Thus, with an increase of ambient air temperature 
by 10°C (within the limits of 0 to 100°C), the winding 
overheating of the relay (without cover),for one and the 
same value of power, is decreased by about 2%. 





Measurements carried out on different types of relays 
showed that this formula {s also valid for higher tem- 
peratures of the ambient air (up to 140°C). 


2. Heat Transmission Coefficient 





The magnitude of the average (nominal) heat trans- 
mission coefficient of windings of electromagnetic relays 
and dc apparatus can be determined from the experimen- 
tal data by means of the formula [1] 

[ee 
ead aS,’ (8) 
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where P is the power consumed by the winding, # is the 
overheating of this winding, and S;, is the computed 
cooling surface of the relay winding. 

Figure 6 gives the experimental curves which show 
the mean heat transmission coefficient of the winding 
of a type RMU relay as a function of its overheating for 
various ambient temperatures. W ithin the limits of an over- 
heating temperature of 20 to 70°C, these curves differ 
very little from straight lines. The large scattering of 
points for overheating values less than 20°C is explained 
by the inaccuracy of measurements of small temperature 
drops. 
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Fig. 7. a shows the mean heat transmission coefficient 
of relays as functions of winding temperature; 1 is for 
relay type RMU and 2 is for relay type RKN; b shows the 
same thing for relay type RES10 (without cover). 
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To elucidate the relationship between relay winding 
heat transmission and temperature, Fig. 7 shows the cons- 
tructed relationship of the average heat transmission coef- 
ficient of relay windings and the temperature of these 
windings. Within the limits of 30 to 210°C, this relationship 
is virtually a linear one which may be expressed by the 
formula 


= Jo + €8 = go + € (0, + 9), (9) 


where % is the magnitude of the heat transmission coef- 
ficient, intercepted on the axis ordinates by the extension 
of the straight line, and e is the tangent of the slope 
of this line. 

From Fig. 7 we find the values of q@ ande for type 
RMU relay : d = 1.422 x 10 w/cm?/*c and 
e = 0.00507 x 10° w/cm*/*c’. 

Thus, within the limits of 30 to 210°C, the mean 
heat transmission coefficient for type RMU relay is ex- 
pressed by the approximate formula 


gi = (1.422 + 0.00507 6) 107%, 


The greatest deviation of the individual values of 
heat transmission coefficient of type RMU relay from the 
straight line does not exceed 2.2% . The values of dg 
and e for relay types RMU, RES10, and RKN are given in 
Table 4. 


(9a) 











TABLE 4 

Relay type dw/eonF-*c | *& w/em? °C 
REs10(withoutcover)| 3.24 -40-° | 0.01965-40-8 
RMU 1,422-10-% 0.00507 - 10-* 
RMN 1:198-10- | 0.00567- 10-9 
RES10 (in cover) 2.77 -10-3 0.0045 -10-% 











The magnitude of the mean heat tramission coef- 
ficient for different temperatures 6, and 9, of the relay 


winding will be g’ = q, + e, andq, =q, + e:, 


from which we find the following expression for the 
mean heat transmission coefficient of the relay winding 
for temperature 05: 
q, = q, + e (6, — 4,) = q, [1'+ 7 (8. — 9,)], (10) 
where 
e 
{~=-7. 
% 
Figure 8 shows the curves of the transmission coef- 
ficient of the mean heat transmission coefficient y as 
a function of winding temperature for relay types RMU, 
RES10 , and RKN, computed by means of Formulas (9a) 
and (11). 
It follows from Fig. 8 that the curves of y = y(@) 
for small-size relay types RMU and RES10 (without cover) 
differ very little from each other, and that the value of 
coefficient y for type RKN relay is about 25% larger. 
Type RES10 relay, protected by an aluminum cover, has 
a value of coefficient y which is about half the size. 


(11) 
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Fig. 8. 1) Relay type RMU; 2) 
relay type RES10 without cover ; 
3) relay type RKN; 4) relay type 
RES10 in cover. 


In the limits of 20 to 160°C for ©, these curves differ 
very little from straight lines. 

In the aforementioned limits, the coefficient as a 
function of temperature for relay type RMU may be 
expressed by the formula 


+1 = (3.45 — 0.0079 8,) 107°. 


Thus, if we know the magnitude of the heat trans- 
mission coefficient of small-size type RMU relays (without 
cover) for winding temperature 9; , we can, with an 
accuracy sufficient for practical purposes, determine the 
heat transmission coefficients of these relays for other 
winding temperatures 6, (in the limits of 20 to 160°C) 
from the following approximate formula ; 


If the initial winding temperature is 6, = 50°C, then 
92 = qi [1 + 0.00306 (8, — 50)}. (13a) 

Formula (13a) shows that with a 10°C increase in 
winding temperature of small-size type RMU relays 
(without cover), the heat transmission coefficient is 
increased by about 3% . 

We find analogously that the increase in heat trans- 
mission coefficient with a 10°C increase in temperature 
is 2.87 % for relay type RES10 (without cover) and 1.5% 
for relay type RES10 in its cover. 

Thus, our work has established that overheating tem- 
perature of windings of small-size electromagnetic 
relays (weights up to 300 g) without covers, and also air- 
tight relays with good heat elimination from the windings, 
is decreased by about 0.2 % with 6»; = 20°C for a °C 
increase in ambient air temperature. 

The temperature coefficients of mean heat trans- 
mission coefficient of these relays is about + 0.3 %/*°C 
for winding temperatures of +50°C. The winding over- 
heating temperature of relays protected by thin covers 
having poor thermal contact with the windings, is de- 
creased by about a 1%/*C rise in ambient temperature. 


(12) 
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Examples 
1. We determine the winding overheating temper- 


ature of a small-size relay (without cover) for ambient 
temperature of +125°C if, for the same power and for 
ambient temperature of 425°C, the winding overheating 
was 50°C. 

We determine the relay winding’s overheating with 
ambient temperature of +125°C from (7) : 


S195 = 50 [1 — (1,94 + 0,0036 x 25) 10-9 (125 — 25)] = 40°C. 


2. We determine the mean heat transmission coef- 
ficient of a small-size relay for winding temperatures 
of 40 and of 100°C, if , for temperature 50°C, the heat 
transmission coefficient equals 1.5x 10° w/cm*/*c. 


We determine the relay’s heat transmission coef- 
ficient by (13); 


deo = 4.5 x 10-8 [1 + (3.45 —0,0079 x 
x 40) 10-% (40 — 50)] = 1.45 x 10-8 w/cm?.*c 


and 
Jioo = 1-5 X 1078 [1 + (3.45 — 0,0079 x 


x 100) 10-* (400 — 50)} = 1,7 x 10? w/cm**¢ 
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Equations are derived which characterize the basic modes and peculiarities of operation of a saturated-reactor 


Introduction 

















The operation of a magnetic amplifier loaded by a 
de motor (MA=M) has its peculfar features. They stem 
from the fact that the motor, as an energy consumer, 
presents itself as an impedance and a variable counter- 
emf whereby, with amplifier control, there fs an action 
on this emf. The presence of the amplifier leads, natur- 
ally, to a change in the drive's characteristics. In[1], 
individual questions in the analysis of such schemes were 
touched on for one special case. In [2], certain design 
approaches were adduced. At the same time, there is 
interest in establishing general regularities characterizing 
the operation of an MA=M system when its elements are 
connected in accordance with the scheme of Fig. 1. In 
the present paper we attempt to fill this void partially. 
We note that the analogous problem, as applied to a 
rotary amplifier, and also for the case when control of 
asynchronous motors is implemented by means of mag- 
netic amplifiers, received consideration in [3-5]. 























Excitation 





As the basis of our investigation, we take the linear- 
ization method described in[6, 7]. This method turns 
out to be particularly effective if one may make the 
two following assumptions; a) the characteristics of 
simultaneous core magnetization by constant and 
variable fields may be approximated by straight line 
segments ; b) the load curve, expressing the dependency 
between the variable component of the saturated reactor 


magnetic amplifier loaded by a dc drive. A simplified method of designing such an amplifier is given. 


emf and the load voltage can be replaced by a straight 
line. 

On the strength of the first assumption, we shall 
suppose that, for amplifier operation, we shall use a 
linear domain of the magnetization characteristic of the 
cores, which are made of material of sufficiently high 
quality. Under these conditions, the increment of the 
mean value of the working windings’ ac current AI,, 
can be found from the expression [6] 

AE, 
Al, =SAlo +R a) 


Here, ly is the control current, E,, is the mean value 
of the ac windings’ emf, S is the slope of the characteristic, 
and R, fs the amplifier's internal impedance. The last 
two parameters are defined as follows; 


S=c sh) : : (2) 
Op; we Q (3) 





R= ] 


Wg is the number of turns in the control winding, W,, is 
the number of turns in the ac winding, Q is the active 
cross section of the two cores, / is the mean length of 
0H OB 

PH, B= oH. , Hg is the 
constant (monotonic) component of the magnetic-field 
strength, H,, is the “mean” value of the variable com- 
ponent of the magnetic-fleld strength ; 


Wil 
H,=—-*, (4) 


the magnetic circuit, ¢= 





and B,, is the “mean” value of the variable component 

of the magnetic induction; 
E., 

5 

B,, = oW ., a) . ( ) 

The second assumption derives from the following 

considerations. It was established in [8] that, with am- 

plifier operation on the counter-emf, the load curve, as 

a function of material quality, can deviate from the 

straight line on both sides. 
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By taking this circumstance into account, as well 

as the fact that it is desirable to obtain simple graphic 
expressions, in the first instance mirroring the qualitative 
side of the phenomena, we shall assume that 

E., = U, ~ (E, + Rly) and, consequently, 


— AE,, = AE, + RAI. (6) 


Here , U, {is the mean value of the supply voltage, 
E, is the motor's counter-emf, R 1s the equivalent active 
impedance of the ac circuit. Figure 2 gives a geometric 
interpretation of these assumptions. The family of 
magnetization curves H,, = f,(Ho, By) is replaced by the 
dashed-line parallel segments. The linear domain is 
marked by the coordinates of points 1 and 2 (B;, H, and 
By, Hg). In correspondence with (6), the increments of 
voltage and emf dropsoccurring in the system are re- 
presented by segments (ab, bc, ce, etc.) whose projections 
on the axis of abscissas sum arithmetically to emf E,, - 


Hwy, ly, Io, M 






q 


sy 





Coefficients of Control and Rigidity 

We adopt the nomenclature that K, = 0n/0 ly for 
M = const is the control coefficient, and K, = -0M/dn 
for ly = const is the rigidity coefficient. Here, n is motor 
speed and M is the torque developed by the motor. 

We first consider the case when there is no feedback, 
and we assume that each conducting arm of the rectifier 
bridge can be replaced by a constant resistance r / 2. The 
back impedances will be considered infinite. 

If we take into account the generally known expres- 
sions for dc motors with independent excitation, 


M = Cul. 





(7) 


E. = Cen, (8) 





we can write (6) in the following way : 


— AE, = AnC, + RAM 
Cy 
Here, I, = l, is the motor armature’s current, R = 
=to +1_g+%%,1, 1s the impedance of the armature winding 
and motor brushes, r; is the impedance of the amplifiers 
ac winding. 
By substituting (7), (8), and (9) in (1), we obtain 


AM 
CyC ¢ 


(9) 


An= — AI, — 
E 





(Ri -+ R) = 


. i 
K. Al, = K Ad, (10) 
r 


where 


k= SR. 





It follows immediately from (9) that, all other 
things being equal, the range of control will be deter- 
mined by the limiting value of AE, , which we denote 
by AE%, - According to (5), this latter quantity de- 
pends on the properties of the core material (AB, in 
Fig. 2), its cross section,and the number of turns of the 
working winding. Since ABw is a constant for a given 
material, an increase of AE, may be attained by 
increasing Q and W,,- Simultaneously, as is clear from 
(3), Ry will also be increased, which implies an increase 
in K,, and a decrease in K,. In view of these circum- 
stances, itis frequently not possible to create an ampli- 
fler which, operating in the scheme under consider- 
ation, can satisfactorily provide the necessary values 
of control range and coefficient, and also the rigidity 











r's 











of the mechanical characteristics. We now show that 
a system in which the amplifier is used with feedback 
is free of such a disadvantage. 


Effect of Feedback on K, and K, 

a) Current feedback. Since, in addition to the 
ampere-turns Wy Aly, the ampere~turns WAI, (Fig. 1) 
also take part in core magnetization, we then obtain, by 
adding the proper term to (10), 

An = E- ( AZo + a Al< ) — Tez (fi + R), (11) 








where R=1,+1!,+% +1. By taking (7) into account, 
we obtain, after some simple transformations, 
An = - AlI,— toe [Ri(1— BS) + Rj. (12) 

Here, By = Wy / Wg is the feedback coefficient, W; 
and ry are, respectively, the number of turns and the 
impedance of the current feedback winding. 

The expression Rj (1~8 4S) can be considered as the 
equivalent internal impedance of the amplifier with 
positive current feedback, the magnitude of which can be 
varied within broad limits by the proper choice of 8; . 
This statement is well known from the general theory 
of amplifiers with feedback [9]. 

Thanks to this, we can obtain the necessary values 
for K including the value which provides (at least 
theoretically) absolute rigidity of the mechanical 
characteristics. The control coefficient does not depend 
on the current feedback. 

b) Voltage feedback . Change of the control coef- 
ficient can be obtained by introducing voltage feedback. 
For this purpose, we use a winding with number of turns 
W,, connected in parallel with the motor via resistor ry. 
The increment of feedback current will be 
Al, ‘ AE. t- Alar. 

re 
turns are AJ,W. = Wo8.(AE, +- Al,r,). Here, 
B, = re — is the feedback coefficient. 

The increment of total current passing through the 
amplifier’s working winding equals AL, = Al, + Al, - 
Magnetization will be implemented by the positive 
feedback of ampere-turns WpAly + W2AI,. Since the 
current Al, does not flow through the motor's armature 
winding, the component r,AlI, should be subtracted from 
the right member of (6), Le., 


— AK, = AE, 4- RAL +(R— r.) AlI,. (13) 


By substituting this last expression in (1) and by 
taking (7) and (8) into account, we find that 








and the corresponding ampere- 


me kAl, 
es Pe. 





AM [(R, + R)—(@ak —-2) rg] 





CeOu (1 — Bak + a) ‘ 





R,+R—r, 
6 as eee 


r? 


’ R=r; + Pa 4 Tc. 


Frequently, t¢( Byk-ay® Ry +R. Then, we can 
conclude from (14) that the introduction of positive 
voltage feedback leads to an increase in K, and a de~ 
crease in K.- The stability of motor speed will be even 
worse than with the absence of feedback.. With negative 
feedback, the inverse phenomena will be observed. 


, 1 
If | Bs |» co, then Ky — S" 


c 

c) Combinea Feedback. It is clear from an analysis 
of the two forms of feedback considered that their pro- 
perties complement each other to a certain extent, and 
that their joint use must give the very best results. This 
is an immediate consequence of the following expression, 
obtained on the basis of the same considerations and 
methods which we used previously: 





kAlo al 
e[t—e (m+n) +4] 


, Pe . 
AM {r, + R a k E (1 +e) +- te |+ arch 


Creu [1 — (8, P) a| aan 


Here ,R=%+l,+%+lo-ln deriving (15), it was 
taken into account that the increment of the total constant 
component of magnetizing force equals WeAly + W,Al, + 
+(W, + W,)Al,, and that AE,, is defined by (13). 

The relationships tp «rp, By<< t2B, ( By kar. << Ry + 
+R are frequently valid. In this case, (15) can be 
simplified to 


An : 





An = 





. kAl, ~ 
Cy (1 — Bak + @) (16) 


AM R,(4—BS)+R, 
CeCu (1 — Bok + a) 








By the proper choices of By and 8», one can attain 
the necessary values of both the control coefficient and 
the rigidity coefficient. From this point of view, com- 
bined feedback is the most advantageous. 


Transient Responses 

We first consider the case when there is no feedback. 
The law of variation of motor speed during the transient 
response may be found by jointly solving (9) and (10) 
(written for instantaneous values) and the following 
equations: 





Aug = roAig + WoQ £(nee) ’ (17) 


* The terms 6,/rz and t¢/t, will be missing if one 
connects winding W, as shown by the dashed line in 
Fig. 1, which is advantageous for B,< 0. 
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W ,VAB, = 


(18) 


d (An) 


AM = AM, + mAn -!- —— 2 (19) 


Here, uy» is the voltage applied to the input of the 
controlling winding, tp is the control circuit's impedance, 
By is the constant (monotonic) component of magnetic 
induction. Some explanation should be given a propos 
of the last equations. Expression (17) derives from 
Kirchhoff's second law for the control circuit; (18) was 
borrowed from [6, 7]. Henceforth, to simplify the pro- 
blem, bearing in mind that the cores are of good quality, 
we shall assume that W. = 0, and we shall denote 
vW,y /wWy by T. Equation (19) expresses Newton's law 
of motion with the assumption that the static impedance 
torque of the productive mechanism has two components, 
one of which is constant (AMg) and the second of which 
is directly proportional to the speed (mAn). On the basis 
of the equations enumerated, we get 














= of d? (An) , md ; 
leh AE + fea(t+ 2) + 
cs t 
(20) 
: 7 d (An)! m ) a AM, 
“mec | dt ait (i+ a An = K,Al, — Re 
Here, - 
oe K,TC, ; Le 
‘el fi Po ’ mec— K, , 
_ AU, _ RK; 
ere edna” 


Finding the solution of (20) poses ne difficulty. For 
a qualitative estimate of the possible results, we consider 
the following particular cases: Tey = 9, Tme-#0, and 




















Telt 0, Tiree = % Correspondingly , we obtain the 
equations 
eee aa in’: Kel, — 2% 

and ? 

salar d(An) | 

‘el ( K, * dt 

’ AM, 
+/(14 K- ) An = K,Aly— 


In the first case, the time to establish the number of 
revolutions is determined by the quantity 
Tmec J 


. eg K, +m : 
1+ K 
since, in the absence of feedback, K, is small. 

In the second case, the duration of the transient 
response depends principally on the amplifier lag, and 
will be characterized by the quantity 


, which will be quite large, 
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In investigating transient responses with feedback 
present, we should, instead of (10), use (16), adding to 
the latter the terms W,Qd(A By)/ dt and W,Qd( A By)/ dt, 
taking into account the appearance of the corresponding 
emf, in each of the feedback windings. It was proved in 
[6] that these emf's can be neglected in the case of 
current feedback and active load. We assume that this 
may be done in the present case, bearing in mind the 
qualitative nature of our estimate. Under these condi- 
tions we may, to elucidate the effect of the feedback, 
use (20), in which the values for K, and K; are chosen 
from (16). It is easily comprehended that an increase in 
K, due to positive current feedback or to negative voltage, 
feedback gives rise to a decrease in T,., and, conse~ 
quently, to an acceleration of the transient response. 
Positive voltage feedback increases K_ and, consequently, 
Tez as well. However, if ty is also simultaneously in - 
creased, maintaining the condition Usk, = const, then 
T ej Will bedecreased. Thus, with the assumptions we have 
made, all forms of feedback are advantageous, in the 
sense of accelerating the transient response. 


We note that if K, < 0, then Tp... < 0 and the 
term with d(An) / dt may become negative. In this case, 
the system exhibits a static instability , the assumption 
of linearity is unrealistic, and autooscillations will arise 
in the system, as we have observed in experiments. 


Several Remarks Regarding Nonlinearity of 
Mechanical Characteristics 

When low-quality cores are used, the two assump- 
tions made will not hold. In this case, it is impossible 
to consider the coefficients S and R; (1) as constants, 
where their maximum corresponds to a linear domain. 
With such cores, the load curve cannot be approximated 
by a straight line, as the second assumption required. In 
this case, a more suitable approximation is E’,, = 
= U’, “(E. +RI,,) , by virtue of which -<dE,,, = 
=(dE, +RdI,,,)(E, + RI,,)/ E,,- By denoting (E, +RI,,)/ 
/Ew=%,R;/ a =Riy and k/a = ky , we can find that 








dn = G, vo— a 


ia 





+ R) = Kdl, a (20) 
r 





Thus, K, and , will depend on the values of a, §, 
and Rj which, in their turn, are determined by the choict 
of the working point. The common structure of (10) and 
(10°) allows the basic propositions of the analysis providel 
earlier to be extended to the given case also, of course, 
in the neighborhoods of each working point 

We now explain the influence of S and R; on the 
linearity of the mechanical characteristics ; For this 
purpose, we determine the relative variation of the 
rigidity coefficient due to the variation of the internal 
impedance and of the slope,respectively on dR, and dS- 
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Using the expression for K, from (10), we find that 
dK dk, 
auptihe die b> Stee 
K, R, 
In the absence of feedback, the linearity of the mecha- 
nical characteristic will be determined by the constancy 
of R;- In the case of positive current feedback 








dK; (1 — BrS) dR, — BR's 
— Ridi—ps)+R ° 
If 
[45 (1— B.S) + RI - then 
: dK 
K. — oo and K-. 
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Increasing the rigidity by means of current feedback 
leads to a sharp manifestation of nonlinearity in the 
mechanical characteristics. Finally, as was shown earlier, 
with negative voltage feedback, Kr converges monoton- 
ically to its limit 1/r, if |Bg|+a@. In this case, 
dK;/ Ky will converge to zero. Consequently, negative 
voltage feedback allows the mechanical characteristic 
to be linearized and, therefore, is frequently used in 
combination with positive current feedback. This 
assertion is confirmed, in particular, by the data of the 
experimental investigation given in [2]. 


Design Considerations 

The order followed in designing amplifiers will be 
different, depending on the problem posed and the data 
originally given, which must be sufficient and which 
must be compatible. The consideration of all the ques- 
tions related to design goes beyond the limits of the 
present work. Here, we limit ourselves to illustrating how 
the expressions given earlier may be employed in design 
work. 

Let there be given the nominal parameters and 
operating mode of the motor, M*, I*., Cy, n® , the 
voltage on the brushes U* cr Tq and Cp , the limiting 
variation of torque and number of rotations AM® and 
An® , the core parameters 14,0, By, By, Hy,and Hg or 
the corresponding family of magnetization curves from 
which these parameters may be determined, and the 
necessary values of the coefficients K, and K, - 

We now determine the required volume of cores. 
This will be used in the best way if, in the process of 
Operation, there are attained the largest admissible 
values (in the limits of the linear domains) of magnetic 
induction Bz and magnetic-fleld strength Hy. By taking 
this circumstance into account, we find, on the basis of 


(5) and (13), that 


AE* + RAI +(R—r)y Ar 
Cc c  ——- 
oW,, AB, =@, 


= B,—B;. The components of the voltage 








where AB na 


drop given inthe left member of (21) are represented in 
Fig. 2 by the following segments; the increase in motor 
emf AE* 


, by segment ce; the voltage drop across all 








the een impedances due to the increase in drive current 
RAI..*, by segment bc; voltage drop due to the current 
of the voltage feedback (/? — r,) AJ, ~(R — re) J, — 
by segment ed, Since j° _ i + i (Fig. 2) then 
from (4), we find that 
Uo+I)w, 
It, 
It is assumed for this that the inequality (/c 
- Al:)W.> Hl always holds. 
Equations (21) and (22) permit the calculation of 
the active volume of the cores V = Q/ If one knows, in 
addition to the original data, the values of R and Iz , Le., 
essentially, rj, tg], andr,. Resistance r, depends on 
the rectifier properties and may be found by well-known 
methods. Ampiifier efficiency and volume of copper 
depend on resistances 1,1 , andt,. Their values can 
be varied within well-known limits. bt assume roughly 


a? (22) 





~ c 
that Dry ST, 2 = (5— 10) cae and, conse- 


quently, ,~=(0,2 te 0,1) k. We note that the core 
volume, found in accordance with (21) and (22) is 
minimal from the point of view of guaranteeing the 
magnetic mode. The possibilities of spacing the 
windings must be verified by additional well-known 
methods. 

By knowing the value of V, we can determine Q orl 
by giving either the cross section or the mean length 
of the magnetic circuit. For different ratios of Q/ J, we 
obtain different-size apertures and, to a certain degree, 
different operating conditions for the cores. The question 
as to the optimal ratio of Qf! requires a special investi- 
gation. As a first approximation, it is advantageous to 
start with the data given in the table below, obtained 
on the basis of the recommendations found in [10] : 





V/2,cm | heated 2.0 10 20 30 60 100 300 


+-——--- — 


fal ,cm \o.02 0.035 0.05 0. 075 0.1 0.15 0.2 0.3 

















By selecting the value of Q, we find from (21) the 
number of turns W,, , and from (3) the magnitude of Ry ° 
Then, using (15), we should determine 6,, 6,, and S 
for the given K, and K,. Since there are three un- 
knowns, one of them should be given, since the necessary 
control coefficient can be obtained by virtue either of 
the proper value of Wy or of W,. When internal feedback 
is used, 8yS 1. Then, 8, andS are determined uniquely. 
Finally, we calculate the supply voltage by the formula 
U, = oW.QB, + (I; —AIQR+ (Ue— AE), 
the meaning of which is clear from Fig. 2, on which the 
particular case is shown when AJ. ~ U7, 


SUMMARY 
The analysis of the equations obtained for the MA-M 
system leads to the following conclusions: 
1. When there is no feedback, it is not possible to 
provide a satisfactory operating mode of the MA~M 
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system. Introduction of current and voltage feedback can 


have an effect on the rigidity of the mechanical charac- 
teristics and on the control coefficient. The best results 

can be obtained with the concurrent employment of the 

aforementioned forms of feedback. 


2. Feedback permits the acceleration of the 


transient response and the linearization of the mechanical 
characteristics. 


3. The expressions obtained may be used for a first- 


approximation design of an amplifier in accordance with 
the methodology presented. 


1. 





LITERATURE CITED 
R. Kh. Bal*yan, “Push-pull magnetic amplifiers with 
de outputs,” Avtomatika i Telemekhanika 17 , 2 
(1956). 
A. B. Rosenbauli and I. N. Selivokhin, “Construction 
of the characteristics of a de drive with choke con- 
trol,” Elektrichestvo No.10 (1958). 
S. Ya. Dunaevskii, “On the range of speed control 
in a motor-generator system with a rotary amplifier,” 
Vestnik elektroprom No. 8 (1951). 
A. A. Sirotin, “An engineering method of calculating 
the characteristic of an independently excited motor 


8. 


10. 





in systems for the automatic maintenance of speed," 
Elektrichestvo No. 12 (1954). 

A. G. Ivakhnenko, Automatic Speed Control of Low 
Power Asynchronous Motors [in Russian] (Izd AN Ule 
SSR, 1953). 

L. V. Safris, "Designing a magnetic amplifier with 
linearized magnetization characteristics,” Trudy 
Rostov. Inst. Inzh. Zh.-d. Transporta No. 20 
(TransZhelkorizdat, 1956). 

L. V. Safris, “On the question of transient responses 
in magnetic amplifiers inductively loaded via 
rectifiers,” Avtomatika { Telemekhanika 19, 3 (1958), 
N. A. Kaluzhnikov, "On the design of choked mag- 
netic amplifiers, connected on a single-phase 
rectifier bridge,” Avtomatika { Telemekhanika 19, 
3 (1958). he 


G. S. Tsykin, Negative Feedback and Its Applica- 
tions {in Russian] (Svyaz*izdat, 1940). 


M. A. Rozenblat and O. A. Sedykh, “General prin- 
ciples in the construction production planning of 
toroidal cores for magnetic amplifiers," Otchet 
Instituta Avtomatiki i Telemekhaniki AN SSSR 
(1957). 








ith 











THE USE OF HALL PICKUPS IN MAGNETIC HEADS 
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This paper reports on the development and experimental investigation of magnetic reading heads based on the 
Hall effect. The heads have the property that their emf is proportional to the residual carrier flux, in contradistinc- 
tion to the usual induction-type reading heads in which the emf is proportional to the velocity of flux change. The 
heads can operate both with magnetic tape and with carriers in the form of drums or disks. 


In ordinary magnetic heads for reading or reproducing 
magnetic recording, there occurs a differentiation of the 
signal stored on a ferromagnetic carrier. With this, there 
appear distortions in the amplitude and phase frequency 
characteristics of the reproduction channel which hinder 
correct transmission of the signal. Exceptions are sinu- 
soidal signals which, in the given case, do not undergo 
changes in form if the system is linear. However, even 
for them there are serious limitations in the range of low 
audio and infraaudio frequencies. 

It is not always possible to tolerate this disadvantage 
of induction heads, despite their great superiority in other 
respects. Therefore, there were developed a number of 
designs of the so-called flux~sensitive heads of the servo 
system type, whose output emf*s are proportional, not to 
the velocity of change, but to the magnitude of the 
residual carrier flux. The best known of these are the 
magnetomodulation head whose operating principle is 
analogous to a magnetic amplifier (sonde) of the second 
harmonic type, and the electrobeam head, based on a 
tube with magnetic deflection of the beam. 

The flux-sensitive heads to be considered below, 
based on the Hall effect, are further developments in 
this direction. Their operating principle amounts to this, 
that the residual carrier flux, by means of a magnetic 
circuit, is applied to a Hall pickup, with the result that 
at the pickup’s output there is generated an emf which is 
proportional to the magnitude of this flux. 


One design of such a head is shownschematically on 
Fig. 1, a; an over-all view of it is shown in Fig, 1, b. 
Used in it is an ordinary ring-shaped magnetic circuit, for 
example, from standard head V-01. The dimensions of 
core 1 are as follows; length of the leading (reading) 
slit d; = 45 microns; length of the trailing slit d, = 0.28 
mm; outside diameter of the core D, = 25 mm; inside 
diameter D; = 18 mm. Height of the core assembly is 
h=7mm. The given design is intended to operate with 
magnetic tape 2 of width 6.35 mm, although it can also 
operate with a magnetic drum. The Hall pickup 3 is 
established in the leading slit of the magnetic circuit. 
Semiconductor pickups of n-type germanium (n-Ge) were 
used, with dimensions of 6x 3 x 0.18 mm. 

The Hall pickup has a de supply. The selection of 
the magnitude of the current turned out to have an 
essential influence on the head's signal-to-noise ratio. 
We used the following method to determine the optimum 
current magnitude. On the carrier there was recorded a 
sinusoidal signal with a frequency of 125 cps (carrier 
speed of motion was 38.1 cm/sec) and with a magnitude 
of residual flux @¢ = 300 millimaxwells. For the carrier 
used, this level was about 10 db lower than the saturation 
level. The pickup's Hall electrodes were supplied 
directly at the high-ohmic input of the amplifier circuit. 
Since the pickup's inherent noise is independent, within 
broad limits, of the magnitude of the signal, the signal-to- 
noise ratio was determined as the ratio of the indication of 


























the instrument at the output of the amplifier circuit at 
the time of carrier motion to the instrument's indication 
for a motionless carrier. This ratio was taken off as a 
function of the magnitude of the current supplied to the 
pickup. As a result, it was found that the optimum 
magnitude of this current for the head under considera- 
tion was 8 ma. With such a current, and for an amplifier 
circuit passband of 20 to 6000 cps, the signal-to-noise ratio 
was 20.4 db; for a passband of 75 to 150 cps, this ratio 
was 33.2 db and, for a passband of 25 to 100,000 cps, it 
equaled 17 db. To eliminate pickup from the constant 
poles, the current supply circuit for the pickup was 
carefully screened. A thermal regimen corresponding to 
the pickup’s optimal supply current and to an ambient 
temperature on the order of 25°C provided prolonged 
stable operation of the head. For example, the pickup 
was supplied with current continuously for 48 hours. With 
this, the operation of the head was periodically verified. 
No operating changes were observed. 

The magnitude of the emf generated by the pickup 
when it is acted upon by the carrier's magnetic field fs 
related to the head's sensitivity. As a preliminary 
measure , we found the character of the calculated 
dependency of the head's sensitivity on the parameters 
of its magnetic circuit and on the so-called voltage 
sensitivity of the Hall pickup [1)}. 
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Fig. 2. Simplified equivalent 
schematic of the head's magnetic 
circuit. Rz is the reluctance of the 
leading (reading)slit,R, is the re- 
luctance of the trailing slit,and 
%c 1s the carrier's residual flux. 


In a first approximation, the equivalent schematic 
of the head's magnetic circuit has the form shown in 
Fig. 2. The emf at the pickup's output may be given as 


E = B, (1) 


where y is the Hall pickup’s voltage sensitivity (micro- 
volts / gauss) and B is the induction acting on the pickup. 
In our case, B is the average value of the induction in 
the trailing slit due to the residual carrier flux $c. 

We determine the head's sensitivity from the 
formula 


oe Ae I Ring 
ee —— (2) 


where v is a coefficient defined by the parameters of 
the head's magnetic circuit; 
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v= 7 = i 
©.S5, 5, +a/R;' 


S, 1s the cross section of the flux in the trailing slit, Ry 
{is the reluctance of the leading slit, i.e., 


7 
1=3-yaR, 
S,+a/R, 





In the general case, ®c is, with certain limitations, 
an arbitrary function of the space coordinates related to 
the carrier, and B and E are functions of time correspond- 
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Fig. 3. Scheme for measuring the Hall 
pickup’s sensitivity, with the pickup lo- 
cated in its operating position between the 
halves of the head's core. 1 is a generator, 
2 is a vacuum tube voltmeter, 3 is the 
head's core,and 4 is the Hall pickup. 


The pickup's voltage sensitivity was determined 
experimentally as shown in Fig. 3. The head's core 
had two windings: magnetizing winding Wy, and measur- 
ing winding W, The measure of magnetic induction in 
the core was taken to be the magnitude of the voltage 
induced in winding W,. The ratio of the mneasured 
magnitude of the Hall emf E to the corresponding value 
of induction B gave the mean value of the pickup'’s 
voltage sensitivity y = 116 microvolts/gauss from whence, 
for the parameters of the magnetic circuit previously 
enumerated, we found the computed sensitivity of the 
head: 
4 = 204 pv/ usec 
The experimental value of the head's sensitivity 
was determined by measuring the emf at the pickup's 
output and seeking the dependency of the emf on the 
magnitude of the residual flux for an 8 ma supply current. 
The relationship obtained was close to a linear one. It 
was found that, for the aforementioned intensity of 
magnetization, the Hall emf (the equivalent value) at 
the head’s output was E = 52 microvolts. From 
this was obtained the experimental value of the head's 
sensitivity: 
52 
Nexp ~ 03 
The magnitude of induction B in the trailing slit 
can be found from (2): 


= 173 uv/ usec 





B= + On = 0.45 gauss 
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On the basis of the dynamic range of 20.4 db, found 
earlier, the head's sensitivity threshold corresponds to a 
yalue of induction of 


0,45 


Bmin= Jo-we0 = 4.3-107* gauss 


The head's frequency characteristic, obviously, is 
determined by the frequency characteristics of the two 
factors in (2). The frequency characteristic of the 
voltage sensitivity y was measured by the circuit of Fig, 3 
for a constant value of induction B = 3 gauss within the 
frequency limits of 20 to 200,000 cps. It was found 
that, within these limits, the voltage sensitivity does not 
depend on frequency. 

In a pickup placed in a variable magnetic field, 
there is induced an induction emf which is added geo- 
metrically to the Hall emf. However, the induced emf 
can be easily compensated by the proper placement on 
the core of one of the conductors joining the Hall 
electrodes to the amplifier circuit. The optimum com- 
pensation is chosen as that giving the minimum deviation 
of the instrument in the circuit of the Hall electrodes at 
the highest operating frequency. On the basis of the 
system's linearity and the applicability of the principle 
of superposition, one can use still another method of 
isolating the Hall emf from the geometric sum of the 
two emf's ; with the current supplying the pickup 
switched in, one measures (without compensation) the 
total emf, after which, with the current switched off, one 
measures the induced emf. Their geometric difference 
gives the Hall emf. By such a method, results were 
obtained which coincided closely with the previous ones 





until the induced emf exceeded the Hall emf by a factor 
of more than 3 or 4. This latter method is effective 
when the Hall pickups are used as pickups for geometric 
sums or differences of stimuli, one of which is propor- 
tional to the frequency. 


The frequency characteristic of coefficient v is 
determined by the magnetic and geometric parameters 


of the core such as, for example, the slit and the frequen- 


cy-dependent losses, and is not characteristic of the 
given head. 


Since, in the range considered, the pickup sensitivity 
does not depend on frequency, the head's frequency 
characteristic is simultaneously the frequency characteris- 
tic of the induction in the trailing slit. 


The head's experimentally obtained frequency 
characteristic is shown in Fig. 4. The signal was recorded 
at a constant writing current but at different frequencies, 
with a speed of carrier motion of 38.1 cm/sec. Some 
lifting of the characteristic in the low-frequency range 
is related to the length of the head contact with the 
carrier. The fall of the characteristic in the high-fre- 
quency range is explained both by slit phenomena and 
by certain other peculiar characteristics of magnetic 
recording. 


On the other hand, since the Hall emf may be 
defined in practice as the geometric difference of the 





total emf and the induced emf, this is confirmed by the 
form of the proportional dependency (1) between the Hall 
emf E and the induction B, {.e., the head possesses the 
property of flux sensitivity. 

This latter was illustrated by the following experi- 
ments. Recording was read from an immovable carrier. 


The magnetization intensity of the carrier corresponded 
to; 


1) a sinusoidal signal with a recording wavelength 
of 16 mm; 


2) a pulse signal with a dipole recording length of 
10 mm (recording from zero to saturation; a mark-space 
ratio of 2); 

3) a rectangular signal with a dipole recording 
length also of 10 mm (recording from saturation to 
saturation; mark-space ratio of 2). 

The carrier, by means of a micrometric device, 
was drawn past the head for a definite distance and 
then stopped, after which a sensitive galvanometer 
measured the voltage on the head's Hall electrodes. 
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The point of reading initiation coincided with the 
beginning of the recording, for which the initial portion 
of the recording was made visible by precipitating on the 
carrier a powder of carbonyl iron from a heptane suspen- 
sion. As a result of this, there was shown the synphasing 
of the carrier's surface residual flux and the voltage on 
the head's Hall electrodes. 





whose length is about 0.35 mm. The design has extended 
pole shoes (the length of each shoe is 25 mm) and is 
intended to read recording of long wavelengths. A sim- 
plified equivalent schematic of the head's magnetic 
circuit can be obtained if, in Fig. 2, one sets R; = @. 
As should be expected, with the same pickup, this design 
has greater sensitivity and a correspondingly better 
dynamic range than the first design. At a frequency of 
40 cps, the increase in sensitivity is about 4 db. 

Other heads based on the Hall effect which were 
investigated differed from the two types considered here 
only by the core configurations and materials; for 
example, we used ferrite cores whichpermitted us to do 
without insulation of the pickup when it was placed in 
the slit, thereby decreasing the length of the slit and 
increasing the head's sensitivity. 


SUMMARY 
A head based on the Hall effect possesses the fol- 
lowing advantages over the most widely disseminated 
type of flux-sensitive transformer—the magnetomodula- 
tion head. 


1. It permits simplification of the electrical reading 
scheme. 

2. It possesses sensitivity to the direction of the 
carrier field, i.e., to the polarity of the signal. 

3. It can operate in a broad frequency range, in 
particular, for frequencies of the order of 100 kcps and 
higher. In this case it is advantageous to use a ferrite 
magnetic circuit. 

4. With heads based on the Hall effect, the danger 
of random erasure or magnetization of the carrier is 
excluded. 

5. Phase distortion of the signal is virtually elimi- 
nated. 
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Fig. 5. 1 is the core, 2 is the carrier, 3 is the pickup. 





We should mention one further trial head design, 
shown in Fig. 5, a and b. The head's magnetic circuit 
is open, and the Hall pickup is placed in the reading s}jt 









































6. The requirements on head shielding are reduced, 

7. Simplicity of design and small dimensions are 
provided. 

This last quality effectively distinguishes the head 
based on the Hall effect from the other well-known 
flux-sensitive head—the head of the electron-beam type. 

The Hall pickup is an inductionless transformer, 
which eases the shielding of such a head in comparision 
with an ordinary induction head of equal sensitivity. 

A head based on the Hall effect may be useful in 
many cases of special uses of magnetic recording, the 
importance of which is continuously growing; for 
example, in various devices of automation and remote 
control, where it is required to read pulses with minimal 
frequency and phase distortions; in geology , particularly 
in oil prospecting, for reading processes which change 
at infraaudio frequencies; in medicine and biology, for 
studying slowly varying processes related to brain acti- 
vity, heart activity, etc., in particular, for reading 
magnetic cardiograms. 

























































































The most efficient use of heads based on the Hall 
effect is in the frequency range from fractions of a cps 
to about 100 cps, and for carrier speeds up to 4 to 5 
cm/sec, i.e., in those cases where induction heads are 
practically ineffectual. The basic disadvantage of a 
head with a Hall pickup for sound reproduction is the re- 
latively small dynamic range. There exist the following 
ways of improving the dynamic range and other indicaton 
of head quality: 










1. Improvement of the technology of preparing Hall 
pickups, the quality of processing the surfaces and contac 







2. The use of Hall pickups with non-current-carryil} 
contacts. 
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3. The use, for pickups, of materials with higher 6. Operation with high levels of carrier magneti- 

Hall constants. zation, for example, with levels corresponding to the 
4. Reduction of the length of the slit in which the magnetic saturation of the carrier, as occurs in the 

Hall pickup 1s established, which introduces the particular ™ajority of special applications of magnetic recording. 


requirement to reduce the width of the pickup. There is tt t , th 
interest, in this connection, in film pickups of materials in conjunction with this, it should be mentioned that, 


with high voltage sensitivity. These latter can be mounted aoe jena Betig 2  he npahtong onty 
in the form of a thin film of thickness of 10 to 100 micro- 
meters on the face surface of a ferrite core. LITERATURE CITED 

5. The use of magnetic recording carriers with high 1. A.R. Regel’ , Semiconductor Measurers of Magnetic- 
values of residual flux. Field Strength{in Russian] (LDNTP, Leningrad, 1956). 
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A description is given of a voltage deflection pickup implemented by semiconductor elements — silicon base 


diodes, transistors, and thermoimpedances. A design methodology is proposed, and calculated relationships and ex- 


perimental results are given. 


General Remarks 

The voltage deflection pickup considered in this 
paper is designed for use in systems for stabilizing 
voltage or any quantity which is transformed to a voltage. 
The special feature of the voltage deflection pickup, 
shown in Fig. 1 and 2, is its implementation by semi- 
conductor elements. This provides its small weight and 
size (Fig. 2), reliability, heightened resistance to 
vibration and shock, and length of useful life . 

Type D 810 silicon base diodes D; and Dj, con- 
nected in opposite arms of bridge Dy-R_~Ryy-Dj “Ry, 
execute the function of voltage calibration. The voltage 
to be regulated, E,,.,. 1s applied to one of the bridge's 
diagonals via ballasr resistor Ry. Connected to the 
other diagonal of the bridge is a preamplifier of silicon 
transistors T,-P102, T}-P102, whose output current con- 
trols the relay power amplifier, implemented by german- 








ium transistors T,, Ty-P13A, Ts and T';-P4D. Thermo- 
impedances Ry ; and RT, are used to compensate the tem- 
perature variations in the scheme. 

To obtain, in load impedances R, and R, , a signal 
whose average value would be a continuous function of 
the input signal for a limited range of variation of the 
latter, one should connect, instead of resistors Rs and 
R's, condensers C, and C", (shown by dashed lines in 
Fig. 1),which would provide self-excitation, i.e., the 
principle of vibration linearization is used. Pulse~width 
modulation is used in the scheme for deviations of the 
controlled voltage. 








* The results of the investigation described in this paper 
were discussed in the April, 1959 seminar of the Laboratoy 
of Automation and Remote Control Elements of the IAT 


AN SSSR . 





Fig. 1. Simplified schematic of the voltage deflection pickup . Ry = R'; = 1.2 kohm; 


Re = R'p = 3.3 kohm; Ry = R's = 100 ohm; Ry = 
= § ohm; R, = 5 ohm; Rg = R's = 1 kohm; Rg = 0 ; Ry = 200 ohm; Ry = 1920 ohm; Ryy= 
= 2100 ohm; R,, = 542 ohm; Ra i-ap = 210 ohm; Ra, = 40 ohm ; Rg = 395 ohm. 
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', = 5.6 kohm; Rg = R's = 15 kohm; Rg = 
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Posing of the Problem 








The problem raised in this paper is the presentation 
of an engineering method of designing the voltage 
pickup. Relationships are derived which, for given 
regulated voltage Ereg, source voltage E, and load 
resistances R, and R';, allow one to determine the re- 
maining circuit parameters for which feasibility of the 
circuit and its stability under temperature variations 
would be guaranteed. 

In addition, a generalized method is propounded for 
the analysis of transistor switching circuits. 


Method of Investigation 

Thanks to the fact that silicon transistors in the cut- 
off region have insignificant current flowing through the 
base-collector junction (I,g = 1 4),it is possible to 
design the bridge together with the thermoimpedances 
Ry, and Rr, , and to implement the additional impe- 
dances Rg and Ry» separately. 





In designing the bridge, the nonlinear base diodes 
D, and Dj , after the piecewise-linear approximation of 
their volt-ampere characteristics, are given as “condi- 
tionally autonomous" dipoles with emf's E, and E° , and 
with impedances R, and R", (Fig. 3). 

The high-ohmic thermoimpedances Ry . and Rr, 
are presented as current sources Ip, and Ir, . 

As a result, we obtain the replacement scheme, 
shown in Fig. 4, for the bridge. It is calculated by the 
methods of the linear theory of electrical circuits. 

By considering the outputs of transistors T; and T‘; as 
current sources (in the equilibrium state the currents 
virtually equal zero), one can compute the relay stages 
of transistors T,, T’,, Ty and T'; separately. The design 
of the latter is implemented from the conditions of 
transistor cutoff and saturation, with account taken of 
temperature variations. 

For the computation of the circuit's sensitivity (cf. 
the Appendix), we introduce the concept of a “candi- 


Fig. 2. Outside view of the voltage deflection pickup. 
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Fig. 3. Volt-ampere characteristics for 
determining parameters. a is for the 
silicon base diodes, b is for the silicon 
transistors . 














Fig. 4. Equivalent replacement scheme of 
the nonlinear bridge. 


tionally autonomous three-pole? allowing one to take 
into account the nonlinear properties of the transistor 
and, at the same time, to use a generalized method of 
designing linear electrical circuits. 


Designing the Nonlinear Bridge 





Use of type D808-813 base diodes shows that, for 
the greatest stability, their operating current ly must be 
of the order of 5 ma. 











The magnitude of currents I... and Iy, on the circuit 
of Fig. 4 are determined, with sufActent accuracy, from 
the equation 
E 

l= zi, (1) 
where E, is the emf defined by the diode characteristics 
(Fig. 3), Ey < E', ; Ry is the magnitude of the thermo- 
resistor (the uheotiet is omitted, since the equation is 
valid for both diodes and thermoresistors). 

The thermoresistors are chosen large enough so that, 
first of all, operation of the thermoresistors on the linear 
portion of the volt-ampere characteristic is provided [1] 
and, secondly, the thermoresistors Rr, and Ry can be 
replaced by current sources. This latter is realized for 
the condition 


To > 101. (2) 


By neglecting the components of I-7 in the bridge, 
we can find the resistances Ryg, Ry and R = Ry + Ra + Rajt 
+Ra, (Rp =R,+ Ra; - Ra, , Fig. 1) from the equations 





; E’. (3) 
Rie= Rat 7 , 
; E. E 
R= 2 (R444) —A,— 48. (4) 


R= reg? mene Pd) 
alo 





(5) 


where R', is the diode's dynamic impedance. 

The remaining unknown resistances, Ra» Ray Rag 
and Ra, are determined from the conditions of equi- 
librium of the circuit far normal and for elevated 
temperatures. For this, it is necessary to determine the 
effect of temperature variations of the characteristics 
of the diodes and thermoresistors on the displacement 
of points ay and a, with respect to point b (Fig. 4). 

For Rq « R and R'g & Rg, the temperature variation 
of the diodes’ characteristic is equivalent to the intro- 
duction of a voltage source (Fig. 3) 

(6) 
ed= ed, + ed, = AE, + Tyra, 
where 


Iy=Ip+Ir+ir, eg,= AE, ra=ARg ip = Aly. 


Using the principle of superposition, we determine 
the effect of temperature variations of the diodes’ 
characteristic from the matrix equation 


A ne Rs — Ry ig be (7) 
in| —Ry Ris + Re + ed}’ 
where 1, and 1, are the increments of currents I, and I, 
(Fig. 4). 

Solving (7), we get 


(Ris + Ris) ey — Riseg 


(8) 
R (Ris + Ris) + Rishi 





i, = Al, = 





By means of (8), we determine the displacements of 
points ay and a,; 
* AUba, = ed + i: (R — Ra— Ra), 
AUba, = eat i, (R— Ra). 
Analogously, we find the displacements due to 
changes in current I, with the condition that R * R 


(variation in IT, has no practical effect on the displace- 
ments of ay and ay) : 


R(R, + R,,) (Ria + Ris) , ; 






(9) 





















AU ja, = — 
be R (Riz + Ris) +- Riokis ty 
a RRiwkis—R,R(Riz+ Ris), 
6 a R (Ria + Ris) + Risks tn (20) 






By taking (8)-(10) into account, one can set up the 
four equilibrium equations for the two temperatures, by 
solving which we obtain 
I, (AUep, + AUep,) + és, Yen, + Yor) 

T,,41 + Loi,, . 
(Io —1,,) (AUeb, + Ay) — lin + i,,) Yer, + Ye) 


loa so Toi, 
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(rio — Poly.) ig, + (yh + Poh) I, 









Ry oil (r3/,, — rel) U, “shZ (rvis + rots) U, (11) 








(rilo—r3l,,) i, (r2is+r3i,) ‘. 





where 





r? = R(Ry + Ris) + Riki, 
r2 = R(Ry. + Ris), 







U, = 1,R — Ey —Uey, — Ih, a Ra, 


1 






Ri2Ris 
4 
Ver, > U'eb, » AVen, and AU",), are the threshold 
voltages and their increments , defined for the cor- 
responding triodes T, and T‘,; from Fig. 3, b. 

As is obvious from (11), for certain special values of 
the parameters on theright, the magnitudes of resisters 
Rg, Ra, and Ry may turn out to be zero. 








Uy = —i,R—e4+AU,, — in, Ras 








Design of Relay Stages (Switching Circults) 
The design of relay stages from transistors T,, T's, 
T, and T', leads to the determination of resistances 
Re = R's, Rg = Rg, Rg = R's . Rg, Ry, and Rg = R', for known 
load resistances R, = R'; and supply voltage E,- The six 
unknown resistances are found from the cutoff and 
saturation conditions of transistors T,, T';, T; and T's 
[2] and from two auxiliary relationships, the first of 
which expresses the condition that there be no shunting 
of the relay amplifier’s input 














(12) 
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Rj, can be set at 500 ohm for P13A), the second of which 
is the condition of insignificant losses in bias resistors 
Rg and Ry: 
10(R, + Rk) < Rs. (13) 
The cutoff condition for transistors T, and T', is 
written in the form 


leodte S ERR’ 


(14) 


where Ico, is the temperature collector current of tran- 
sistor T,- 


The cutoff condition for transistors T,; and T's, 
neglecting the emitter current Ie, of triode T, in resistors 
Rg and Ry , is defined by 

3+ Ret RB; 

The saturation condition for triodes T, and T'», 

neglecting Ico, and le, in Ry, is written 


Po (3 _— ris) Toon 
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RoE < RR, + Rg (Rs + Re) — 
R, , Rs+Re\ 
— Ro + eR): (16) 


where 
ro= RR + A+R, 
= RR, + RR, + RRs, 
rg = R.R;Rg, 
8, is the current transmission coefficient of the base of 


transistor T,, and I,, is the saturation current of the 
collector of transistor T>. 


The saturation condition for transistors T, and T} 
is defined by the inequality 


KR 
By (Rat Re Ry ~ £0. Sb (17) 


is the base current of transistor T,. We give this 
condition in expanded form; 





where I 


E \ 
i (= cat cova } SE, (Ass + Relts) — 


— I _ RR,fre.- 
CO, 2“*5’ a (17°) 
By taking , in (14)-(17), the value of I. which cor- 


responds to the heightened temperature, we can arrive 


at equalities. A huge increase in the right members of 






= 24+0.15v. 


Fig. 5. Oscillograms of pickup output voltages with self-excitation. 
a is for E.., = 240.15 v, b is for Ereg * 24 v,andc is for Ereg = 





(14)=-(17) leads to a lowering of the circuit's sensitivity 
(cf. the Appendix). 

When the pickup is used in the autooscillatory mode, 
the half-period of the autooscillation, with the conditions 


Ru=Ru=, R;,=0, R= R; = ce 
may be found approximately from the formula 


CEs RoR, + RoR, + ahs) (18) 
(Ry + Ry) (Rs + Re) : 
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where Ico, is the collector current of transistor T,(T*;) 
at equilibrium, the deviations from which are approx- 
imately proportional to the relatively small deviations 
of the regulated voltage ELeg" 


Experimental Investigation 

The voltage-deviation pickup, designed in accordance 
with the formulas given above for the case of E.., = 24 v 
and E, = 27 v, was subjected to experimental investiga~ 
tion. The pickup's basic parameters are given in the 
caption of Fig. 1. Testing of the pickup in the tempera- 
ture range from +20 to +60°C showed that it reacted 
reliably to deviations of + 1 volt. The output current 
in loads Rg and R'; in the cutoff state for a temperature 
of + 60 °C did not exceed 1 ma. The circuit operated 
reliably on an inductive load (a differential electro- 
magnetic mechanism). With output stages of transistors 
P4D, the circuit was capable of a switching power of 200 w. 
The results of an investigation of the pickup in the pro- 
portional mode for a self-excitation frequency of 50 cps 
are shown in Fig. 5. 





Conclusions 

The theoretical and experimental investigation 
provided gives a foundation for concluding that the voltage 
deviation pickup constructed of semiconducting elements 
was efficiently developed. The essential advantages of 
the pickup described —its small weight and dimensions, 
the increased resistance to shock and vibration, the long 
useful service life, inherent in semiconducting 
elements--permit the expanding of the domain of 
applicability of voltage deviation pickups. 

The author wishes to express his gratitude to B. S. 
Sotskov for his valuable advice on the execution of the 
present work. 











APPENDIX 
Calculating the Sensitivity of a Relay 
Schenre 

The sensitivity of a relay circuit is defined by the 
threshold values of controlling current for which the 
circuit transfers (operates). 

In relay schemes constructed of transistors, the 
currents and voltages on the electrodes of these latter 
vary within such wide limits that there appear essential 
nonlinearities in the transistors’ properties (cutoff region, 
saturation and nonlinear active region). However, for 
a broad class of transistor relay schemes, for definite 
threshold values of controlling current, it is not neces- 
sary to know in detail the nonlinear characteristics of 
the transistors. It suffices to know the characteristics in 
the neighborhoods of isolated points, corresponding to 
the moments of circuit transition points on the boundaries 
between the cutoff and saturation regions on the one hand, 
and the active region, on the other. 

A nonlinear transistor can be presented as a linear 
“conditionally autonomous” three-pole (a four-pole with 
a short-circuited side) for the regions of cutoff and satura- 
tion, and also for two (or more) segments of the active re- 
gion, for example, for the initial (close to the cutoff re- 
gion) and for the final (close to the saturation region) seg- 
ments. 

Indeed, the behavior of a transistor for constant 
current in any form of connection (with common emitter, 
base, collector) is completely described by two nonlinear 
equations. 





For the correspondingly chosen dependent and indepen- 


dent variables, there are six possible variants (systems) of 
equations. For example, for a circuit with a common 
emitter, the equations for the dependent base and 
collector currents have the form 


1 =Y¥i (Upp: Use 1p = ¥2 (Ugh, U,,)» 


where I}, is the base current, I, is the collector current 
Ugp 1s the emitter~base voltage and U, , is the emitter- 
collector voltage. 

On any of the aforementioned segments (in the 
neighborhood of a point,), one can expand both functions 
in a Taylor series for two variables, and drop the terms 
for higher-order derivatives. For example , for the 
initial portion of the active domain, 


OY, oY; 
1h =—Gteb Yer—Yenr > + S¥ee ee — Yect ) + 





oY 
U, I= s — 
+ ¥i(9 247, Yect 4, aU Yen ) + 
of U, U 
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or 

Ih = Vy Yep + Vig 7 ect Yol ; 

Fe = Va Ven + Yeo Ue, + Yai: 
where Y4I , Yyot, Yogt » and Yg97 are the nonautonomous 
transistor parameters forthe transistor as a four-pole 
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system of conductances for the initial active segment, 


Yior=¥1(Ueyy Yeoy —WnqVeb 1 + Yuer%er ): 
Yoon =¥2 eng Vect)—Mapeny + Yerq%, 1 


are the “conditionally autonomous” parameters of the 
transistor as a four-pole for the initial active segment. 

By using the well-known relationships for nonauto- 
nomous and autonomous parameters [3], one can rewrite 
the parameters obtained in any system of parameters and 
for any scheme for connecting the transistor. The para- 
meters for the remaining segments are determined 
analogously. 

A “conditionally autonomous” three-pole can be pre- 
sented in the form of a nonautonomous three-pole (core) 
with external autonomous dipoles. Such a representation 
permits the use of the well-developed methods of analy- 
zing schemes with nonautonomous multipoles [4]. 

By experimental measurements or by determining 
from static characteristics the aforementioned transistor 
parameters, one can replace the circuit of Fig. 1 by an 
equivalent linear replacement circuit including , together 
with dipoles, multipole elements and sources. The order 
of circuit design is as follows: For the first step, one 
determines the transmission coefficient of the portion of 
the circuit from Ereg to the output of the preamplifier 
R, and R', may be grounded). The second step consists of the 
determination of the value of collector current of triode 
T,(T';) which is necessary to operate the relay output 
stages. For this, the transistors are replaced by “con- 
ditionally autonomous” three-poles. As a preliminary 
step, one determines the collectors’ saturation current, 
and removes outside the three-pole core the autonomous 
two-poles. Thereafter, the replacement circuit is 
constructed. 

To determine one (negative) threshold value of the 
current applied to the input of triode T, , it is sufficient 
to determine the magnitude of it for which the base 
current of cutoff triode T, becomes equal to the boundary 
base current at the beginning of the active portion; the 
value of this last current is determined from the triode's 
forw ard transmission characteristic. 

To determine the second (positive) threshold value 
of current applied to the input of triode Ty, it is first 
necessary to find the value of base current of triode T; 
for which saturation of triode T', occurs. For this, after 
disconnecting and grounding the left end of feedback 
resistor Rs (Fig. 1), one sets up the matrix of conductance 
[Y] (or of impedances [Z] ) by the method presented in 
[4]. By a well-known formula, one determines the 
current transmission factor from the base of T, to loadRy: 


Ais 


kK, = ——_—__ 
An +R, A 


where A is the determinant of matrix [Y], and Ag, and 
Age ate cofactors of the determinant (subscript 1 refers 
to the input (base of T,) and subscript 2 refers to the 
output J. After this, one can determine the base current 








of triode T', which corresponds to saturation of T's, and 
thereafter one finds the threshold value of the input 
current. 

In our case, the problem is essentially simplified, 
since the circuit is symmetric and 1s controlled only by 
negative values of the input current, {i.e., only one value 
of threshold current is determined. 

We note that the method presented permits one to 
take into account the effect on relay sensitivity, not 
only of the transistors’ nonlinearities, but also the 
temperature variations of their parameters. For this, it is 
necessary to know the transistor characteristics in the 
operating temperature range. 
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REGARDING THE THEORY OF THIRD-ORDER LINEAR SYSTEMS 


Yu, G. Zarenin 
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Translated from Avtomatika 1 Telemekhanika, Vol. 21, No. 3, pp. 417-419, March,1960 


Original article submitted September 28, 1959 


The frequency properties of linear third-order systems are considered. The possible types of frequency charac- 
teristics are defined and the relationships which uniquely determine their form are found. The results of the investi- 


gation are illustrated by diagrams. 


Many automatic control systems can be described 
by linear third-order differential equations with constant 
coefficients. In spite of this, the theory of third-order 
systems is still not completely worked out. In particular, 
there has not yet been obtained an exhaustive solution to 
the problem of investigating the frequency properties of 
these systems. The literature is lacking in statements 
concerning the possible types of amplitude-frequency 
characteristics for these systems, as well as any informa- 
tion as to how one might determine the exact number of 
extremal points of these characteristics (not to mention 
their coordinates), if the coefficients of the system's 
differential equation are known. Here, we attempt to 
fill this void. 

The investigation of the amplitude-frequency 
characteristics of third-order systems is divided into 
three stages. 

a) Determine the possible types of these characteris- 
tics. As the basis for classifying amplitude-frequency 
characteristics, it is convenient to take the number of 
their extremal points. 

b) Elucidate the relationships by means of which one 
might analytically determine, from the form of the sys- 
tem's differential equation, the type of its frequency 
characteristic. 

c) Define the abscissas of the extremal points of the 
system's amplitude-frequency characteristic in terms of 
the coefficients of its differential equation. 

We write the differential equation of a third-order 
linear system in the form 

(Ap + Bp? + Cp+D)9=4, @) 
where p is the differentiation operator, g is the instan- 
taneous value of some physical quantity at the system's 
output which is of interest to us, # is the instantaneous 
value of the signal applied to the input and A, B, C, D 
are constant coefficients. 

To obtain the simplest possible relationships in the 
sequel , it is convenient [1] to divide both sides of (1) 
by the coefficient of the highest-order derivative. 
Denoting the new coefficient thereby obtained by 
a,b,c ,andd ,wecan write 
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(p? + ap? + bp + c)e = dy. (2) 


With this, the equation of the amplitude-frequency 
characteristic has the form 
d. 


Gas | 
| Il (jw) | V (c _— aw")? 4 (bw — w)? : (3) 


K (w) 

The number of its extremal points equals the 
number of positive real roots of the first derivative with 
respect to w of the radicand in the denominator of (3). 
After a single differentiation of this expression, we ob- 
tain a fifth-order polynomial which contains no free term; 
corresponding to this one may assert that the first ex- 
tremum of a third-order system's frequency characteristic 
lies on the axis of ordinates. In addition to the root 
w = 0, the polynomial under consideration has four pair- 
wise conjugate roots. By considering all these roots to 
be real, we arrive at the conclusion that the maximum 
possible number of extremal points of the system's 
amplitude-frequency characteristic equals three. The 
minimum number of extrema, realized in the case when 
all four of the aforementioned roots are imaginary, 
equals one. 

Thus, for a third-order system only three types of 
frequency characteristics are possible (Fig. 1). The 
dashed lines show frequency characteristics which are, 
hypothetically, theoretically possible but which, in prac- 
tice, are not realizable. 

Let us consider in more detail the roots of the first 
derivative of the radicand expression. By differentiating 
this expression, and introducing the notation 


b? — 2ac = Nan@b -- a? = R, (4) 
we obtain the polynomial 
30> — 2Ru* + Nw = 0, (5) 


the roots of which equal 


1 er 
w, = 0, o545 = Va VR +V hh? — 3N. (6) 
Discarding the roots less than zero, we obtain the 
abscissas of the extremal points of the system's amplitude- 
frequency characteristic; 
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The following conclusions derive from Expressions 
(1). 

1. If N < 0 then, independently of the sign of R, 
frequency w, is imaginary and wy fs real. The system's 
amplitude-frequency characteristic is a curve with two 
extrema—at the frequencies w = w, = 0 and w = wy 
(Fig. 1, b). 

2. ForN> 0 andR< 0, both frequencies w, andwy, 
are imaginary. The amplitude-frequency characteristic 
has only one extremum —at the frequency w = w, = 6 
(Fig. 1, a). The amplitude-frequency characteristic has 
the same form in the case when N>0 and R>0 but R*<3N. 

3. Frequencies wy, and w, are both real only if the 
conditions N> 0,R> 0, andR*> 3N all hold. With 
this, the amplitude-frequency characteristic has the form 
of a curve with extrema at the frequencies w, , w,,and 
ws (Fig. 1, c). 

The conclusions just drawn, together with (7), exhaust 
the problem of investigating the frequency properties of 
third-order systems. 

For the practical use of these results, it is very 
convenient to employ a diagram with coordinates N and 
R (Fig. 2), on which the various hatchings indicate the 
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MAH. 
Fig. 2. I is the region of one extremum, 


II is the region of two extrema, and III is 
the region of three extrema. 


Fig. 3. I is the region of one 
extremum, II is the region of 
two extrema and III is the re- 
gion of three extrema. 


regions corresponding to systems with various types of 
amplitude-frequency characteristics. Similar to those, 
like the well-known Vyshnegradskii diagram [1], which 
permit one, without solving equations, to judge of the 
stability and character of the transient response of a 
third-order system, the diagram of Fig, 2 gives a represen- 
tation of its amplitude-frequency characteristic. 

It is interesting to note that the juxtaposition of the 
diagram proposed by us with the Vyshnegradskii diagram 
(Fig. 3) graphically illustrates the relationship between the 
character of a third-order system's transient response and 
the form of its arnplitude~frequency characteristic. 

This relationship is implicitly contained in the 
relationships obtained by Vyshnegradskii and in the re- 
lationships given in the present work and, if necessary, 
could easily be elucidated analytically. 
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DETERMINATION OF THE INTERVAL OF POSITIVENESS 
OF THE REAL FREQUENCY CHARACTERISTIC FROM THE GRAPH 


OF THE TRANSIENT RESPONSE 


L. P. Shiniberov 
(Leningrad) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 420-421, March, 1960 


Original article submitted April 20, 1959 


If one knows the transient response of a closed auto- 
matic control system to a unit step stimulus, then one 
may find from it, by a simple graphical construction , 
the boundary frequency of the positive interval of the 
real frequency characteristic w,, without a detailed 
determination of the entire characteristic. 
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Fig. 1. Determination of Wp from 
the graph of the transient response. 


We assume that the curve of Fig. 1 is the transient 
response referred to for the output coordinate a, = f(t). 
On the figure we draw two auxiliary lines; line DE 
parallel to the axis of abscissas from the axis of ordinates 
to the maximum point of the transient response, and line 
AB, coinciding with the majority of points on the ascend- 
ing (increasing) part of the transient response curve, from 
the axis of abscissas to line DE. We bisect line AB and 
find the value of time t,, which corresponds to the mid- 
point of this line (point C). Then, the boundary fre- 
quency of the positive interval can be determined from 
the formula 


=. 1 

“=a (1) 
In the example given, Wp = n/ 2(0.5) = 3.14 sec™?, 

which is close to the true value of Wp = 3.20 sec™4, 
We now prove the validity of the result obtained. 
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Fig. 2. Typical trapezoidal 
transient response. 


We divide the curve of the transient response, given 
in Fig. 1, into three parts; 1) the unit step transient 
response +1(t)i 2) the positive trapezoid DEFL lying 
above the function +1(t) with altitude Ky = DL: 3) the 
negative trapezoid DBAO with altitude K, = -OD. From 
these trapezoidal characteristics it is possible to find 
the real frequency characteristic of the automatic control 
system. If we denote the elements of the trapezoidal 
characteristic as was done in Fig. 2, then in the given 
case the system's real frequency characteristic R(w) will 


be expressed by the formula [1] 2 


. cos tw sin A.@ 
Rio) = 3 | K,— okt (“2 ) (")] i¢ 
t=1 

For the value w = w,, this expression must vanish. 
Since the positive trapezoid DBEFL is always several 
times wider than the negative trapezoid DBAO, the 
maximum of the frequency characteristic constructed 
from it will lie in the region of low frequencies, and it 
will affect the form of the total frequency characteristic 
only for frequencies less than Wp while for frequencies 
close to wp or larger than Wp» it will have a constant 
value, equal to Ky. 

This means mathematically that, for frequencies 
w> 5/t, the term wKyty{cos t, w/t, w)(sinA,w/Ayw) 
in (2) can be ignored, since it is significantly smaller 
than 0.1K,;. Thus, in the neighborhood of frequency 
W,, the real frequency characteristic, when account is 
taken of the fact that K; + K, + 1 = 0, will be approx- 
imately expressed by the formula 











cos tow 
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R (@) = WK oto ( 
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The curve for R(w) constructed in accordance with 
(3) will pass through zero for the first time for cos %w/t, w= 
= 0, which occurs for the value tw =" /2, Le., for 


t:=lm 


tow 


o=-, since 
2tm 

As is clear, the frequency w = m/ 2t,, is the boundary 
frequency of the positive interval, f.e., (1) 4s valid. 

The value of w, thus obtained may serve as an 
estimate of the cutoff frequency w,, of the system's 
logarithmic amplitude characteristic since, as is well 
known, the logarithmic amplitude characteristic intersects 
the frequency axis for the values [2] 


@, = (0.55 — 0.7) a,. (4) 


The suggested method of determining wp gives good 
results when one is considering transient responses, either 
determined experimentally or stated as desirable, of 
power servo systems which are distinguished by compar- 
atively small oscillatory indices; M= 1.1 to 1.4. 

For other transfer functions, the error reaches 15 to 
20 %, and this method can be used only for an approximate 
estimate of Wy 
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PROCESSING EXPERIMENTAL FREQUENCY CHARACTERISTICS 


G. I, Monastyrshin 
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Translated from Avtomatika i Telemekhanika, Vol. 21, No. 3, pp. 422-428, March, 1960 


Original article submitted May 20, 1959 


For recording forced oscillations, loop oscillographs 
are ordinarily used, recording both the stimulus and the 
reaction of the element. In those cases where the oscil- 
lations do not hve a perfectly sinusoidal form, oscil- 
logram processing is carried out by means of a harmonic 
analyzer which separates out the harmonic component of 
the reaction which has the same frequency as the stimulus; 


for brevity, we shall henceforth cali it simply the reaction’ 


harmonic component. 

Figure 1 gives a set of oscillograms of a servo drive's 
oscillations, established under the action of sinusoidal 
stimuli of various frequencies and constant amplitudes. 
On these oscillograms, U is the stimulus and V is the 
servo drive's reaction. The oscillograms also bear the 
reaction's harmonic component. The cross on each 
oscillogram marks the maximum deviation of the servo 
drive's oscillation from the given harmonic component. 


s 


We note that, in the example given, not only the sinusoidal 


character of the oscillation is disrupted, but, as is clear 
from Fig. 1, its periodicity is also disturbed to a certain 
extent, so that in the cases of Fig. 1,c andd the 
harmonic analysis {fs carried out over an interval of 
several (two) periods of the stimulus. 





The results of the measurements are collected tn the 
table, where wy is the angular frequency of the harmonic 
stimulus, A; {s the stimulus amplitude, B, {s the amplitude 
of the reaction's harmonic component, ¢; 1s the phase 
shift between the stimulus and the reaction'’s harmonic 
component and AB; fs the maximum deviation of the 
reaction from its harmonic component. 

We remark that the table can be obtained without 
the intermediary of an oscillograph by means of direct 
measurements using the method cited in [1]. 

In what follows we shall not deal with any concrete 
form of the deviations of the forced oscillations from 
their harmonic components, but shall deal with the study 
abstractly, assuming that any character of these deviation 
is possible with limitations on the modulus of the AB; . 
Corresponding to such a treatment, the forced oscillation 
are characterized, not by lines, but by the bands included 
between two sinusoids, equidistant from the harmonic 
component of the oscillation,and separated from it by the 
maximum deviation marked on the oscillograms. 

By considering the reaction's harmonic component 
as the projection of the polar vector By , we can interpret 
the deviations of actual oscillations from their harmonic 























i 20sec -—| Aj, mm| Bj-mm| % degreg ss, mm 
1 3.123 31.0 20.95 — 47°20’ 3 
2 6. 148 31.0 13, 83 — 74°03’ 4 
3 12.45 31.0 7.36 —101°46’ 3 
4 25 13 31.0 3. 61 — 138°08’ 5) 














components as the projections of an auxiliary vector ABj, 
limited in modulus. 

If we abstract from the concrete definiteness of the 
rotation of this auxiliary vector we obtain, instead of the 
end of the vector B,, a circular smear with radius AB, . 

The vector interpretation allows one to note that 
when ABj« B, (for our purposes, it suffices in practice if 
AB, < 0.5 By) the ratio AB, : B, is determined approxima- 
tely both by the maximum relative variation of the 
harmonic component's amplitude and by the maximum 
variation of its phase (measured in radians) for which the 
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harmonic oscillations will remain within the limits of 
the smear defined earlier or, what amounts to the same 
thing, within the limits of the forced oscillations’ bands 
In this sense, instead of the lines of the frequency, 
amplitude and phase characteristics, one can consider 
the bands of the amplitude and phase characteristics. 
Figure 2 shows such bands, constructed from the oscil- 
lograms of Fig. 1. In the neighborhood of the frequency 
4 cps, the bands of the phase characteristic are not bowie 
due to the fact that the condition on the smallness of 
AB, does not hold. 
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Fig. 1. 


In using an element's experimental frequency 
characteristic, it may turn out to be advantageous to 
define the system by the differential equation correspond- 
ing to these characteristics [2-6]. For example, the 
problem arises in mathematical modelling of setting up 
the scheme which reproduces the frequency characteristics. 

We assume that the linear differential equation 
defined by the experimental frequency characteristics 
has the form 

n 


™m 
> 2 7%= > oy, (1) 
k=0 k=0 
where the U(K) are the stimulus and its derivatives, and 
the v(k) are the reaction and its derivatives. 


As represented on the oscillograms, the stimulus U; 
and the reaction’s harmonic component V, may be 
written in the form 


U,=A,sin wf, (2) 
V, = B, sin (wt + 9). (3) 


To determine the numerical values of the coefficient 
a, and by of (1), we use the approximate equation ob- 
tained by substituting (2) and (3) in (1): 
n 


> a,A,o* sin (w + >t) = 


k=0 


2 he wtsin(oe +t FA) (4) 


To provide a graphic representation, we denote the 
left and right members of (4) by polar vectors Ry and S; 
(Fig. 3). 
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To find the unknown coefficients aj, and by, we start 
from the condition 


min ri ( | ): 


|”, | 





(5) 


where €; = $1-Ry and the py are weight coefficients. 

It is easily seen from Fig. 3 that if, as a result of 
seeking the values of the coefficients of (1), the vectors 
€, Me within the limits of certain circular smears whose 
radii are in the same ratios to the moduli of vectors Ry 
as the ABy are to the By, then the forced oscillations 


defined by (1) , with the values found for the coefficients, 


lie within the limits of the bands considered earlier. 


In correspondence with this, by setting the requirements 


on the accuracy of the approximation of the frequency 
characteristics, and omitting the arbitrary coefficient of 
proportionality, we establish the weight coefficients: 


a (a) (8) 


If in (6), instead of the quantities AB; , which depend 


on the frequency of the forced oscillations,we use the 
quantity max AB, and omit the constant factor thereby 
entailed, then 


7 
P, = Bi. © 
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With this, if the stimulus amplitude is taken to be 
constant, then the width of the corresponding band of 
the amplitude characteristic ceases to depend on fre- 
quency, and the width of the band of the phase charac- 


teristic turns out to be inversely proportional to the value 
of the amplitude characteristic (for sufficiently small 4B) 
It is obvious that the bands corresponding to (6) lie inside 
the bands defined by (7). 


We note that (5) for | g |<<| R;,|, as may be observed 
from Fig. 3, is equivalent to the condition 


i 


where the A; are the absolute deviations of the phase 
characteristic of (1) from the experimental one, measured 
in radians, and the 6, are the relative deviations of the 
amplitude characteristic of (1) from the experimental. 

For what follows, it is advantageous to give (5) in 
the form 


“Dal(sar) +(e) | @ 


where xj and y; are the projections of vector €;. 
It is easily seen that 
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k=0 


With (10), (11) , and (12) kept in mind, we shall have the following equations for determining the coefficients 


sought ; 








Since the quantities following the partial differentia- 
tion signs in (13) do not, in the final analysis, depend on 
t, one can set t = 0 in (10), (11), and (12) as a first step, 
thus simplifying the computations. 

If (1) does not contain derivatives of U or V then (13) 
will be linear since, in the first case, |Ryl = Ay (we set a 
equal to unity) and, in the second, |R,| a) 1 S;| = By (we 
set by equal to unity). 

Otherwise, the coefficients being sought can, to a 
first approximation, be determined, for example, from 


the condition S| 
min )) |, |. 
i 


To render more exact the solution thus obtained one 
might determine the quantities | Rj| graphically from the 
vector diagram or by (12) for t = 0 and the numerical 
values of the coefficients of the first approximation, 
where for the differentiation of the fractions in (13), one 
differentiates only the numerator (as in the case when 
the denominator does not depend on the parameter with 
tespect to which the differentiation is being performed), 
which is possible if the numerator is sufficiently small 
with respect to the denominator. 

Omitting the detailed computations, we give the 


final expressions determining the coefficients of the 
equation 


(14) 


U aia a,U + aU + a,U ®) a a,u™ a a,U) — 


= bV +b,V + bY + dV + dv + by, (15) 


In terms of this example, we consider the structure 
of the formulation which permits one to obtain the analo- 
gous formulas for any more complicated case. The case 





(13) 








» minthecase ap = :) 


e+, ™ inthe case bp = 1 


considered apparently encompasses all practically 
possible particular cases. All these particular cases may 
easily be obtained from the one considered by simply 
omitting the proper terms of the expressions provided. 

The normal equations defining the coefficients of 
(15) are given by (16). 

In these equations,p; kj? = 1 (it is also possible to 
set k= 1) for the first approximation in correspondence 
with condirion (14), and may thereafter be determined 
graphically in the rendering more accurate of the solution 
by the coefficients of the first approximation by means 
of the vector diagram (Iy” = |Ry| “*) or by 


k? =A~? [(a,@, — agw? + agw?)? 


a7) 
+ (1—aqu? + agot)3}-1, 


For simpler cases it suffices, in (16), to delete the 
columns, and the rows symmetric with them about the 
principal diagonal in the matrix of the unknown coef- 
ficients, for the coefficients of the higher derivatives 
which are lacking in (15). The corresponding terms in 
(17) are also omitted. In the particular case when there 
are no derivatives in the left member of (15), there 
remain only the last six columns of the left side of (16) 
and, correspondingly, only the six last equations k;"=A4~*, 
and the problem {s significantly simplified. If there are 
no derivatives in the right member of (15), it is then 
advantageous to change the problem slightly, setting 
bp = 1 and taking a» as an unknown. In this case, the 
simplified normal equations obtained from (16) are 
analogous to the preceding ones if one first makes the 
following changes in (16) ; the by, are replaced by the 
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Fig. 4. 


a, and the a, are replaced by the by and, in the ex- 
pressions for the coefficients cj, dj, ey, and f;, Ay fs 
replaced by By B, is replaced by Ay % is replaced by 
-g; and the ky are defined as By~2. One should then 
retain the six last equations and, correspondingly, the 
six last columns of the left members of these equations, 
determining the p, by (6) or (7). 

If, as a result of the computation of the coefficients 
of (15), the difference between its frequency characteristics 
and the experimental ones, defined by the harmonic 
components of the actual oscillations, turns out to be 
sufficiently small-for example, significantly less than 
the possible deviations inside the previously defined 
bands of the forced oscillations -we shall then consider 
that the problem of determining the linear equation 
which approximates the experimental frequency charac- 
teristics has been successfully solved. 

In the concrete example considered, given by Figs. 
land 2 and by the table, we are investigating an equation 
of the form 

U+ aU = beV + dV, 
for which we shall have the normal equations 


(18) 


a >)P 
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2 42,2 
iki Ajo; — 


a ° P 
— bo ps pki A, Bo; sing, — by >) Pik24 Bw? cos 9, = 0. 


i i 


. — ‘ 
= a; >) PRA, Bio, sin 9; + bo > p,Be= 
i . 


: 


= a pki A,B, cos ¢;, 
i 


—a 9) pjk{A,Bywi cos 9; + (19) 
i 
+b, p : p,ki Be i=—)) PATA, B,wo,sing, 
i i 
Setting the p, by (6) and taking, in correspondence 
with (17), 
k? = A~® (a2w? + 4)-1 (20) 
we obtain , as the result of computing the first and 
second approximations and the corresponding rounding, 


a, = — 0,034, bo = 1,11, 6, == 0,325. 
The corresponding frequency characteristic is given 
in Fig. 2. The circles on the figure mark the experimen- 


tal points, while the heavy lines show the results of the 
approximation . 
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It should not be overlooked that the frequency 
characteristics considered can not determine all the 
properties of the actual element. Due to this, the 
question as to the use of sucha type of characteristics, 
and the equations approximating them, can form the 
subject of a concrete study. 


For the example considered, Fig. 4 shows the 
reaction to a unit step function stimulus, obtained from 
(18) by computing the values of the coefficients, and 
the analogous reaction of the actual servo drive. 

It may happen that, in practice, one does not 
succeed in obtaining a sufficiently small error of 
approximation of the experimental frequency characteris- 


tics. For example, these errors may remain beyond the 
limits of the experimentally defined bands of the forced 
oscillations. In this case, the question as to the possibility 
of using the results of the approximation can be decided 
by fixing the band of the forced oscillations in cor- 
respondence with the approximation errors obtained. 


We note that a treatment analogous to that presented 
here is also possible in the case when the experimental 
results are free from nonlinear deviations, and only 
measurement errors distort them. 
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INFORMATION . 


THE EMU-8 ELECTRONIC ANALOG COMPU TER 


As the result of the work on improving electronic 
analog computers and their elements which has been 
carried on at the Institute for Automation and Remote 
Control of the Academy of Science of the USSR in con- 
junction with a collective under the direction of G. N. 
Grigor’yan, there has been created a new electronic 
analog computer, type EMU-8, production of which has 
already begun. 


Domain of Application of Type EMU-8 
Electronic Analog Computer 

The type EMU-8 electronic analog computer is 
designed for the investigation of linear and nonlinear 
dynamic systems. Structurally, it is a general-purpose 
computer capable of real-time operation (joined with 
elements of actual apparatus) with accelerated and 
retarded solution processes within the limits of its passband. 
The computer can be used to greatest effect in the 
solution of technological problems where the error in 
giving the original parameters is sufficiently large (from 
5 to 10%). 








Among these problems may be listed: 

1) investigation of the dynamics of automatic 
control and regulation systems (automated systems, 
missile control, boiler aggregate regulation, steam 
turbine control, various automated electric drives, etc.); 

2) certain problems of nuclear physics and technology; 

3) investigation of combustion processes in internal- 
combustion engines; 

4) investigation of transient responses in complicated 
electrical circuits; 

5) investigations in the area of nonlinear mechanics. 


Brief Technical Characteristics of the 





EMU-8* 

The EMU-8 computer consists of individual building 
blocks each of dimensions 450 x 460 x 320 mm. Each 
building block can provide the solution of linear and nor 
linear differential equations up to second-order, inclusive, 
* For more detailed technical characteristics of the 
computer plus its functional schematic , see Vesmik 
Akademii Nauk No. 7 (1958) . 
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while the standard package of building blocks provides 
solutions of differential equations up to tenth-order, 
inclusive. Plug-in units permit the execution of the fol- 
lowing mathematical operations; multiplication and 
division of two independent variables, multiplication by 

a constant factor, differentiation, integration, reproduc- 
tion of various nonlinear functions and standard nonlinear- 
ities. 

The character of the mathematical operations exe- 
cuted {s determined by the type of plug-in units, which 
allows operational amplifiers to be used economically. 

By choosing several building blocks with the neces- 
sary variants of plug-in units, one can solve problems of 
various complexity. 


In the standard package, five building blocks are 
used. 


Six block indicators are connected to each package 
to provide visual observation of problem solution. 


Each building block is designed in the form of a 
completely independent instrument, and may be con- 
nected directly to a 220v line. The power requirement 
of each block is about 160w. 


Line production of the computer was organized in 
1960. The wholesale price of a package is 135,000 
rubles. 


Bureau of Factory Information 





